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ABSTRACT 


Volume 1 of this report discusses the radiative and collisional 
processes which determine the recombination rate in a plasma. A gen¬ 
eral expression is derived for the recombination rate, and numerical 
values are given for a hydrogen plasma in the temperature range 250°K 
to 64,000°K and for all electron densities. These results are applied to 
an expanding hydrogen plasma and to the expanding debris from a nuclear 
explosion, and indicate that in an explosion occurring in a vacuum, sig¬ 
nificant recombination may take place at times later than a millisecond 
after the detonation. 

Volume II contains seven papers dealing with the physics of recombi¬ 
nation in plasmas. Included are the radiative transition probabilities for 
a number of atoms and ions obtained from a theory of atomic screening 
constants. 
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THE STATIONARY PLASMA* 


2.1 The Basic Equations 

In order to investigate the non-equilibrium properties of a 
partially ionized plasma, it is necessary to consider the electronic and 
atomic processes which affect the populations of the individual atomic 
energy levels. A particular level may be populated by spontaneous and 
stimulated transitions from other discrete levels, by impact induced 
transitions from bound states, by radiative recombination of free elec¬ 
trons and by three-body recombination processes involving free electrons. 
It may be depopulated by the inverse processes of radiative transitions 
out of the level, impact induced transitions out of the level, photo¬ 
ionization and impact induced ionization. 

The increase in population per second of a particular level (p) is 
given by the differences between the rates of population and depopula¬ 
tion. The equations for the populations in each of the levels will 
have the form 

dt n(p) = /, a jp n(J) + b p • a PP n(p) (2_1) 

J Vp 

In this equation n(p) is the population of the p tb level. The expres¬ 
sion a^ n(j) appearing in the summation represents the rate at which 
atoms are entering the p tb level from the J tb level. The terra b^ 
represents the rate of capture of free electrons into the p tb level 
and the expression a^ n(p) represents the rate at which atoms arc 

*This chapter is repeated in Volume II to make it available in an un¬ 
classified document. 
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leaving the p th level for other levels and the ionized state. The a's 
and b's include the effects of both radiative and collision processes, 
and they are functions of the density and temperature of the free elec¬ 
trons. Explicit expressions for the coefficients are developed in 
Appendix A. 

2.2 Comparison of the Radiative and Colllslonal Processes 

Tables giving the coefficients for the radiative and colllslonal 
processes are collected in Appendix A. The tables for hydrogen are 
essentially complete. For other elements the tables represent a start 
on the complete assembly of data needed for estimating recombination in 
complex plasmas. For the present they serve as an indication of the 
compilation and are useful in showing the pattern of variation of the 
coefficients. 

An examination of the tables shows that the spontaneous transition 
probabilities decrease rapidly with increasing excitation, in contrast 
to the colllslonal excitation and ionization probabilities, which in¬ 
crease rapidly. Thus, from the tables for hydrogen in Appendix A, it 
may be observed that even if the electron density is as low as 10*^ cm ^ 
and the electron temperature is as low as 1000 U K, electron impact ioniza¬ 
tion of level n = 10 is more rapid than the loss by spontaneous transi- 
20 -3 

tion. For an electron density of 10 cm and an electron temperature 
of 1000 a K, electron impact ionization of the first excited level is more 
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rapid than the loss by a spontaneous transition. Electronic excita¬ 
tion and de-excitation processes are still more efficient so that a 
Boltzmann distribution of the highly excited states is to be expected 
except at low electron densities and temperatures. 

2.3 Previous Calculations of Recombination Rates 
( 2 - 2 ) 

Baker and Menzel have applied equations of the form of 

Equation (2-1) to an investigation of the relative intensities of 
spectral lines in gaseous nebulae in which the ground state of atomic 
hydrogen is ionized by stellar radiation. In these circumstances, a 
steady state prevails so that 

57 n(p) = 0 

and Equation (2-1) reduces to a set of algebraic equations, the solu¬ 
tion of which is equivalent to the Inversion of a matrix. Baker and 
Menzel do not distinguish between the individual sub-levels corres¬ 
ponding to a given principal quantum number and consequently assume 
implicity that 

n(n,l) = \ n(n,i) (2-2) 

n V 
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/ 2 - 3 ) 

Burgess' ' has repeated the calculations taking account of the 
different sub-levels and he gives the equilibrium populations for electron 
temperatures of 10,000°K and 20,000°K. He finds that the ^-levels are 
far from being populated according to their statistical weights and that 
for an optically thin plasma, the p-levels (4«1) are much under-populated 
due to the large probability of np - Is transitions. 

In both investigations, all collision processes are ignored. Al¬ 
though this is a permissible approximation for the low-lying levels of 
hydrogen when the electron density and the electron temperature are not 
high (as is the case in most gaseous nebulae), it will ultimately be 
invalid as n increases. 

Since collision processes which redistribute the sub-levels corres- 


ponding to a given n, such as 


e + li(2p) - e + H(2s) 

(2-3) 

H + + H(2p) - H' ! ' + H(2s) 

(2-4) 


are very efficient, the effect of including collision processes will 
tend to ensure Equation (2-2), The more refined treatment of Burgess is 
therefore actually less accurate than that of Baker and Menzel when the 
electron temperature is high and when the principal quantum number n 
is large. 

Giovanelli^"^ has included electron collision processes in his 
treatment, which again assumes statistical equilibrium. In his first 
paper in which he assumes Equation (2-2jjio shows that Saha's ionization 
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formula applies to the populations of the high-lying levels, provided the 
electron temperature is employed. This occurs because the most efficient 
depopulating mechanism is electron impact ionization and the most effi¬ 
cient populating mechanism is the inverse process of three-body recom¬ 
bination and these processes are very rapid. The same situation must 
apply even in non-equilibrium cases. 

(2-5) 

In a second paper Giovanelli discusses the populations of the 

low-lying levels. The study contains an inconsistency, which may be 
serious, in that while he distinguishes between the individual sub- 
levels, he ignores redistribution processes such as Equation (2-3). 

His results are therefore invalid at high electron densities and tempera¬ 
tures. We note in passing that his results are also Invalid at low elec¬ 
tron densities and temperatures in that he ignores the decay of the 
metastable 2s state by the two-quantum process 

H(2s) - H(ls) + hv t f hv 2 (2-5) 

Discussions of the hydrogen plasma in non-equilibrium situations 

(2-6) (2-7) 

have been given by D'Angelo' 1 and by Bates and Kingston' . In 

both investigations, heavy particle collisions are ignored and it is 

assumed that the plasma is optically thin. Both investigations assume 

Equation (2-2) and ignore processes such as Equation (2-3). This is 

equivalent to supposing the Equation (2-2) is very fast. 

D'Angelo obtains an approximate solution of Equation (2-1) by attemp¬ 
ting to follow through in detail the reaction paths of any individual 
electron, making the assumption that recombination is effective only for 
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Chose paths which terminate in the ground state. He ignores the effect 
of electron impact excitation and de-excitation. He gives results for 

three electron temperatures, 1000°K, 3000°K and 10,000°K and for electron 
12 13 -3 

densities in the range from 10 to 10 cm only. His investigation 
reveals the importance of three-body recombination processes (as did the 
earlier work of Giovanelli) but his method of solution is laborious to 
apply and the accuracy of the numerical values of the effective recombina¬ 
tion coefficients is uncertain. 

The method employed by Bates and Kingston is simple and more accurate. 
They argue that the rates at which atoms in a particular excited state 
are produced or destroyed are very much greater than the rate at which 
the very small number density of atoms in that state changes. It follows 
that for an optically thin plasma except for the ground state 

n( P ) =0 (P = 2, 3, ...) (2-6) 

so that recombination is effective only into the ground state. Then 
writing 

3? "(» ' + We) ". 2 < 2 ' 7 > 

(where n(l) signifies the ground state) 

defines the effective recombination coefficient “gf f( n e * T e )» which will 
depend upon the electron density as well as the electron temperature. 

Their method reduces the set of coupled differential equations, Equation 
(2-1), to a set of coupled algebraic equations. The number of levels re¬ 
quiring individual examination is restricted since the populations 
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of the highly excited levels are given by the Saha formula, the level 
at which this becomes true being provided by the calculations. 

2.4 The Number of Discrete Levels 

In none of the investigations we have described has the question of 
the number of levels that should be included in the equations been exam¬ 
ined. It is clearly not infinite because ultimately the discrete levels 
will overlap the continuum. 

Ivanov-Kholodnyi, Nikol'skii and Gulyaev' ' have discussed the 
lowering of the continuum due to pre-ionization. The treatment is not 
entirely satisfactory but should suffice for our purposes. They find 
that the number n Q of levels of different principal quantum number which 
remain discrete is given by 

21.65 - £n(n ) 

fn(n Q +1)- ^ -— (2-8) 

Values of n Q have been computed from Equation (2-8) for a wide range 
of electron densities and they are presented in Figure 2-1. The values 
are quite large and it seems probable that Stark broadening of the lines 
will cause a greater effective depression of the levels. The Stark 
broadening has been discussed by GrJem.Kolb and Shen^ 2 " 8 ^ and it may 
be necessary to carry out detailed computations based on their theory. 
However, the fact that the highly-excited levels are in essentially 
thermal equilibrium and play only a small role in determining the effect¬ 
ive recombination suggests that the solution of the equations is insensitive 
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(°u) 31 VIS QNnoa 1S3H9IH 

Figure 2-1. Highest Bound State of an Electron j.n Hydrogen 
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to the number of levels included. The Stark broadening will be impor¬ 
tant, however, in determining the opacity, which in turn effects the 
recombination. 

2.5 Computation of Recombination Rate in Optically Thin Hydrogen 

(2-7) 

Adopting the method of Bates and Kingston , we first change 
the variables in Equation (2-1) from n(p) to 0(p) by introducing the 
relationship 

n(p) - n fi n(X + ) o> <h 2 /2nmkT e ) 3/2 exp(I p /kT e ) 6(p) (2-9) 

(where • is mass of electron, k is Boltzmann's constant, h Planck's 
constant and l p the ionization potential of this p th level.) 

The new variables 0(p) are the ratios of populations in the p th levels 
to the populations derived from the Saha relationship. At equilibrium 
0(p) is unity. 


The introduction of (2-9) into Equation (2-1) gives an equation of 
the form* 


ii 8(|>) -X ‘ip n(j) 

Up 


b 

P 


a 

PP 


0(P> 


P = 1, 2, 3, ... 


(2-10) 


Assuming with Bates and Kingston that 

d|M = 0 p » 2, 3, 4,... (2-11) 


*the a's do not have same meanings in Equation (2-10) as in Equation 

( 2 - 1 ) 
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we obtain 


0 a JP n(J> + b p ’ a pP 9<P) 

j*P 

p = 2, 3, 4,... (2-12) 

with the equation for the ground state where p = 1 treated as a special 
case: 


^ ®(1) = y aj! 0(J) + b x - a u 0(1) (2-13) 

jVp 

If the series of equations in (2-12) are truncated at some level k, the 
Equations (2-12) are in the form of k-1 linear equations with k unknowns, 
the unknowns being taken aB the 0's. By inversion of the matrix, a 
solution may be obtained for all d's in terms of one unknown, 3(1). Be¬ 
cause of the linearity of the equations, the solution will have the form 

0(p) = c q + c L 0(1) 

p «= 2, 3... k (2-14) 

The right side of Equation (2-13) can then be evaluated also in terms of 
the single unknown 0(1) and taking account of Equations (2-9) and (2-7), 
one My obtain the effective recombination coefficient in the form: 

a ,f£ ( "JV " *0 ' "l •<»> «' l5 > 

Values of a Q and a^ in Equation (2-15) have been obtained for an 
optically thin hydrogen plasM by machine methods for electron tempera¬ 
tures ranging from 250°K to 64,000°K and for all electron densities. 
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11 


indices give the power of 10 by which the entries must be multiplied. 












M i 
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Indices give the power of 10 by which the entries must be multiplied. 



levels up to n = 10 being Included in the computation. Preliminary 

values of a and a. hare been provided by Bates and are given in Tables 
o i 

2-la and 2-lb, respectively. 

There are certain irregularities in the values as given in the 
tables which suggest there might be an error in parts of the tabulation. 
Until a further check of the computations has been completed the values 
in the table should be regarded as uncertain by a factor of 2. 

2.6 Dtacjtf3ij.pn.of Compilation flP j Rg.gu.Ks 

In order to illustrate the significance of the differential 
Equations (2-10), the algebraic Equations (2-12), and the solution as 
given in Tables 2-la and 2-lb, it will be worthwhile to consider the 
changes which take place in the populations of the various excited levels 
of a hydrogen plasma. Suppose that for an optically thin hydrogen plasma, 
the temperature and number density of the free electrons are maintained 
constant by some arbitrary means. Under these conditions, the rate of 
change of the populations of all the bound states are given by the 
differential Equations (2-1), and for given initial conditions which 
assign a value to the population of each of the bound states at the out¬ 
set, a solution of the Equations (2-1) will give the populations of each 
of the states for all subsequent time. The general character of this 
solution is illustrated schematically in Figure 2-2. 

In the figure, the populations of the several states are described by 
the variable 0 , which is the ratio of the number in a state to the 
number in that state under conditions of thermodynamic equilibrium. 

Since we are supposing that the number density and temperature 
of the free electrons remains constant, 
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RATIO OF POPULATION TO EQUILIBRIUM POPULATION 


010 * 40 - 1310 



TIME -- 


Figure 2-2. Scheaatlc Dlagraa of Population Changes During 
Recoablnatlon 
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0(p) will be proportional to n(p). In Figure 2-2, it is supposed that 
the ground state is unpopulated at the outset, but the initial popula¬ 
tions of the various excited states are taken at values which have no 
particular significance. As indicated in the figure, at the outset there 
will be a very rapid change in the populations of the upper states. In 
particular, the value of 6 for the upper states will change much more 
rapidly than that of the ground state. After a short time, the rapid 
changes in 8 of the upper states will decline, and a pseudo-steady state 
will occur such that, if viewed on a short time scale, all the populations 
remain constant. There will, however, be a non-zero rate of Increase in 
the ground state, which will make itself felt over a longer period of 
time. This pseudo-steady state condition is representated in Figure 2-2 
as occurring just before the first time break, and it is this condition 
which is described by the algebraic Equations, 2-12 and by their solu¬ 
tion obtained taking the recombination coefficient cr * a Q . At an appre¬ 
ciably later period of time, as represented after the first time break 
in Figure 2-2, the increase in 8 for the ground state will be appreciable, 
and there will be a corresponding change in the pseudo-steady state values 
of the populations of the various higher states. This situation is re¬ 
presented in the middle time-period of Figure 2-2, and is also described 
by the solution to the algebraic Equations, 2-12, but now since 0(1) is 
not equal to zero, it is necessary to take the effective recombination 
rate, as or - a Q - 3^(1). 

The populations of the various states continue to change as recomb¬ 
ination proceeds, and the right hand section of Figure 2-2 indicates 
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the general nature of these changes on a much different time scale than 
the portions of the figures at the left and the center. The populations 
will eventually approach asymptotic values at which they are constant, 
and at which time the net recombination ceases. This condition, which 
is indicated at the extreme right of the figure represents the steady 
state under conditions of radiation loss. It can be seen that for this 
condition, the population of the ground state exceeds that which it would 
have for thermQdynamic equilibrium. The population of the ground state 
can be obtained readily from the Tables 2-la and 2-lb by setting a - 
a Q - a x 0(1) - 0. 

In simplifying the differential Equations 2-10, to the corresponding 
algebraic Equations,. 2-12, one gives up any possibility of obtaining 
information about the initial variation shown at the extreme left of 
Figure 2-2. Except for this initial response, which is unimportant for 
the applications considered here, the effective recombination rate 
and the populations of all the states can be obtained throughout the 
period illustrated in the second and third time blocks in the Figure 
from the algebraic equations. 

It should be observed that in the proceedure followed here, one 
supposes that the temperature and number density of the free electrons 
and the conditions on light absorption as described by the X"Parameters 
are given conditions of the problem. To do otherwise would make the 
equations (either the differential Equations 2-10 or Lhe algebraic 
Equations 2-12) non-linear and would vastly increase the difficulties of 
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obtaining a solution. The solution obtained here can nonetheless be 
applied to important problems in which the temperature, density and 
light absorption in the plasma change with time because the pseudo¬ 
steady state to which Equation (2-12) applies is achieved in times 
short compared to those for significant variation in the plasma properties. 
This matter will be discussed further in Section 3 of Volume I. 

2.7 Plajma, Otfhqr. Than .jfrdyofflin 

The general plasma has several ionized components each of which 
can exist in several states of ionization. Nevertheless, the general 
plasma may be treated in a manner similar to that employed for the hy¬ 
drogen plasma, but with certain modifications. 

Just as for the hydrogen plasma, it is necessary to collect to¬ 
gether all the radiative and colllsional transition probabilities. Un¬ 
fortunately, no systematic tabulation of these probabilities has been 
given even for- ground states of systems and the collection of the data 
is a major task. We have made a beginning in assembling the data and 
the tabulated value are included in Appendix A. 

The basic equations are similar to (2-1) and their solution pre¬ 
sents no problem provided heavy particle collisions can be ignored. 

With this proviso, there is a set of equations of the form (2-1) for 
each component of the plasma and these sets of equations are uncoupled. 

The quasi-steady state method may then be applied to each component 
separately. 
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It does, in fact, appear reasonable to neglect heavy particle 
collisions except for charge transfer processes Involving the ground 
states. This exception arises because, although charge transfer into 
the ground state of a system may be inefficient, it does couple compo¬ 
nents together and it has a permanent recombining effect. It does not 
introduce any difficulties into solving the equations since the ground 
state equation receives special consideration in any event. 


18 



REFERENCES 


2-1 Ivanov-Kholodnyi, G. S., et al . "Ionization and Excitation of 

Hydrogen I. Elementary Processes for the Upper Levels", Soviet 
Astronomy-AJ 4 , 754 (1961). 

2-2 Baker, J. G. and Menzel, D. H., Astrophys. J. 86, No. 1 (July, 
1937). 

2-3 Burgess, A., "The Hydrogen Recombination Spectrum", Mon. Nat. 
Roy. Astron. Soc. 118 . 478 (1958). 

2-4 Giovanelli, R. G., "Hydrogen Atmospheres in the Absence of 
Thermodynamic Equilibrium" 

I. "The Population of the Upper Atomic States, Australian 
J. Sci. Research Al, 275 (1948) 

2-5 Giovanelli, R. G., "The Population of the Lower Atomic States" 
Australin J. Sci. Research Al, 289 (1948). 

2-6 D'Angelo, N., "Recombination of Ions and Electrons", Phys. Rev. 
121 . 505 (1961). 

2-7 Bates, D. R. and Kingston. A. E., "Recombination through 
Electron-Electron Collisions", Nature 189 . 652 (1961). 

2-8 Griem, H. R., Kolb, A. C., and Shen, K. Y., "Stark Broadening 
of Hydrogen Lines in a Plasma", Phys. Rev. 116 . 4 (1959). 


19 



APPENDIX A 

A. Dalgarno, A. Naqvi, E. Ashley and H. E. Stubbs 
ATOMIC PROCESSES IN A RECOMBINING PLASMA 

The time dependence of the population of each state of a recombining 
plasma can be obtained from a consideration of the several processes by 
which the atoms enter and leave a given state. Both radiative and 
collisional processes must be considered. 

A.1 Radiative Transitions Between Bound States 

The rate of depopulation of the state p by spontaneous transitions 
to lower states q, may be written 

n(p )£ A(p, q) (A-l) 

q<p 

where A(p,q) is the probability of a spontaneous transition from the 
state p to the state q measured in sec 

If the plasma is not optically thin, stimulated emissions and absorp¬ 
tions occur. For the rate of upward transitions which depopulate the 
state p, we write 

Y_ w r X(p ’ r) A(r ’ p) <A-2) 

P *>p 

where is the statistical weight of the state p and X(p,r) is a factor 
which takes account of the opacity. In particular for a plasma in ther¬ 
modynamic equilibrium with the radiation field 



X(r,p) 


<exp (hv 



(A-3) 


where hv^ the energy of the r - p transition. For an optically thin 
plasma,Xvanishes. The rate of stimulated downward transitions can be 
written 


n (p) 


X(p,q) A(p,q) 


(A-4) 


The p state can also be populated by bound-bound radiative transi¬ 
tions and we have for the rates of population by the three processes 


Z Mr) A(r,p) 

r>p 


Z «<*> 

!>p 


X(r,p) A(r,p) 


(A-5) 


(A-6) 


Z X(q,p) A(p,q). (A-7) 

q<p q 


The probability of a spontaneous transition from a level p to a 
lower level q may be expressed in the form 


,. 4 3 

A(p >q) - —" S(p,q) 
3hc u 

P 


(A-8) 


where v is the frequency of the emitted radiation, u is the statistical 

P 

weight of the p^ level and S(p,q) is the transition strength, defined 
according to 

s<P,q> -« 2 |<¥ p . >| (*-« 
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where ¥ and ¥ are the configurations of the initial and final states 

p q 

and r is the position vector of the j " electron. In the case of Russell- 
Saunders coupling, S(p,q) reduces to the form 

S(p,q) -/(L) /(M) a 2 (A-10) 

where /(L) is the relative strength of a multiplet within a transition 
array, ^(M) is the relative strength of a line within a multiplet and 
a is a transition integral involving radial wave functions only. Tables 
of and (M) are available for most multiplets (Goldberg^ A , 

White and Eliason^ A Russell^ ^, Allen^ A ^, and if necessary, they 
can easily be extended following Rohrlich^ A ^). The transition integral 
o presents greater difficulties. If R(n^j|r) is the initial radial wave 
function of the active electron and R(n'i'|r) is the final wave function, 

2 

(A-11) 

where is the lower of the azimuthal quantum numbers 2 and 2' which 
characterize the active electrons. 

For hydrogen, values of a^(nX,n , J*) have been computed by Green, 

Rush and Chandler^ A ^ and the corresponding values of A(n£, n'jJ 1 ) for 
n and n' less than or equal to 8 are presented in Table A-l. Of interest 
also are the total transition probabilities A(n,n') obtained by summing 
A(n.C, n’^’) over all possible values of 2 and i' consistent with n and 
n', assuming that the sub-levels are distributed according to their 
statistical weights. The total transition probabilities so obtained are 
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presented in Table A-2. There are ocasional differences from previous 
tabulations of Menzel and Pekeris^ 

For more complex systems, approximations must be made. For many 
cases and especially for the excited states of light elements, the 
Coulomb approximation of Bates and Damgaard^* ^ should be adequate. In 
this approximation the active electron is regarded as moving in the 
field of an effective nuclear charge being determined from the observed 
ionization potentials for the states. Bates and Damgaard have presented 
2 

tables from which a may be determined for transitions involving small 
values of & and £ and we have used them to calculate the transition 
probabilities of the more important lines of Li I, Li II, C I, C II, 

C III, N I, N II, N III, 0 I, 0 II and 0 III. The results, illustrative 
of the information which may be obtained by the Bates and Damgaard 
methods, are presented in Table A-3. Also, Table A-4 lists the results 
of calculations by this method of the transition probabilities of 
resonance lines of a number of atoms of interest. 

For heavier elements it may not be possible to apply the Coulomb 
approximation since the ionization potentials are not always known and 
a less empirical approximation scheme is required A theory using 
screening constants, which is described in the following appendices, 
is capable of dealing with these cases. 


26 





















DROGEN IN 


8 9 10 11 12 13 


» 

») 

►> 

»> 

0 


> 

) 

) 

) 

) 


1.232 (3) 
4.674 (4) 
2.299 (4) 
1.286 (4) 
7.799 (3) 
5.007 (3) 
3.357 (3) 
2.329 (3) 
1.663 (3) 
1.215 (3) 
9.065 (2) 
6.882 (2) 


7.147 (4) 
2.811 (4) 
1.426 (4) 
8.187 (3) 
5.007 (3) 
3.323 (3) 
2.267 (3) 
1.597 (3) 
1.156 (3) 
8.550 (2) 
6.449 (2) 


4.374 (4) 
1.773 (4) 
9.226 (3) 
5.414 (3) 
3.422 (3) 
2.278 (3) 
1.577 (3) 
1.126 (3) 
8.248 (2) 
6.170 (2) 


2.797 (4) 
1.163 (4) 
6.182 (3) 
3.876 (3) 
2.376 (3) 
1.605 (3) 
1.126 (3) 
8.136 (2) 
6.023 (2) 


1.856 (4) 
7.880 (3) 
4.268 (3) 
2.595 (3) 
1.692 (3) 
1.158 (3) 
8.219 (2) 
6.002 ( 2 ) 


1.270 (4) 
5.493 (3) 
3.024 (3) 
1.865 (3) 
1.232 (3) 
8.527 (2) 
6.116 (2) 


14 15 16 17 18 19 


8.928 ( 3 ) - 

3.923 ( 3 ) 6.425 ( 3 ) ..... 

2.191 ( 3 ) 2.863 ( 3 ) 4.717 ( 3 ) - 

1.368 ( 3 ) 1.619 ( 3 ) 2.128 ( 3 ) 3.528 ( 3 ) 

9.139 ( 2 ) 1.022 ( 3 ) 1.217 ( 3 ) 1.609 ( 3 ) 2.679 ( 3 ) 

6.394 ( 2 ) 6.899 ( 2 ) 7.763 ( 2 ) 9.291 ( 2 ) 1.234 ( 3 ) 2.066 ( 3 ) 
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TABLE A-3 


BOUND-BOUND TRANSITION PROBABILITIES A(p,q) FOR THE IMPORTANT 
LINES OF LiI, Lill, Cl, CII, CIII, HI, Nil, NIII, 01, 011 AND OIII 

___ __ Li I Lines _ 


Transition 

J 2 

J 1 


2 

O’ 

A(10 7 sec S 

2p 2 P°-2s 2 S 

1/2 

1/2 

6709.8 

5.424 

3.63 


3/2 

1/2 

6709.6 

5.424 

3.63 

3p 2 P°-2s 2 S 

1/2 

1/2 

3233.6 

0.01641 

0.098* 


3/2 

1/2 

3233.6 

0.01641 

0.098* 

4p 2 r°-2s 2 S 

1/2 

1/2 

2742.0 

0.01392 

0.136* 


3/2 

1/2 

2742.0 

0.01392 

0.136* 

5p 2 P°-2s 2 S 

1/2 

1/2 

2563.0 

3.90xl0' 3 

0.0468* 


3/2 

3/2 

2563.0 

3.90xl0' 3 

0.0468* 


Lill Lines 


Transition 

x(X) 


2 

a 

A(10 7 sec’ 1 ) 

ls2p 1 P°-ls 2 *S 

'.99.3 

6 

3.72.;10* 3 

. ..... . . 

1.90* 

ls3p 1 P°-ls 2 l S 

178.0 

6 

9.647xl0' 3 

691 

ls4p 1 P°-ls 2 l S 

1_ 

171.6 

6 

3.219xl0* 3 

257 


Cl Lines 


2p( 2 P°)3s 3 P°-2p 2 3 P 

1657 

18 

.. “ 

0.01397 

1.24* 

2s2p 3 3 D°-2s 2 2p Z J P 

1561 

30 

00.9563 

101 

2s2p 3 3 P°-2s 2 2p 2 3 P 

1329 

18 

0.01284 

2 .21* 

2p( 2 P°)4s 3 P°2p 2 3 P 

1280 

18 

0.001488 

0.287* 


*below 907. maximum value 
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CIII Lines 


2s( 2 S)2p 1 P°-2s 2 l S 

977 

6 

0.5723 

248 

2s( 2 S)3p l P°-28 2 *S 

386 

6 

0.05579 

392 

2s( 2 S)4p 1 P°-2a 2 l S 

310 

6 

0.006529 

88.5 


NI Lines 


2s 2 2p 2 ( 3 P)3s A P- 





2.V V 

1200 

12 

0.007382 

0.860* 

2p 2 ( 3 P)4s 4 P-2p 3 V 

964 

12 

0.0007751 

0.175* 

2p 2 ( 3 P)3d 4 P-2p 3 4S° 

953 

120 

0.04008 

93.7 

2s2p 4 4 P-2s 2 2p 3 4 S° 

1135 

1135 

0.5631 


2p 2 ( 3 P)5s 4 P-2p 3 4 S° 

911 

12 

0.0003190 

0.0853* 


*Below 907. maximum value 
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TABLE A-3 (continued) 
Nil Lines 
















TABLE A-3 (continued) 


Oil Lines (C.A.) 


Transition 

x(8) 


2 

or 

A(10^sec *) 

2s2p 4 4 P-2s 2 2p 3 4 S° 

834 


0.2670 


2s 2 2p 2 ( 3 P)3s 4 P- 





2s 2 2p 3 4 S° 

539 

12 

0.006351 

8.2* 

2p 2 ( 3 P)3d 4 P-2p 3 4 S° 

430 

120 

0.005050 

128 | 

2p 2 ( 3 P)3d 4 D-2p 3 4 S° 

430 


0,005356 



0111 Line s 


2s2p 3 3 D°-2s 2 2p 2 

3 P 

834 

30 

0.1417 

98.4* 

2s2p 3 3 P°-2s 2 2p 2 

3 P 

703 

18 

0 1538 

179 

2s2p 3 3 S°-2s 2 2p 2 

3 P 

508 

12 

0,1643 

1014 

2s 2 2p( 2 P°)3s 3 P°- 
2s 2 2p 2 3 P 


374 

18 

0 006293 

48.6* 

2s 2 2p( 2 P°)3d 3 D°- 
2s 2 2 P 2 3 P 


306 

135 

0.005453 j 

346 


*Belov 90% maximum value 
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TABLE A-4 

TRANSITION PROBABILITIES FOR RESONANCE LINES 


Element 

Transition 
(Upper - Lower) 

J 2 * J 1 

X(X) 

a 2 

A(10 6 sec' 

Rb 

I 

4p 6 ( 1 S)5p 2 P°- 
4p 6 ( 1 S)5s 2 S 

1/2 - 1/2 

7947.60 

8.82 

35.5 

Sr 

I 

5s( 2 S)5p l P°- 
5s 2 L S 

1 - 0 

4607.331 

4.72 

195 

Sr 

II 

4p 6 ( 1 S)5p 2 P°- I 
4p 6 ( 1 S)5s 2 S 

1/2 - 1/2 

4077.714 

4.92 

147 

Zr 

I 

4d 2 5s(a 4 F)5p 3 F°- 

4 - 4 

6121.49 

3.19 

52.5 



4d 2 5s 2 3 F . 

2 - 2 

6134.58 


49,6 

Zr 

II 

4d 2 (a 3 F)5p V- 5 
4d 2 (a 3 F)5s 4 F 

1/2 - 4 1/2 

3391.96 

3.59 

186 

Nb 

I 

4d 3 5s(a 5 F)5p 6 F°- 
4d 4 (a 5 D)5s 6 D 

5 1/2-4 1/2 

4058.933 

3.39 

410 

Mo 

I 

4d 5 (a 6 S)5p 7 P°- 

4 - 3 

3798 359 

3.35 

123 



4d 5 (a 6 S)5s 7 S 

3 - 3 

3864.115 


117 




2 - 3 

3902.968 


114 

Tc 

I 

5d 5 5s(a 7 S)5p 6 P°- 

3 1/2-2 1/2 

4298.2 

3.15 

160 



5d 5 5. 2 6 S 

2 1/2-2 1/2 

4263.4 


164 




1 1/2-2 1/2 

4239.4 


167 

Ru 

I 

4d 7 (a 4 F)5p 5 D°- 
4d 7 (a 4 F)5s 5 F 

Entire 

Multiplet 

3799.124 

3.27 

120 
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TABLE A-4 (continued) 


Transition J o ' J i n .1 

Element (Upper - Lower) A(A) a AflO sec ) 


Rh I 

4d 8 ( 3 F)5p V- 

3 1/2-4 1/2 

3692.357 

3.21 

115 


4d 8 ( 3 F)5s 4 F 

2 1/2-3 1/2 

3657 987 


1082 

Ag I 

o 

w 

•v 

1 1/2- 1/2 

3280 682 

2.95 

169 


4d 10 ( 1 S)5s 2 S 

1/2 - 1/2 

3382.890 


154 

Cd I 

4d l0 5«O 2 S)5p *P° 
4d 10 5 , 2 ‘ S 

- 1 - 0 

2288 7 

1 99 

671 

Tn I 

5s 2 ( 1 S)6s 2 S- 

1/2- 1 1/2 

4511 310 

0 890 

39 2 


5s 2 ( 1 S)5 P 2 P° 

1/2 - 1/2 

4101 764 


26 l 

Sn 1 

5s 2 5p( 2 P°I)6s 3 P C 

’-1-2 

3175 046 

0 202 

10 6 


2 2 3 ^ 

5s Z 5p 2 3 P 

0 - 1 

3034 120 


29.2 

Cs I 

5p 6 ( 1 S)6p 2 P°- 

1 1/2- 1/2 

8521 10 

9 94 

32.5 


5p 6 ( l S)6s 2 S 

1/2 - 1/2 

8943 50 


28 1 

Ba 1 

6s( 2 S)6p 1 P°- 
6s 2 l S 

1 - 0 

5535 484 

5 78 

138 

Ba 11 

5p 6 ( 1 S)6p 2 P°- 

1 1/2- 1/2 

4554 033 

5 96 

127 


5p 6 ( 1 S)6s 2 S 

1/2 - 1/2 

4934.086 


100 

La I 

Sdesfa^Hp 2 D°* 

2 1/2-2 1/2 

7270.11 

5.10 

50.0 


5d6s 2 2 D 
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TABLE A-4 (continued) 


Element 

Trend t Ion 
(Upper - Lower) 

J - J 

2 J 1 

x<X) 

o 2 

A(10 6 eec’ 1 ) 

Rh 1 

4d 8 ( 3 P)5p V- 

3 1/2-4 1/2 

3692,337 

3.21 

115 


4d 8 ( 3 F)5e 4 F 

2 1/2-3 1/2 

3657,987 


1082 

Ag I 

4d 10 ( l S)5p 2 P°- 

1 1/2- 1/2 

3280 682 

2.95 

169 


4d l0 ( l S)5e 2 S 

1/2 - 1/2 

3382,890 


154 

Cd I 

4d 10 5eO 2 S)5p l P° 

- 1 - 0 

2288.7 

1.99 

671 


4d l °5. 2 l S 





In I 

5e 2 ( 1 S)6e 2 S- 

1 /2- 1 1/2 

4511.310 

0 890 

39 2 


5e 2 ( 1 S)5p 2 P° 

1/2 - 1/2 

4101 764 


26,1 

Sn I 

5e 2 5p( 2 P°l)6e 3 P°- 1-2 

3175,046 

0 202 

10.6 


5e 2 5p 2 3 P 

0-1 

3034.120 


29.2 

Ce I 

5p 6 ( l S)6p 2 P°- 

1 1/2- 1/2 

8521.10 

9.94 

32.5 


5p 6 ( l S)6e 2 S 

1/2 - 1/2 

8943.50 


28.1 

Be I 

6 e( 2 S)6p l P°- 

1 - 0 

5535.484 

5.78 

138 


6 e 2 l S 





Be II 

5p 6 ( l S)6p V- 

1 1/2- 1/2 

4554,033 

5.96 

127 


5p 6 ( l S)6e 2 S 

1/2 - 1/2 

4934.086 


100 

U I 

5d6e(e l D)6p V- 

2 1/2-2 1/2 

7270.11 

5.10 

50.0 






A. 2 Radiative Bound-Free Transitions 


In order to calculate the coefficient a(p) for radiative capture of 
electrons into state p 

X + + e - X(p) + hv 

we must make some assumption about the electron velocity distribution. 
Since the electrons exchange energy readily in elastic collisions, they 
may be described by a Maxwellian distribution characterised by an elec¬ 
tron temperature T g . The recombination coefficient is accordingly 

“(HV ■ ^ (; yVy 3 %, exp(i p /kT e ) 

x J Q V (P) (hv) 2 exp(-hv/kT e ) [l +X(v)] d(hv) 

P (A-12) 

where 1 ^ is the ionisation potential of the p C ^ level, hv is energy of 
the radiation photon and Q v (p) is the photoionization cross section 
corresponding to the process 

X(p) + hv - X + + e (A-13) 

The rate of population by radiative capture is 

a(p|T.) n # n(X + ) (A-14) 

where n # and n(X + ) are respectively the electron and X + ion number den¬ 
sities. 
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The rate of depopulation by photoionization may be written 


n(p) 0(P) 

where ^ 

P(P) - J Q V (P) hv (hv) d(hv) (A-15) 

X P 

The quantal description of bound-free transitions is entirely similar 

to that of bound-bound transitions except that the radial wave function 

R(n'f'|r) in (A-ll) becomes R(e'C'|r), describing an electron of positive 

energy e' moving in the field of the ion. Thus, the cross section for 

the transition of an electron from the bound state (ni!) to the continuum 

state of energy €' is - 

4naa 

Q v (ni!;€’) - 

X'-2±l 

where is the ionization potential, C^, are simple algebraic quantities 
which depend upon the initial and final atomic configurations and 
g(nX;e'X') is a radial transition integral 

g(ni;€’i') - \ij r 3 R(n.eJr) R(e'i|r)dr (A-17) 

o 

Although exact analytical expressions are available for the matrix 

elements describing the photoionization of the excited states of hydrogen 
(A-9) 

(Gordon / ), they are complicated and have not yet been given a numerical 
form except for the lowest lying levels. Giovanelli^ A " 10 ^ employs an 
approximate formula due to Gaunt^ A " 11 ^ for the photoionization cross 


< a -i. 



t 
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which in numerical form is 


Q v <p> 


2.82xl0 29 
5 3 


(cm 4 ) 


(A-19) 


P v 

This choice has also been made by Elwert^ A ~* 2 ^ and Ivanov-Kholodnyi 
et al .. but it is of questionable accuracy. Recently, Burgess^" 1 *^ 

has presented tables useful for low velocities of the ejected electron, 
from which more accurate photoionization cross sections may be obtained 
for a hydrogen plasma. His calculations distinguish between the various 
sub-levels and apply to all levels with p £ 12 . 


Using his calculations, the radiative recombination coefficients may 
be calculated according to (A-12) and Burgess has given values of a(nX|T £ ) 
for electron temperatures of 10 , 000 °K and 20 , 000 °K. He do not wish to 
distinguish between the individual sub-levels (except possibly in the 
latest stages of the plasma decay). For a hydrogenic plasma of nuclear 
charge Z, after suitably averaging over the sub-levels, we obtain the 
formula(Seaton ( A * 15 )) 5 

a(p,Z|T # ) = DZ x p 3/2 S p (X) (A-20) 

where 

D - ^ 5 ^ Oj c = 5.197xl0" 14 c » 3 IK 1 (A-21) 

X.“§(A-22) 

and 

S p (x) = S<°> (Xp) + x' l/3 s (1) (X p ) + x‘ 2/3 S (2) (Xp) (A-23) 
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with 

S (o) (x) exp (x^ Ei(x) 


(A-2A ) 


S< 1 > (x) 0.1728 x l/3 (u + l)‘ 5/3 (u - l)e‘ XU du (A-25) 


S (2) (X) - 0.0496 x 2/3 f (u . l)‘ 7/3 (u 2 t | u + l)e" XU du (A-26) 

J o 


u M 

p" 


r 


■ dv (A-27) 


Seaton^has published tables of xS^° 7 (x), xS^.fx) and xS^(x), iroro 
which the radiative recombination coefficients a(p,2|T L ) may be calculated 
readily. Values obtained lor hydrogen (Z 1) in this way have been provided 
by Bales and are shown in Table A-^>. Seaton also presents tables similar 
to those for the xS^ m \x) lrom which the total recombination coefficient 
to all levels 


“> I V 

p'=l 

may be calculated directly. 


(A-28) 


For more complex systems the only method generally applicable is an 
extension of the Coulomb approximation. This yields formulas identical to 
(A-16) and (A-17) except that n is replaced by an effective quantum number 
n* z./\j I n £ > x being the residual charge on the atom when the electron 
has been ejected. 
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The indictee give the power* of 10 by which the entries oust be nultiplied. 














According to Burgess and Seaton 

g(n*';eV) • G(n * / ; ^ cos n [n* + ( i(e') + X(n*2;e'2')] (A-29) 

r(n*,2) 


G(n*2;t'0 




a 

l+t 'n* u xr 


-TP,i,(n*) (A-31) 


lffc'n* 

n(e') is the continuum phase shift obtained by extrapolation of the quantum 
defects t i(I .) of the energy levels of the bound states. 


£(n*» i) 


2 \iU) 
*3 Se 


(A-32) 


b /2* e :t 


■ (n* > l + 2) 


C 22'’ '* ^22' have been tabulated b X Burgess and Seaton for 1 £ n* £ 12 

and f, }' £ 3. For values of n* < 2 + 2, additional tables are presented 
which permit the calculation of g(n*i;€'-f') for transitions from specific 
bound states to the continuum. 


As an example of the application of Burgess and Seaton's^ A ~ 
formulas and tables to the treatment of complex atomic systems, they have 
been used to calculate the photoionization cross sections and radiative 
recombination coefficients for some of the observed states of aluminum. 
The ground state configuration of Al I is 3s 2 3p and the observed excited 
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2 

states, with the exception of one 3s3p configuration, are configurations 
of the form 3s^nfwhere 3 < n < 11 and 0 < l < 3 (cf. C. MooreThe 
latter configurations consist of a single (nf) electron outside a closed 
shell for which the Cassume the particularly simple form (selection 
rules on the azimuthal quantum number allow only transitions where 
l' -- l 1 l)i 


Ul l 

C jU1 w 2f+1 * C f-1 * 2f+1 


(A-34) 


Table A-6 gives the Burgess and Seaton coefficients for the allowed 
transitions in A1I. These coefficients were used to calculate the func¬ 
tional parameters G^^,(n*), , (n*) for individual transitions 

The results, together with the ionization potentials and quantum defects 
of the bound states are listed in Table A-7. The £(n v -,£), also 

given in Table A-7, were determined from the curves of Figure A-l which 
shows plots of the continuous phase shifts for each(nl) series. 


For small positive electron energies e'» the phase shift may be taken 
as a linear function of e'. The extrapolated quantum defects determined 
from the curves of Figure A-l, lead to the following approximate expres¬ 
sions for u(e'): 


H(€') - 1.755 


(ns series) 


H(€') = 1.264 

- 0.065c' 

(np series) 

(A-35) 

H(c') = 0.985 

+ 3.3c' 

(nd series) 


H(c') - 0.044 

+ 0.26c' 

(nf series) 
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Photoionization cross sections were calculated for all but the 3s3p 2 
and nf levels, for ejected electron energies e' lying in the interval 
(Rydbergs) 0 5 e'— 0.5 for the nd levels, 0 1 «'l 0.3 for the ns levels 
and 0 1 £'l 0.2 for the np series. For electron energies in excess of 
the maximum values indicated, it was found that the approximate formulas 
led to incorrect results. This may in part be due to the assumption of 
a linear dependence for the n(c'); however, it should be noted that the 
approximations inherent in the derivation of (A-29) limit the validity 
of the calculations to small ejected electron energies. 

The photoionization cross sections obtained are listed in Table A-8. 

2 2 

The 3s 3d and 3s 4d cross sections were calculated under the assumption 
that n* ^ f+2 since explicit expressions are not tabulated by Burgess and 
Seatonfor nd - e'p and nd “* e'f transitions when n* < 1+2. The 
resulting cross sections are therefore of lesser accuracy than the remain¬ 
ing cross sections shown. 

The radiative recombination coefficient may be calculated from Bq. (A-12) 
once the photoionization cross sections are known as functions of the ejected 
electron energy. Thus, in terms of electron energy for a thin plasma 

v■ 7(- **> 

Instead of evaluating the integral of (A-36), the recombination coefficient 
can be determined approximately by use of the expression 

a(n/|T # ) * v Q c (nf) (A-37) 


4f 




be Multiplied . €* is la Kydbsrgs. 



















































where Q e (n£) the cross section for radiative recombination Is given by 
detailed balancing 

(I ,+0 2 

Q e (ni) - (2i+l) —“- Q v (n t) (A-38) 

me £ 

and v is the average electron speed 

/.« Y ' 2 

(A-39) 

Substituting (A-39) and (A-38) in (A-37), we find 

«MV- (efr ) 1/2 ^ V *> 2 Q v (n/) 

6 . Q <»/> 

- 2.046x10° (2/+1) (I p +€) ~i / 2 


The radiative recombination coefficients calculated from (A-40) are 
given in Table A-9 for electron temperatures of 6190°K and 24,760°K. 



RADIATIVE RECOMBINATION COEFFICIENTS 
FOR A1 I a(t»i|T)ln cm 3 /»ec) 


T @ - 6190°K T g - 24,760°K 


5 

s 

P i 

d 1 

1 _ i_ 1 

s 

P 

d 

3 


2.11(-13) 

6.48(-14) 



8.04(-15) 

1.94(-14) 

4 

9.80(-15) 

3.90(-14) 

2.36(-14) 


1.99(-16) 

1.73(-15) 

1.53(-14) 

5 

1.86(-16) 

1.32(-14) 

1.04(-16) 


5.34(-17) 

7.68(-16) 

2.56(-15) 

6 

7.28(-16) 

5.62(-15) 

9.39(-16) 


9.24(-17) 

2.68(-16) 

6.80(-16) 

7 

4.94(-16) 

5.06(-15) 

1.33(-15) 


1.0K-16) 

2.96(-16) 

2.21(-16) 

8 

1.74(-17) 


1.28(-14) 


8.35(-17) 


1.31(-16) 

9 

8.35(-18) 


1.10(-16) 


7.63(-17) 


1.03(-16) 

10 

2.22(-18) 


8.66(-16) 


B.04(-17) 


8.46(-16) 

11 



7.33(-16) 




7.52(-16) 


The brackets give the power of 10 by which the 
entries must ',<3 multiplied. 






A.3 Electron Collision Induced Transitions Between Bound States 


Let q £ (p,q) be the cross section for the excitation to the q state 
by the impact of an electron of energy c with an atom in the p state. 

The rate of depopulation of the p th state by electron impact is 


" e n<P> ^ C(P,r) 

r > p 


(A-41) 


where the collision probability C(p,r) depends upon the electron temp¬ 
erature and is obtained by Integrating the product of the electron 
velocity and the cross section over the Maxwellian distribution accord¬ 
ing to 

. h oo 

C(p,r) - (KTg)”^ 2 q € (P.r) € exp(-e/kT # ) dc (A-42) 


The probability of de-excitation to a lower level is similarly 


q 


(A-43) 


3 




The direct and inveree collision probabilities are related through the 
principle of detailed balancing: 

w f I ■ I 

C(p,q) - ^ exp ^ C(q,p) (A-44) 


The p state may also be populated by electron impact induced 
excitation and de-excitation and the appropriate rates are 


-.1 
q < p 

.. I 


n(q) C(q,p) 


n(r) C(r,p) 


(A-45) 


(A-46) 


r > p 

There are various simple ways of computing the probabilities of 
excitation and de-excitation transitions induced by electron-atom 
collisions (c.f. Allen**’ 4 * none of which is entirely satsifactory. 
Those quantum mechanical treatments in which some reliability may be 
placed (e.g. Mott and Massey** *®*, Bates et al** 19 *, Milford and co 
workers**' 20 *, McCrea and Me Eirgan**" 21 *, and Me Carrol**’ 22 * have 
been based on the Born approximation, and therefore, are valid only 
for high electron impact energies. 
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A simplified formula describing the cross-section for excitation from s 
q level to a p level induced by the form (c.f. Ailen^ A "^) 

Q_(q»P) - 1.44 x 10 ll *a 2 (Ch) 2 f(q,p) ~( — - j) (cm 2 ) (A-47) 

€ ° € € « « € 

q.p 

where e denotes the impacting electron energy, and ^ is the ex¬ 
citation energy, both expressed in Rydbergs (e^ ^ - c^- . f(q,p) represents 

the oscillator strength for absorption, and is related to the spontaneous 
transition probability A(P»q) by 

3 b> 

f(q.p) SV? z^ A <P*q) (A-48) 

8* « V q 

Equations (A-47) and (A-48) may be used to calculate colllslonal transition 
probabilities for the atomic systems for which the Spontaneous transition 
probabilities are given in Tables A-l, A-2, and A-3. If Q ( (p,q) denotes 
the corresponding colllslonal de-excitation cross section, then (A-47) * 
and (A-48), with the aid of the following expression (c.f. Alien^ 4 " 4 ^): 

Vp Q c (p,q) “ Vq Q e (q,p) <A " 49) 

«*» p *<P.q) * U qf(q*P) (A-50)' 


yields 


Q (p.q) - 5.78 x 10“ xa 


£(q 


"(r) 


(A-51) 


Expressions similar to (A-47) and (A-51) with values of A(p,q ) obtained 
by the formulas of Mensel and Pekerie^ A " 7 ^ for hydrogen were employed 
by Ciovaneili^ A " 10 ^ and Ivanov-Kholodnyi et al^ A ~*^. 
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Perhaps the most accurate method of computing colllalonai transition 
probabilities is the classical method of Gryzinski^’^, which represents 
an Improvement over the original formulation of Thomson^, in that 
it takes account of the fact that the atomic electrons are moving. The 
cross-sections so desired fall off too rapidly at high electron energies, 
*» t * instead of the correct (e 1 In e), but their accuracy probably 
suffices, the effect of excitation and de-excitation not being a pri- 
mary one. 


Considering in some detail the classical transfer of momentum during 
a coulomb collision between an incident electron of energy e and an atomic elec¬ 
tron having binding energy Gryzinakl derived the following approximate ex¬ 
pression for the differential cross-section cj(/*) corresponding to a 
loss in energy of the incident electron of magnitude : 


o(Afc) . 




_\ 3/2 




1 - 


-f -£)[(»■ + r)f- 


(A-52) 

The total cross section for a collision involving an energy loss /2c subject 
to the restriction c < c < c is then given by 

p » - q 

P 

Q-<€ .« ) • fo(<*)d(AE) - Q<€ ) - Q(« ) (A-53) 

€ p q J p q 

<1 

where Q(c^) the cross section for a collision with electron loss greater 

than c t s defined by Equation (A-4V). Thus when Ae - c , as in the case 
P P 

for hydrogenic orbits. 


1/2 
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<»’> 

J a* * « 

p 


/ *B^ /2 2 

H 1 " TV (** ) 12* * « 


(A-54) 


The collision transition probability is given by 

\l/2 

rt>T I ft <, . _ . 

“ q’ P 


«>•<» - (& 12 <“,)‘ 3/2 J P W %> « •' ,/k \c 


(A-55) 


Inserting Q f («^; *p) lesda to the following collision prob¬ 
abilities for hydrogen, where * ■ 1/p 2 and ■ 1/q 2 

o(p.o . V ft *■* x 105 y} «y 

e q 1 qp T # 


(A-56) 


(A-57) 


(x 


*• 4 


5 [2 -(g + y(l + 4g)] (y - l) 1/2 (g + 1) 


* 1/2 (g ♦ y) 3/2 


^C-3g ♦ y(3 + 4g)] y 1/2 


(g + y) 3/2 


X < 1 ♦ g 
(A-5(0 


y > l + 8 

CA-59) 


C(p,q) has been calculated for hydrogen by Bates using equation 
(A-56). The results are given in Table (A>10) for p & q less than 
or equal to 10. 
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its give the powers of 10 by which eh* entries wist be suitiplied. 





























TABU A-10 (continued) 



The bracket* give the power* of 10 by which the entrie* mist be miltiplied. 

























A-10 (continued) 



Ihc bracket* give the power* of 10 by vbich the entries mist be multiplied. 
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For non-hydrogenic systerna having known energy level* € p , the more 
general expressions of (A-S2) and (A-53) must be used for the calcula* 
tion of collision probabilities. 


60 


.i .4 Electron Collision Induced Bound-Free Transitions 


Let q(p,e) dc be the cross section for the formation of electrons 
with energy In the range e, c + de by collision of an electron having 
energy E with an atom in the p state. The total Ionisation cross section 

is 

E " I p 

Q e (p) • \ J q E (P.e> (A-60) 

and the probability of ionization by collision is 


D(P> "(to) (kT e >’ 3/2 J q E (p)E «xp<-E/kT e ) (A-61) 

The rate of depopulation of the (p) state due to electron impact ioniza¬ 
tion is 

n(p) n g D(p) (A-62) 


The inverse process is three-body recombination: 


X + + e + e - X(p) + e. (A-63) 

The efficiency of the three-body process may be written in the form 
00 • 

F < p > “(^) (W e )“ 3 j J * 2 exp(-« 1 /kT e ) exp(-e 2 /kT e ) 


x Y(p|« 1 ,* 2 )d« 1 de 2 


(A-64) 


where Y(pI*i* 2 ) is th « probability that an atom is found in the (p) state 
as a result of the collision with a positive ion and two alectrons with 
energies « 1 and c 2 respectively. The probability is relatad to tha in¬ 
verse process of electron impact ionisation by detailed balancing. Thus 


61 



16 rrnu). 

<I p + «! + « 2 > <l E (pl e i> " 1- *1*2 Y(pI« 1 . « 2 ) < A ’ 65 ) 

P 

Where ou^ is the statistical weight of the continuum electron and 

E " X p + *1 + *2 ' (A-66) 

The rate of population of level (p) by the three-body recombination 
process (Equation A-63) is then 

n* n(X + ) P(p) (A-67) 

Alternatively, the probability of three-body recombination F(p) may 
be obtained by a simpler method noted by Ivanov-Kholodnyi, et al ^ A 
They remark that in thermodynamic equilibrium, the number of particles 
leaving a given quantum level via any specific process is equal to the 
number of particles entering this state via the inverse process. It 
follows that 

F(p) - <P(P) D(p) (A-68) 
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where 


, /2nmk T \-3/2 

cp(p) , p 2 Q— t) exp(I p /kT e ) (A-69) 


As in the case of excitation, the cross sections decrease too rapidly 
as the energy of impact increases. Since ionization and its inverse 
are the controlling processes, the error may be serious and we have 
therefore given some attention to the problem of how we may obtain im¬ 
proved accuracy, whilst still retaining the universal applicability and 
the simplicity of the Gryzinski method. A new universal method that 
has resulted is described below. 


Again for reasons similar to those mentioned above for bound- 
bound collision probabilities, the most accurate method for computing 
electron Impact ionization probabilities, which is of universal ap¬ 
plicability is that of Gryzlnaki^“ 23 \ Giovanelli^” 1 ^ and 
Blwert^*" 12 ^ used the classical Thomson formula which can be expected 
to give somewhat less accurate results than that of Grysinski. 
Ivanov-Kholodnyi et ai^ A ^ employed approximate expressions simi¬ 
lar to those for the colllslonal excitation cross-sections. 


i 3 


The classical Thomson cross-section (cf. Thomsonor 
Fowlerf 0 r the ionization of an atom having n Q valence elec¬ 
trons by an electron with energy €, with the corresponding ejection 
of an electron from a p-state into an energy Interval £, £ + d£ la 


4 "o* V 

% (P,C> d£ = - dC 

‘ *« ♦ «/ 


(A-70) 


The Thomson total ionization cross-section is then 


4 nx a p€-e f 

% <p) * / p - 

* * J « «+ « r 


illii 


(A-71) 


a Q is the radius of the first Bohn orbit in hydrogen. 
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As pointed out in the previous section, the Grysinski formulas 
represent an improvement over the Thomson derivation in that account 
is taken of the motion of the bound atomic electrons. According to 
Grysinski the total ionisation cross-section for an electron-electron 
collision in which the impacting electron suffers an energy loss AE 
greater than the threshold ionisation energy s is given by 

Q(« ) . J »«*) i (AE) (A- 72) 

* *P 

where AE ^ ^ is the maximum permissible transfer of energy consistent 
with conservation of momentum and energy, and a (AE) is the differ¬ 
ential cross-section described by Iquation (A-S2). (The expression 
for Q(c) given in (A-54) was obtained from the integration indicated 
by (A-72)). Thus the cross-section for ionisation from the p level 
is just (A-54). The ionisation transition probability for a Maxwel¬ 
lian electron energy dietrlbution is given by 

^ CA-7* 


Making the substitutions x - e/« p we find from (A-54) and (A-73) that 



(cm 3 /sec) x > 2 




x < 2 


In the case of hydrogen, « p ■ 1/p 2 so that 



Values of D(pf T^) in hydrogen (provided by Bates) are given in Table 
A-ll for p leas than or equal to 10. For the more general case, Equa¬ 
tion (A-74) must be used to compute D(nl | T e >. 

As in the case of excitation, the cross-sections decrease too 
rapidly as the energy of impact increases. Since ionisation and its 
inverse are the controlling processes at the higher energy levels, 
the error may be serious and we have therefore given some attention 
to the problem of how to obtain improved accuracy, whilst still re¬ 
taining the universal applicability and the simplicity of the Grysin- 
ski method. 


Rewriting Equation (A-71) for the Thomson ionisation cross-section 
in the fora 


(«) 


, 2 1 




and introducing the dimensionless parameter X - s/I 


CA-76) 
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(A-47) becomes 


1 <M XI > . iu.n 

-_ -■ An a_ 7 ( 1-7 ) 


(A-77) 


The right-hand side of (A-77) is a universal function of X, independ¬ 
ent of the particular ionising process. 


It is known from the asymptotic form of the Born approximation 
(cf. Kott and Massey^ A " 13 ^) that for large values of X, behaves 
asymptotically according to 

Q t - X ' 1 In X (A-78) 

so that (A-77) is incorrect in the limit of high impact energies. 

The modification of (A-71) by Grysinski^ A " 23 ^ also has this defect. 


We shall attempt to improve (A-77) by introducing a correcting 
factor f (X) such that 




f(X) 


(A-79) 


4* a g(X) 


and determine f (X) by appeal to experimental data. 

In figure (A-2), all the available experimental data^ A-2 ^ thru 
' on the electron impact aingle ionisation of atoms has been 
used to determine l2 V°o “ a function of X, adopting for n Q the 
number of electrons in the valence shell. If f (X) exists such that 
(A-79)is true, this procedure would produce a single curve. Instead, 
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there results distinct curves which are similar in shape but have an 
appreciable spread in magnitude of the cross-scctions. This spread 
represents the possible error that would be a consequence of the use 
of the mean of the curves of Figure (A-2) as a universal function. 

He may reduce this factor by taking partial account of the elec¬ 
tron structure of the atomic systems, it being an oversimplification 
to adopt for n simply the number of valence electrons n o> The signifi¬ 
cant quantity is actually the equivalent number of dispersion elec¬ 
trons for transitions into the continum. 

two alternative ways. He may 
the same magnitude at large 
that the maximum value of 

(A-80) 

(A-81) 

The latter procedure is the simpler and we adopt it. It leads to the 
curves shown in Figure (A-3), all of the normalised experimental data 
for ionisation of neutral atoms lies within the shaded area. No 
curve now departs from the mean of all the curves by a factor of more 
than about 10% (cf shaded area of Figure A-3) and we shall select the 
mean as a universal curve. The normalised curve for ionisation of 


To determine n, we may proceed in 
choose n to ensure all the curves have 
energies or we may observe from (A-71) 
is given by 


and determine n from 
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Figure A-3. Universal Curve for Electron Impact Ionization 
Cross Sections 



of He* ion falls just outside of this area and is also shown in Figure 
A-3. 

The corresponding values of n are given in Table (A-12a). We see 
that there is a pattern to n which allows us to predict it with some 
reliability at least for the lighter elements. Thus, since n appears to 
be an increasing function of atomic number, we can assign integer values 
of n to appropriate ranges of the atomic number as shown in Table (A-12b) 

From the universal curve of Figure (A-3), we can proceed to the de¬ 
termination of a universal ionization probability curve as a function of 
electron temperature T^. Thus from Equation (A-58) we find 

D( V = (m) l/2 ( kT e )" 3/2 /°° Q t <M> Xe” XI/kT e I 2 dA (A-82) 

Using Equation (A-77) 

D(T e ) - (^) 1/2 ( k V“ 3/2 na o 2 f^ 10 k 8(A)e” X/kT ® dX (A-83) 

where T^ = T g /I is the electron temperature measured in units of the 
ionization potential. Finally 

D(T e )T e 3/2 = ^ 1/2 4ji na ^2 J °°Xg(X)e" X/kT e dX (A-84) 

The right hand side of which is a function of the electron temperature 
T^ only and may be computed once for all atomic systems for any given 
value of T^, using the universal curve g(X) given in Figure (A-3). 
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TABLE A-12a 

EQUIVALENT NUMBER OF DISPERSION ELECTRONS, 
n DETERMINED BY EQUATION (A-81) 








DT 3/2 

The function - * ■ — is shown as a function of T' in Figure (A-4). 
n e 

The ionization probability for any element for any temperature may be 
read off the graph, the appropriate abscissa being obtained by dividing 
T^ by the ionization potential I. 


I 

I 

I 
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Figure A-4. Universal lonisetioa Probability Curve 
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A.5 Heavy Particle Collisions 


Consider now transitions produced by collisions with heavy particles 
which we label Y^. Let Q* (p,r) be the cross section for the excitation 
to the r** 1 state by the Impact of a heavy particle Y^ with an atom 
X in the p** 1 state, the energy of relative ontion being e. The rate of 
depopulation of the p C ^ state by heavy particle, impact is 


«(P) Z n(Y t ) Y. P i 
1 r>p 


where the collision probability P^(p,r) is obtained by multiplying the 
cross section by the relative velocity of the collision and integrating 
over ion velocity distribution f(v) to give 

p l<P, r ) -J v Q* (P.r) f(v) dv (A-86) 

If f(v) is the Maxwellian velocity distribution corresponding to temp- 


!(p,r) € exp(-€/kT) dc 


where p is the reduced mass. 


The probability of de-excitation to a lower level is similarly 


n(p) J n(Y^) £ P t (p,q) (A-M) 

i q<P 
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The rates of population ere 


n (Y 1 ) n(q) P t (q,p * (A-89) 

l q<p 


Y n <v y n<r) p i (r,p) <A ‘ 90) 

i r>p 


The ionization and recombination involving collision with ions may 
be formulated in a manner similar to that for electron collisions. The 
rate of depopulation by ion Impact induced ionization may be written 


n(p) (n<X + ) Q(p) + Z n(Y+) Q^p)} (A-91) 

OD 

Q(P)- i-x (W)' 3/2 J QgiP) exp (-E/2RT) dE 




(A-92) 


Q E (p) Is the total ionization cross section of the state (p) by ions of 
energy E. 

The rate of population by three-body recombination 

X + + Y* + e -* X(p) + Y* (A-93) 

may be written 

n(X + ) n(Yj) R(p) (A-9A) 

where R(p) is an expression exactly similar to Equation (A-64) 
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Since for allowed transitions the cross sections for the impact of 
an electron and a singly charged positive ion of the same velocity are 
of much the same sice, excitation and ionization by ion impact is usually 
very inefficient compared to electron impact. For example, three-body 
recombination 

X + + X + + e -* X + + X (A-95) 


will be negligible compared to 

X + + e + e -» X ♦ e 


(A-96) 


Exceptions may occur when the energy defect of the process is 
small as in 


H + ♦ H(n J) - H + + H(n,/') 


(A-97) 


and for transitions between pairs of highly excited states. The effects 
of ion impacts is to strengthen the argument that a Boltzmann equili¬ 
brium prevails for the excited states. The position regarding charge 
transfer reactions is analogous, charge transfer being rapid at low 
velocities only when the energy defect is low. 


We merely wish to provde a few maabers in order to show quantitatively 
that, apart from re-dlstrlbutlng processes, (and transitions between 
pairs of highly excited states), electron collisions are more efficient 
than heavy particle collisions. 
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(A-34) 

Seaton' has calculated the rate coefficient for the re-distribution 

process 

H + + H(2s) ■* H + + H(2p) (A-98) 

o -4 3 

at 10,000 K, and he obtains the very large value of about 5 x 10 cm 

sec Re-distribution may also occur by charge transfer 

H + + H(2s) ♦ H(2p) + H + (A-99) 

but Boyd and Dalgarno^ A ^ have shown that (A-99) Is much less rapid 
than (A-98). 

Excitation and ionization proceed much more slowly and may be 
neglected. Thus, the maximum values of the cross sections for the 
direct processes 

H + + H(ls) - H + + H(2s) (A-100) 

H + + H(ls) - H + + H(2p) (A-101) 

are respectively about 5 x 10 ^ cm^ at about 20 kev impact energy, 
(Batea^ and 1 x 10 ^ cm^ at about 30 kev (Bates^ and 
the cross sections are decreasing rapidly with decreasing energy. 

Similarly, it may be seen from the calculations of Bates and 
(A-38 A-39) 

Dalgarno' that excitation accoispanying charge transfer such 

as 

H + + H(ls) •* H(2p) ♦ H + (A-102) 
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Is extremely Inefficient for temperatures below several kev. 

Thus, for a hydrogen plasma or other singly ionized plasma of a 
single component, processes due to ion impact may be ignored. When more 
than one component exists in a plasma, however, it may be necessary to 
consider heavy particle processes, for charge transfer processes are 
unique in coupling together the two components. This fact may be seen 
from Equation (A-106), by noting that if the heavy body coefficients 
P^, Q i< and are set equal to zero, all the terms containing n(Y^), 
the i.umber density of other ions vanish. Account must be taken not 
only of ordinary charge transfer 

X 4 * + Y + - X + + Y 4 * (A-103) 

but also of radiative charge transfer 

X 44 + Y 4 - X + + Y* 4 + hv (A-104) 

which may be more rapid than the previous reaction. 

An investigation of the possible rates of these charge transfer 
reactions, one of the participants of which is an ionic system in its 
ground state, is necessary, and as an example of the latter process 

He 44 + H ■* Hs + ♦ He 4 + hv (A-105) 

may be treated with high accuracy. 

Special features also arise in the case of a single component plasaui 
in which the component may exist in more than one state of ionization. 


A.6 Tht Basic Equations 


The increase in population per second of a particular level (p) Is 
given by the difference between the rates of population and depopulation 
that we have written in the earlier sections. Thus 

[ Z A(P,q> 0 + X(p,q) ) + "" £ w r X(p,r)A(r,p)] 

L q<P P 

+ n e[? C(p,r)+ q* C(P.q) + D(p)J 

+ Zn(Y 1 ) ^F^p.r) + Zp^P.q) + Q t (p)J 

+ P(p)| n(p) 

-Z|A(r,p) [l + X (r,p)J + n # C(r,p) + ^ n(Y t )P t (r,p)| n(r) 

+ : (q o^ A(pq) + n e C<q i p) + ? n < Y 1 > P i* , » P ‘j} " (q) 

♦ CKP) n # X(n + ) + n* X + y f (p) 

+ n e ^ "i <Y i ) n(X+) R i (p) (A-106) 

He arrive at a set of coupled, differentiel equetions of the form 

MEl + a n(p ) . V a n + b (A-107) 
dt pp VK/ Z- PJ j p 
j * ? 

where the a's and b's can be identified by comparison with the 
preceding equation. 
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APPENDIX I 
A. Naqvi 


ATOMIC TRANSITION PROBABILITIES 
B.l INTRODUCTION 

Consider two energy lsvsls, 1 end 2, of any atom, and lat thair 
nnergies, E^ and E 2 , be such that E^ > Ej. Transitions from 2 to I result 
In enlsslon, and the reverse transitions, from 1 to 2, are, among other 
processes, caused by the absorption of radiation, of frequency v 2l given 
by 

E z -H l = lr\ (b-d 

The amount of energy (In erg cm' 3 sec' 1 ) emitted In a spectral line 
is given by 


E zl - Tz\l (b-2) 

where n 2 Is the number of atoms per cm 3 In level 2. A(2,l) and B(2,l) 

are, respectively, the probabilities (in sec* 1 ) of spontaneous transitions 
and the transitions stimulatod by the presence of an external radiation 
field. I^ is the intensity of this external radiation field at the 
frequency v 2l> 

Similarly the amount of energy (in ergs cm' 3 sec' 1 ) absorbed in a 
spectral line is given by 
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Fi2.- a, B0> 2 )l2\la24^ 


(B-3) 


where B(l,2) is the probeblllty of the upward transition, from level 1 to 
level 2, with the absorption of radiation of frequency v 2 ^« 

The three transition probabilities ere related to each other by 
the following relations due to Einstein. 

(B-4) 


kM- «'•« 


(B-5) 


where end g 2 ere the statistical weights of the two levels. It Is 
therefore sufficient to determine one of the three transition probabilities, 
usually A(2,l). 

The theory for calculating atomic transition probabilities for 

electric dipole radiation Is well known (Condon and Shortley^®* 1 ^, 

( 3 , 2 ) 

Slater ). Several reviews have boon written la recent years also 
(Garstang (B ~ 3) , Uaqvi (B " 4) ). The probability A(2,l) of a spontaneous 
radiative treasltlen Is given by 
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(B-6) 


whare g 2 la the statistical weight of level 2, X is the wavelength of the 
spectral line enltted, S(2,l) is known as the line strength and for LS 
coupling can be written as 

sti.o.so.xj-T&ortnu* 

where Tm and fw are factors which arise as a result of integration 
over the angular part of the wave function. ]f(L) is known as the rela¬ 
tive strength of a line within the aultlplet to which it belongs and 
vf (M) is known as the relative aultlplet strength. Both of these quanti¬ 
ties can be calculated in a closed fore for any given line and any given 
aultlplet. j*(L) has been tabulated by Mbits and Ellason (B-S) and by 
Ruesell (B * 6) . J(M) has been tabuUted by Goldberg**" 7 * nd B " 8) . These 
tables have also been reproduced by Allen* B "*\ Recently Rohrlich* 8 "* 0 ^ 
has derived general expressions, baaed on Racah's aethod, fr<m which 
aultlplet strengths for transitions not covered by Goldberg can be obtained. 

The electric-dipole transition integral a* is given by 

where the P's ore r tines the radial wave functions, and i, is the greater 
of the two aslauthal quantun nuabors and J^. The calculation of the 
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absolute transition probabilities tharafora dapanda upon tha calculation 
2 

of tha transition Integral <x , whieh in turn dapands upon an appropriate 
cholca of tha radial wave functions. 

B.2 THE RADIAL HAVE FUNCTIONS 

The radial wave functions, P(ni;r), arc tha solutions of tha radial 
part of tha non-relativist1c Schroedlnger aquation, 

{ jF * + t l E - V ( r ) ]-^ 7 T L) } P (' 1lt ' ,, ' )ii0 <*-»> 

where E is the eigen value of tha Hamiltonlan for a given stationary atata 
and V(r) is the central field in which tha electron moves. Note that 
atonic units era used throughout this report. 

Tha importance of a knowledge of tha radial wave functions cannot 
be over emphasised. Tha radial wave function, in a nutshell, gives us 
an insight into all other physical properties of the atom. A knowledge of 
tha radial wave functions Is required not only In tha calculation of tha 
tranaltion prokabllltlaa but also photoIonisation cross sections, elastic 
scattering croaa-sectIons, collision cross-sections for excitation of 
atoms, x-ray scattering and numerous molecular and solid stats properties, 
among others. 

Tha most commonly used method Is tha self-consistent field method, 
developed by Hartree and extended by Fock to include exchange effects (for 
a recent review, see Hartree(B-ll)^ These are obtained by numerical 
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solution of a systsa of simultaneous, inhomogeneous, non-linear, intagro 
differential equations, the so-called Fock equations. Efforts have also 
bean made to Include the effects due to configuration mixing and correla¬ 
tion of electrons. (See Hartree^ 8 " 11 ^), The procedure is very laborious, 
and has to be performed separately for each configuration of each atom. 
Since numerical methods are used the results are given in a tabular form, 
and cannot be used for further analytical calculations. It is, therefore, 
not possible to develop analytical methods of great generality, but every 
problem must be solved separately. Thua although the theory was developed 
in all essential details over 30 years ago, a list of atoms and their con¬ 
figurations for which the self-consistent field radial wave functions are 
available is pitifully small. In particular it is clear that if one wants 
to calculste atomic parameters, such as transition probabilities, for a 
large number of atoms one must develop methods of considerable generality 
which can, at one time, be applied to a large number of ceses. 

On account of their obvious usefulness, many attempts have been 
made to obtain analytic wave functions (for example Zener^ 8 ” 12 ^; 

Eckert (B ‘ l3) ; Slater (B " l4> ; Morse, Young and Haurwits (B ‘ l5) ; Tubia (B " l6) ; 
Morse end Yilmas**’’ 17 ^; Breene**" 18 and ftoothan, Sachs and 

Weias* 8 " 20 *). Some of these were merely attempts to obtain analytical 
expressions of which the parameters were determined by comparison with 
the tabular wave functions of the self-consistent field method. 

Since the advent of quantum theory, it has been very tempting to 
use the exact results obtained for hydrogen, for the more complex atoms 
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alio, by tha substitution of althar an "affactiva" principal quantum 

number n* for n, or a "screened" nuclaar charga (* Z - $) for tha 

a a 

actual nuclaar charga Z. Here s a denotes the screening of an electron a 
by ail tha other electrons in tha atoai. Tha eereoalag peraamters hove 
bean determined usually by a semi-empirical method (Pauling (B ‘ 2l> i 
Pauling and Sherman <B " 22) ), although a variational treatment haa been 
applied to some simple cases. The intereat in the uee of hydrogenlc wave 
functiona with a screened nuclear charge has been revived through a method 
for calculating the screening parameters developed by Leyter^*" 23 ^ 

1.3 THE HYDROGENIC HAVE FUNCTIONS 


The screened hydrogenlc redial wave functions, briefly referred to 
as tha "hydrogenlc wave functions", are the solutions of Equation 1-9 with 
V(r) replaced Z^/r, the Coulomb potential due to a screened nuolear charge 
(» Z • s # |), where a^j is a screening parameter and Z is the atomic 
number. They are given by 


P(ntZ w j ; r )= 



•«»r (h-10) 


where a* denotes (n-4-1), the number ef nodes ef the function and ^ is 
a confluent hyp crg con a trlc function defined below, l e c e u se ef the way 
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Z ni occurs In Equation MO, we shall find It convsnlont to substltuts In 
ths shove equation, 




(B-U) 


whore ^ equals 2Z/n and « a (« 2s Q /n) Is a reduced screening parameter. 


We now have 




(B-12) 


1/2 

0 = riaSLl 


# 

,F, ^ ^ ^ ]T\ 4hi r 


(B-13) 


(B-14) 


where 



/ > M«V- « 

*• TT 


(B-15) 
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It la to be noted that t* different for different ni and that 
the radial wave functions (Equations B-10 or B-12) for the sane n but dif- 
ferent i are not strictly orthogonal, It Is possible to make corrections 
due to the lack of orthogonality of the wave functions, but they are not 
considered in this report, 

For subsequent use in calculating the Slater F integrals we write 
below an expression for the square of the radial wave function (Equations 
B-12 and B~14)„ 

> r /) ^ ^ 

an* {> |> l, 

where 

f / i M 

5 (B-i7) 

B.4 RESUME OF LAYZER'S SCREENING THEORY 

Layser (B “ 23) has shown that the eigen values 
the correct quadratic and linear Z-dependence if the 
wave functions satisfy the condition that, 

P (nly;r)= P M («ly;r>0("tyO 


of energy will possess 
one-electron radial 


(B-18) 
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where the superscript H denotes s hydrogenIc wave function end 0(n/ y ;r) 
denotes a tern of order one In and 

'ft- - 0 (f) 

>0 

so that, at least, for large "y , 0(n i y ;r) term cen be neglected. Note 
that our notation differs slightly from that used by Layser. 



Neglecting 0(nl y jr) terms, and applying the variational principle 
Layser has shown that the screening parameter which gives the screen** 

ing of the electron a by ell the other electrons in the atom, is given by 


4 

' 1 * 1 


(B-20) 


where denotes the metrlx elements of the Coulomb interectlon energy 
and denotes the number of electron in the eubshell a. We can use 
several different epproxlmatlons for the Coulomb interaction, for exasqple 
we may or mey not include configuration mixing and we mey either conalder 
the everege energy of a given configuretlon or we may add to it the con¬ 
tributions which come from different spectroscopic terms (characterised 
by quantum numbers S and L), belonging to the same configuration. 


In the work reported here, we have not considered configuretlon 
mixing, but we consider the energy of individual terms as wall as the 
average energy of the configuration. When required, the screening con¬ 
stants which are dependent on SL can be labelled expllcitiy as * n £g{ i » 
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whereas those based cm the average energy will be labelled as * 


ni* 

The term energies as functions of the screening parameters ere 

given by 


Wfo) <B - ll) 

where the summation Is over ell the subshells in the atom. This is a very 
simple formula for calculating the term energies from the screening param¬ 
eters. 


B.5 THE COULOMB INTERACTIONS OF ELECTRONS 

The nmst commonly used approximation for complex atoms (see, for 
example, Condon and Shortleyla the aepereble wave function approx¬ 
imation in which each electron la aeeumed to be in a atetlonery state in 
the field of the nucleus and all the other electrons. The wave function 
of an N-electron atoms is expressed in the form of a determinant, made up 
of the wave functions of individual electrons, ea follows. 



^1“’) u,(a') . 

4/a') 



- 

■ 


(B-22) 






I 

I 


where u is the weve function of the aingle electron problem, a 1 denotes 
the eet of quantum numbers (n 1 m* m*) and subscripts of u, from 1 to 
N, denote the electron to which this wave function belongs. The deter- 
mlnantal form Insures the antisymmetry property of the wave functions. 

For a many-electron wave function of the above form, Slater 
has shown that the matrix elements of the mutual Coulomb Interactions of 
the electrons, l/r A j, are sums of the matrix elements of the corresponding 
two electron problem. 


Ignoring configuration mixing, the average electrostatic energy 
of a configuration may be written as 






(B-23) 


where 


= O' T £*?)%( ("If" W 


and the two electron Interaction (a,0) Is given by 

=r/f*(w P ) F h c«,?) 


(B-25) 


where the prime over the summation sign Indicates that k does not assume 
all Integral values, but for F k all even integral values and for G k either 
all even or all odd Integral values. ^ denotes the lesser of and 1^. 
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The coefficients f k (ig /p) «nd g k (ig 1^) arise from the angular 
and the spin part of the wave functions. These can be calculated by the 
well known matrix methods. For details of these calculations the reader 
is referred to Condon and Shortley^ 8 ’ 1 ^ and to the more powerful methods, 
due to Racah^ 8 '^, based on tensor operators. These coefficients are 
listed in Tables B-l end B-2. 

It It 

F and G are integrals over the radial part of the wave functions, 
also known as Slater integrals. They are defined by 


, P(Vp ;r ») 


with 




CB-27) 



TABLE B-l 


COEFFICIENTS OF SLATER INTEGRALS FOR EQUIVALENT ELECTRONS 


f°(i„ V 

_a a 

1 

l 

l 

l 


f 2 (i„ J *> 

_a a 

o 

-2/25 

•2/63 

-4/195 


f4 (i 2 ) 

_flf 8 

0 

0 

-2/63 

-2/143 


f6 ^« V 

_a a 

o 

0 

o 

-100/5577 


COEFFICIENTS OF SLATER INTEGRALS FOR NON-EQUIVALENT ELECTRONS 
(For a tabulated combination (2 , 2.) all non-vanishing coef¬ 
ficients ars listed) p 



0 0 1 - 1/2 

0 110 
0 2 10 

0 3 10 

1 1 1 - 1/6 

12 10 
13 10 

221 - 1/10 

2 3 10 

3 3 1 -!/!<♦ 


1 

2 

3 

4 

5 

6 

8 

8 

8 

8 

8 

8 

0 

0 

0 

0 

0 

0 

-1/6 

0 

0 

0 

0 

0 

0 

-1/10 

0 

0 

0 

0 

0 

0 

-1/14 

0 

0 

0 

0 

-1/15 

0 

0 

0 

0 

-1/15 

0 

-3/70 

0 

0 

0 

0 

-3/70 

0 

-2/63 

0 

0 

0 

-1/35 

0 

-1/35 

0 

0 

-3/70 

0 

-2/105 

0 

-5/231 

0 

0 

-2/105 

0 

-1/77 

0 

-50/3003 
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where r < end r > denote respectively the smaller end the greater of the 
two radius vectors, r^ and r 2> of the two electrons. 

The electrostatic energy for the different terms belonging any 
configuration can now be written as 



where V fc for some of the simpler cases Is as follows: 

(1) All shells completely filled, or there Is a single incompletely 
filled shell with either one electron or one hole. In these cases the 
configuration consists of a single term and, 

V t =0 

(ii) There is a single incompletely filled shell (designated by 
n^ l^) with m electrons 

V t sl)fV,■«,«.) 

Note that the screening parameters are functions of the pertiale 
of the Sleter Integrals F k and G k with respect to the appropriate screened 
nuclear change and F k and G k , which we evaluate using the hydrogenic 
wave function (Equation '-12), are themselves functions of the screened 
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nuclear chargee, and therefore of the screening parameter Jt^’s. Zt Is 
therefore readily seen that the right-hand side of Equation 2-20 Is a non¬ 
linear function of all the screening parameters, it^, » 2 .« N , when 

N Is now Che number of distinct subshells of electrons characterised by 
n and 1 quantum numbers. 


I 

I 
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ft 6 THE SLATER INTEGRALS 


From Equations (B-26), (B-27) and (3-28) it is easy to see that 

tc k 

F and G are symmetric in the two electrons. However, on account of 
the occurrence of r < and t > in this equation, it is necessary to break 
up the integral into two parts, which for the F integral we write as 
follows: 


F k («. WoYlO * 



F* = j r;)dn dn 


(3-32) 


(B-33) 


Note that F^ and F ^ 


are not symmetrical in the two electrons. 


It is easy to prove that 
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(B-34) 


f* (*>,?, >yyO= ff M*-**i f i >^0 

For the G integral, these two parts turn out to be equal and we have 

Jo Jo CB -35) 


The calculations of the Slater F and G integrals, using hydro- 
genic wave functions, is extremely laborious, for electron shells whose 
wave functions possess more than one or two nodes. The author is not 
aware of any published work where these have been calculated as functions 
of the screened nuclear charges of the two electrons Involved, nor are 
any general formulae for calculating them available in literature. We 
have, therefore, developed general formulae which are suitable for hand 
computations as well as computations using digital computers. 


F%i^A-,y,Vi)=Q^ A) 


<S-36) 


-K M > -M 


F- 7 , % A. A* 


( b-37) 


101 


y. 

c^-Ul— 

1 " Cy» 4 Y^"' +a,,1 ~’ 

t- C'' 

rT i 

Anf = (M>*+t4a)i 

•2ZZj. 


(B-38) 


(B-39) 


(B*40) 


(B-41) 


* (>-*>>. t \ (avi, + i«t-\+k-i - 

and ia defined by Equation (B-17). A symbol such as dsnotss 

the lesser of the two quantities inside the bracket. 

I (M1) 


’ (VKf’"*'""'" 


i*j+U+tr Qvant-i-j- 

1 Cf-r-oi 
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♦<*. 


o)tf} 


<B-44) 


<T^ 1,+U-v a-vlM-i^ir (b- 45> 

D - \ -k^U+ir 


For the case of two equivalent electrons F k and G k are Identical 


we have the following simple expressions, 

41 ^ 4*1 <^>9 

where t and D are now given by 

*T = xUa+k+jp 


(B-47) 


(B-48) 


D * ^ ^ \ 


(B-49) 


and ♦ (t,D) Is given by Equation (B-44 ). 

A number of formulas using the above equations are given in 
Table 8-3. Machine computations for t\J*, r^ k ‘ l «*m* Cp have also been 
coexisted for all two electron combinations check marked in Table B-4, 
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TABLE B-3 


SLATER INTEGRALS 

(The formulae listed here are for F* (ap; X^i) and G k (op; X^l), defined 
by Equation (B-50), where - ^a/^. The subscript op are omitted for 
brevity.) 

F°6s,«) = i- [\- (i 4X) r (H3 X 4 l)] 

f*gmp. \[\~ c^rjXx+o] 

ros, 3 d> -l[i - o+x) 7 (^ + ')] 


= i[i - o+x)%u yX x^+ i38$*yW+! 
F*(ip,3<l> ^ [a - (XW4 45V4 !5X4i)J 
F\ jp , 34 >X"[3- 6 +«( 35 o/-373/4151 f 

*+■ 1° gy v 4i7X+3;] 
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TABLE B-3 (continued) 


F(3p,jJ)- -£-[| _ (l-t-y) -t-H^ 

r 

f^p^J) = 1L p -o+xj C}i5~y 1 -\ wx+'wx+JkiX^ 


+-i3oX J -HloV + aiV+a)j 


G°6 v s > a X 3 0+X) 7 (s IS Tk:-Vtv> 


G°6s,3s)= 2 x\ 1+xj ? (5 X*- 3.X + U X-6 2. X 4 ai) 

c 1 ('*»»»- H x 3 (i+xj 7 

G y X* ( 1+ X) 3 (Vfc x - ?H X-> 1-0 
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TABLE B-3 (continued) 


XYl+xTO 01 ^-7oX-f' l f) 

G'(«,}>>> % X J OxK 4 c ^X* ih-sox^i^sX-wX^s) 

a i - 5 

G(is,3d) - is X (t+X) 

G'CMJ- a x 3 (h-xj"(i^hX-4sX+9 / ) 

G*(*/*l> x ~)% w xUs^s isnX-37«x+27) 

G'(if,3s> A x"( i+75"(l X-7°xW6X V -36oXl I 87) 

G°(if,3|>> 3 X S (l+X)"(l«^X- H4E x 4 SIH) 

£Wp> ^f^O+X) (ux^s- y+is) 
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TABLE B-3 (continued) 


\i<> y 3 c»+x5 


G 3 0M> m y 3 (i+^j" 

v 5" —13 v 

1 0+y; 

G*0>vjJ>^- X s (i+7) (azos^-n^X-eUt) 
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The set of simultaneous equations (fe-20) contain partial deri- 
k k 

vatives of F (a,f3) and 6 (a,£) with respect to^ a « well as Al- 

k k 

though F (a, (2) and G (a, ~) are symmetric in the two electrons a and 

their partlals are not symmetric. From the preceeding expressions 
k k 

for F and G the following general properties of these two integrals 
can be derived very simply, 

RVP'AlV = 0 (B-50) 

rTfcftoPi ' V) (MO 

where R stands for either integral F or G and 

= '/fy* = - (*-52) 

The following expression for the partial derivatives are also derived 
easily, 

= vuf' <fr53> 

" ixcw <r ' 54> 

=• P > X«((l ( r -55) 
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It is therefore clear that both partial derivatives of F and 
G integrals can be written as functions of a single parameter X^. 

The case for two equivalent electrons is a particularly simple 
one; the partial of F k (a,a) is merely a number, which is, in fact, the 
value of this integral for ^ = n^"" * t * le ca#e hydrogen atom (with 
atomic number unity) . 

=< ' THE SCREENING PARAMETERS 


We now rewrite Equation (b- 20) explicitly in terms of the 
partial derivatives of the Slater Integrals 






is merely a number, Independent of any of the X^'s, 


56) 

(3-57) 


(B-58) 


(b-59) 


f k ("f" SL) and g k ("l" SL) denote the coefficients in which the SL depen¬ 
dent contributions to the electrostatic energy are explicitly included. 

The sys&ol "i" denotes the orbital quantum numbers of all the incompletely 
filled shells Including the parentages, etc. 
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(B-60) 


1 


We note that, 


%(& - 


;Z-~rcV 


and for the caee Z -* od 


•v -v 00 ^ 

W~* Hi - "mT 


(B-61) 


The method for computing the screening parameters is as follows. 
n B 

We first calculate a set of n's with X. . = . We call these asymptotic 

Q(p n Q, 

n's or x® . With these x® 's, and some suitably choosen large value of 
Z, we compute a new set of X^, and use these to recalculate x^'s by 
iteration. We then decrease Z by successive steps and keep repeating the 
same iterative procedure. Thus we work along an iso-electronic sequence 
and obtain the screening parameters for any desired number of atoms and 
ions belonging to this iso-electronic sequence. Since the procedure is 
very time consuming, a program for computations on a digital computor 
has been developed. 


The asymptotic screening constants for two electron combina¬ 
tions, and using the average energy of a configuration only, can however 
be computed quite easily. They are defined as follows: 

CB-62) 

for a i fi 
and 

•y-'Rrf'O 


ill 


(B-63) 



The asymptotic screening constants for the average energy of a 
configuration can then be constructed from these two electron screening 
constants as follows: 




(B-64) 


B.8 THE CALCULATION OF THE TRANSITION INTEGRALS 


The only case of interest for an electric dipole transition Is 
when the azimuthal quantum number of the active electron changes by one 
unit. Therefore, the calculation of the following Integral is of suf¬ 
ficient generality: 




(B-65) 


Jo 


p- s 


(B-66) 




A similar integral with ^ V i h ** b ** n *v*lu*ted by Gordon^®” 26 ^. 

The necessary formulas for the evaluation of Integrals Involving 
confluent hypergeonetrlc functions have been given, aseng other, by 
Landau and Lifahit*^®" 2 ^ in their mathematical appendix. 
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The p index car be reduced by using the formula 

{ Jr-X "!>?)- j£ (-tf-i, <)} (B. 


67 ) 


whereas the s index can be reduced by using another formula . 

]>(-<- *?)= 

+ S (V-k S+svtOl r ri,t->«X sljr (- «r,-Vu+0 

Successive applications of Equations (B-67) and (B-68) to (b-66) yield 

3jtta ,-Ylx J- —J ~Yl W+ w 0|] 

. Jit w ?>-^[Yi (^ 4 Y* (f ■+ 69) 

jt't can be evaluated in closed form by another formula 
<M> -Ola Y' / —•£»—«?». 

3?K.‘w),c-u * jfaoota) (Yv>0 ( ,. 70) 

- 4 **./(%-&) 

where jFj 1* * hypergeometric function given, In general, by 

c/„ , o.b u o,Ca4i)b(bti) u * 

| t ^ (B-71) 

+ cxCq4 <Xa4a)b(HQ(b4a) 3 . 

c (c+i) (C4J.; i 


113 



The series terminates after a finite number of terms when a and/or b are 
negative integers. On substitution of (B-70) in (B-69) and simplifica¬ 
tion we get -**1 

j&k.-i-d-w wn* la 

•J> (-**')* (v^'-^ + vo(Mrf^ (-<- *) 

where V *V 4 #» 1 

F6*,-M^h-^ p ) 

LlT»-~ JTi -> 


(B-72) 


(B-73) 


The following relation may be used to reduce the number of hypergeometric 
functions from three to two in (B-72) 9 lf desired. 


(c - »■»- - ifc- «•)/'■)£ 60 •+ («■•»)• P-vf, C*i>, 


(B-7*) 


wn«r« 

Jj&O* 

Combining (B-65) with (B-72) we get 






-aCM» 

+(*,W*--v1fO 


j(Y,V 




(B-75) 


(-<-0 

Ox 4 #) 
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which for th« apodal caao^pj^K* raducaa to tha following 
aquation flrat darivad by Gordon® 

lb - - 


w.-Wt-af-a. 

cn,-*o 


l |f,t ^-,-^p 


te-76> 


-&£)* 

apodal caaaa can ba obtainad aaallv 


VI-3 


The following apodal caaaa can ba obtainad aaally 

It( w P>' s .ViV*-) (rf-0 'Jff K (iy L -ny,) 3 cb- 

I. en>» s '>Y,YO •^(w'-oYvf 


•77) 


J1-H 


4 < ^*'TfiV ~ ? 

I.K*p £ 




m-h 


(V^>.> hn 


T <( a r>"‘ *,*%.) [***•*?£ 


(B-78) 


CB-79) 


(B-80) 


- 3 

cv,->,r 

o.-*o w4H 
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B.9 ALTERNATE EXPRESSION FOR TRANSITION INTEGRALS 


For the uk« of generality we shall evaluate the following 
integral, 50 


tj M * j fMiVS) 

‘ .r r +*i-Vl <* r 


(B-81) 


For Coulomb potential wave functions, we shall later put q equal to aero, 
p equals unity for the electric dipole case, and it equals 2, 3, 4, etc. 
for higher oultipole cases. Substituting from Equation (3-12) we get. 


(B-82) 


. jf ( F, 

This integration can be performed using certain standard results In the 
theory of confluent hypergeometrlc functions, (see for example Slater ^ 8 28 \ 
page 54). Performing the integration we obtain,. 


pmge J**/ • rarfond.ua uw UH|[IIUW wm wwtu,. . j 

12 > -9^ fr, 

• (B-83)* 


*A similar result has also been obtained by Prof. D. H. Mensel (unpublished). 
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where the expression in the square bracket following the sunsution sign 
can be readily seen to be equal to d* . 2 ?j Is the hypergeoaetrlc func¬ 

tion defined by Equation (B-71). The conputation of the expression on 
the right hand side of the above equation is quite straight forward in 
principle, but becones very laborious as n| and nj become large. 


We shall now proceed on to the apeclal case we are interested 


in, newly the case where p ■ 1 and q ■ 0. For this case, 

t F, ^ 2-^j+T- 

gt*r ,i On?)'. (»W 

-Z-i (?*-£> i y. ~ 

# T\ Ivv.) 

= li 4 {«■*«.+ /»*«• V 

^ -Oi 

l,(«,W»',y,Vv)= 

»'•* >!" x 

(W0* ,+fv,H 


(B-84) 


(B-85) 


X", £. £ i ■»! 

V V? 

■(VVO^ 


<B-#6) 
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Neediest to say that tha same expressions can also ba obtained by integra¬ 
tion of Equation (b- 81) with the radial wave functions given by Equation 
(B-12) combined with Equation( b- 14). Soma further simplifications of the 
above expression can be brought about as follows. The common detMMd.nator 
of all the terM under the two suaomtlons is <% +%)*? + Taking 
the connon denominator and performing binomial expansion of the resulting 
torn, naaaly (fc + n 5 "J, ua obtain 

Jjf , 

" dfii & 


i-? 






w jw iwr-vwi 


• (**+*£- /V\ 

? ’w-mM-i-w >». ly, ) 




(B-87) 


Replacing p by another index u, such that u ■ p + J, we get, 

), *** f 

~ i21 (-£-) 

1 #, sw-i-u)«(u-i)rT/ 




„ -v^f+Mr-w v « 

-I 


U=o 




(B-M) 
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,th4)='L 


^ *■ i 

-T- Z+ 4 , 4 ? 2 . 


.(U.+v*-m)'. ( 2 !dL^zlii 24 —— 


u takes on all Integral values from 0 to (n| + nj). \|r is a function of 
the quantum numbers of the two electrons under consideration, but it is 
Independent of the screened nuclear charges,^ and^. When n* and 
n* are so large that desk calculations are too tedlus, the calculation 
of \^ u on electronic computors is a trivial matter, giving us analyti- 
cal expressions of 1^ in terms of as follows, 

I, K^vUv.vo« Q>Qx i t,+ ^ 

. *• ® v (B-90) 


<&**> U,._ , . 


Kl.% «?"-«■ , <* 




electrons, t 1* not. 


Table 2-5 gives expressions for 1^ for a selected number of 


transitions. 


TABLE B-5 DIPOLE IMT8GRALS 


Ijfcp.ta) = 

bMifft, 

(fy 

- 

I,(3f>,U) . 

(Aiifa££. 

-( z yh) 

I,(4 f ,U) - 

w 

i 

• 

O^C» 

1 

I/Sf.U) - 



X<6f-U) - 

*'dih£ 

i 

vb 

1 
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TAKE B-5 (eont.) 


1 

I 

I 

I 

I 


5/a. 1/v 

IM - ( i±g t [-3y,+' 5 y.k-*i] 

I to' 2 *) * [4^ *&»*»»*-#] 


LM - 

im* - “ ggy 


<**»> 


kw - “ffy W 
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TAK.E B-5 (cont.) 


IM - 


I 

I IM ■- 

I 

I K2f A) - 

I(3f,5-) - 


■* 9 I‘ 


-»>V 


uFiw? jjvtyXV 
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TABLEB-5 (cont.) 


IM 5 —{(^ - 

+ 2i4(6^‘j v S 

l&ty m a ^y [(r»i 

-W(jf-»*»V "*}•$'tffy'bf 

t224(/1 <*& •*$£»»*$*♦$ 
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TABLE B-5 (cent.) 


- t3#j, t ^* jfXj^T 
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TAB.EB-5 (cont.) 






vj0(w) - |+7e w] 
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V&4 - 


i 

I 

! 


TABLE B -5 (cont.) 


• - f (ty-wwi 

- *($' 4S »V f°&$+ so £)( Vb'f 

" s b$ 

*M s t$ 

-**(*$£ 
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l(6fM * 





TAK.E b- 5 (eont.) 


***$&***$■$'**$$" 
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l(3Uf) * 


llMfa) - 
K&Uf) - 
K&u-p) - 
i(7«Uf) • 

M - 


l(MM) ‘ 


TAK.E B-5 (cont.) 
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- 

Xfaty - 

- 

I(6J,4f) - 


Vz r ^ X 31 ^ ~\ 

2SC17) V' * • f to-y - 5-9V HfiK- l5 y»] 

(V%y v." 

(VYO' 

Z56 (a>? vlj /v c v5\^V' c,7 ^Yx- + ^'^ r ' 5 ^ 

~ly^W r ~ 


W l 




***$)(}**) 


IfrSf) 




TABLE B-5 (cont.) 
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TJULEB-5 (cont.) 


K 5 ^f) * — (H) k ‘ ‘tiytxwf 

'W&'Vby'tostfiy if 

»,*>■ (»•>;<»») 

«*.«■ *3$/ fo-») 

i *^ L -(i>-»(‘j- < i-) 


I«.u) 
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TABLE B-5 (cont.) 


IM- 1 £3 ^- CYrn^ 

X(6f»44) - ^T-V' 1 Vi 

|l ifi 

im- 3 %%P Cv -- Y ' ) 
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w 



um - 


B-5 (cont.) 


-^ M *V' 4 atfXa-*) 

♦W*] 


BIO Results of Computations 


(a) The Normalization Factor, Q n£ 


Q n ^ is defined by Equation (B-13) 



1_ 

tti+i): 


(B-13.) 


The following recursion relations are useful in extending a 
table of q nl : 

I 1/2 

y n,i-i 


q t Q ' ^ (B- 91 ) 


2l(2i+l) 


r n 1/2 

Q , ail mill Q 

H nt |_ n n-i J w n-l,i 

in/ ’ [*? iTTziry Vi. 


A few special cases of are given below: 

*" 1/2 


<*n 


«n,0 " 2 


i-1 

(B- 93 ) 


(B-94) 


<n.n-l * < 2n > ! 


• " 1/2 


(B-96), 


A sl^>le procedure in calculating is to start with Q nQ which 
• 1/2 

is always 2 , and successively calculate Q^, Q^, using Equation fB-91 > 

until one arrives at n- j. Equation (B-96) can be used as a final check 

* The work performed for this section was sponsored by both this contract 
and Contract Ho. AF 33(600)-4S6S7 


134 


The value* of for n up to 9 and all allowed values of l were 
computed in the above manner, and subsequently machine computations 
were also performed. Table B-6 lists the Q^. 

(b) The Coefficients d^ 

The d*^'s are the coefficients occurring in the hypergeometrlc 
function jPj f-(n-i-l), 2i+2;x\ . They are given by Bquation (B-15). 


d D* ■ ( ' l > ! ofeHt 1 ! 

(B-15) 

For i > n-X-1, d^ vanishes. 


The following recursion relations are useful 
the table of d^: 

in extending 

.i (n-i-i) .1-1 

d ni * (21+l+i)i d ni 

(B-97) 

d i . _ d i 

ni (n-J^ (22+0(21+1+1) d n,i-l 

(B-98) 


(B-99) 

d nX * (2i+i)Ki+l+i) d n-l,l-l 

(B-100) 

A few special cases of d^ are of interest 


d ni " 1 

(B-101) 
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.2041241+00 


.9354141+00 


.3333338+00 8d .3622848+00 

.3726788-01 8f .6397088-01 

.4564358+00 8g .6155598-02 

.7905698-01 8h .3494708-03 

.4980128-02 81 .1185408-04 

.5773508+00 8k .2186208-06 

.1322888+00 9p .1054098+01 

.1259888-01 9d .4624818+00 

.5249508-03 9f .9343538-01 

.6972168+00 9g .1046258-01 

.1972038+00 9h .7117658-03 

.2439758-01 91 .3060688-04 

.1515408-02 9k .8244688-06 

.4569118-04 91 .1249778-07 

.8164968+00 
.2738618+00 
.4123938-01 
.3290308-02 
.1465378-03 
.3386858-05 
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(B-102) 


n-l-1 

' 2(4+1) 


d 


1 

n,0 


(- 1) 1 



(B-103) 


A table of for values of n up to 5 was prepared using the above 
recusion relations and subsequently machine computations were performed 
for n up to 9. These are given in Table B-7 

(c) The Screening Parameters 

Equations (B-56) define a system of simultaneous non-linear 
equations in which the right hand side is a function of all the screening 
parameters of a given problem. A Fortran program for use 

on digital computers has been developed for solving these equations. 

We work along an isoelectronic sequence. We use as our starting point 
the case when Z -* oo . The screening parameters for this case can be 
computed easily vith given by Equation (B-61). As an intermediate 
step we perform iterative calculations for Z = Z q + 40, where Z Q is the 
lowest member of the isoelectronic sequence, using n Q 's for Z * CD case 
as our starting point. The Iterations are carried until the relative 
difference between the results of two successive iteration is less than 
a prescribed quantity e (for computations reported here € * 0.0001). 

The process is repeated successively for the cases Z « Z q + 19, 2^ + 18, 
....Z Q . The screening parameters, expressed as s(a)(* n Q ) are given 
in Tables B-8 to B-19. 
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.7142861-01 -.2976191-02 
.4166671-01 -.1388891-02 
.2727271-01 -.7575761-03 
.1923061-01 -.4578751-03 







TABLE B-7(continued) 
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Ion 

TABLE B-8 

SCREENING PARAMETERS 

Be I Iaoelectronic Sequence 

Configuration: ls*2a^ 
Term: l S 

s(ls) 

a(2a) 

Li* 

3.2761 E-l 


1.7902 E+0 

Be I 

3.1861 E-l 


1.8599 E+0 

B II 

3.1369 E-l 


1.8917 E+0 

C III 

3.1184 E-l 


1.9025 B+0 

N IV 

3.1164 E-l 


1.9036 E+0 

0 V 

3.1221 E-l 


1.9006 E+0 

F VI 

3.1312 E-l 


1.8961 EfO 

Ne VII 

3.1414 E-l 


1.8912 E+0 

Na VIII 

3.1517 E-l 


1.8863 EfO 

Mg IX 

3.1616 E-l 


1.8817 EfO 

A1 X 

3.1708 E-l 


1.8775 E+0 

Si XI 

3.1794 E-l 


1.8736 B+0 

P XII 

3.1873 E-l 


1.8700 E+0 

S XIII 

3.1946 E-l 


1.8668 E+0 

Cl XIV 

3.2013 E-l 


1.8639 B+0 

A XV 

3.2074 E-l 


1.8612 B+0 

K XVI 

3.2130 E-l 


1.8587 B+0 

Ca XVII 

3.2182 E-l 


1.8565 B+0 

Sc XVIII 

3.2230 E-l 


1.8544 1+0 

T1 XIX 

3.2274 E-l 


1.8525 B+0 

Z ■ oo 

3.3353 E-l 


1.8078 B+0 
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TABLE B-9 

SCKEMXM6 PARAMETBIS 
Be I Iaoelectronic Sequence 
2 

Conflguretion: Is 2s2p 
Tens: *P 


Ion s(le) s(2s) s(2p) 


Li' 

3.1832 

1-1 

1.5178 

BfO 

3.0030 

BfO 

Be 1 

3.2321 

B-l 

2.1229 

BfO 

2.0640 

BfO 

B II 

3.2444 

1-1 

2.1965 

BfO 

1.9749 

BfO 

C III 

3.2590 

E-l 

2.2165 

BfO 

1.9390 

BfO 

N I? 

3.2766 

E-l 

2.2186 

BfO 

1.9180 

BfO 

0 V 

3.2949 

E-l 

2.2138 

E+O 

1.9037 

BfO 

P VI 

3.3126 

B-l 

2.2065 

BfO 

1.8932 

BfO 

Ne VII 

3.3289 

E-l 

2.1987 

BfO 

1.8849 

BfO 

He VIII 

3.3437 

E-l 

2.1911 

BfO 

1.8783 

BfO 

Mg DC 

3.3570 

E-l 

2.1840 

BfO 

1.8728 

BfO 

A1 X 

3.3690 

E-l 

2.1775 

BfO 

1.8682 

BfO 

81 XI 

3.3797 

E-l 

2.1715 

BfO 

1.8642 

BfO 

P XII 

3.3893 

E-l 

2.1662 

BfO 

1.8608 

BfO 

8 XIII 

3.3980 

E-l 

2.1614 

BfO 

1.8578 

BfO 

Cl XIV 

3.4059 

1-1 

2.1570 

BfO 

1.8551 

BfO 

A XV 

3.4130 

E-l 

2.1530 

BfO 

1.8528 

BfO 

X XVI 

3.4195 

1-1 

2.1494 

BfO 

1.8507 

BfO 

Ca XVII 

3.4254 

B-l 

2.1461 

BfO 

1.8488 

BfO 

Sc 8VIII 

3.4308 

E-l 

2.1430 

BfO 

1.8471 

BfO 

Ti XIX 

3.4358 

1-1 

2.1402 

BfO 

1.8456 

BfO 

8 

■ 

N 

3.5507 

1-1 

2.0773 

BfO 

1.8129 

BfO 
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TABLE B-iO 
SCREEN DIG PARAMETERS 
Ba X Xaoalectronlc Saquanca 
2 

Configuration: It 2t 3p 
Tara: 1 P 


Ion 

•<10 

»(2p) 

• <3p) 

Li” 

3.1396 E-l 

1.9799 BfO 

2.5502 BfO 

Ba I 

3.1969 E-l 

1.8484 BfO 

2.7110 BfO 

B II 

3.2529 E-l 

1.7602 BfO 

2.7775 BfO 

C III 

3.2957 E-l 

1.7074 BfO 

2.7953 BfO 

N IV 

3.3280 E-l 

1.6774 BfO 

2.7864 BfO 

0 V 

3.3527 E-l 

1.6598 BfO 

2.7677 BfO 

F VI 

3.3720 E-l 

1.6490 BfO 

2.7471 BfO 

Na VII 

3.3876 E-l 

1.6419 BfO 

2.7276 BfO 

Na VIII 

3.4003 E-l 

1.6370 BfO 

2.7102 BfO 

Mg IX 

3.4110 E-l 

1.6334 BfO 

2.6950 BfO 

Al X 

3.4200 E-l 

1.6307 BfO 

2.6818 BfO 

Si XI 

3.4276 E-l 

1.6285 BfO 

2.6703 BfO 

P XII 

3.4343 E-l 

1.6268 BfO 

2.6603 BfO 

S XIII 

3.4401 E-l 

1.6253 BfO 

2.6514 BfO 

Cl XIV 

3.4453 E-l 

1.6240 BfO 

2.6437 BfO 

A XV 

3.4498 E-l 

1.6230 BfO 

2.6368 BfO 

K XVI 

3.4539 E-l 

1.6221 BfO 

2.6306 EfO 

Ca XVII 

3.4575 E-l 

1.6212 BfO 

2.6251 BfO 

Sc XVIII 

3.4609 E-l 

1.6205 BfO 

2.6201 BfO 

Ti XIX 

3.4639 E-l 

1.6199 BfO 

2.6156 BfO 

Z - 00 

3.5261 E-l 

1.6072 BfO 

2.5266 BfO 
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tMLI B-ll 

SCtESflMG PABAMETEBB 
0 1 Iaoalactronlc Saquance 
2 2 4 

Configuration: la 2a 2p 
Tara: *8 


Ion 

■ (!•) 

a(2a) 

a(2p) 

■" 

3.3288 E-l 

2.9181 BfO 

3.6713 BfO 

0 X 

3.4078 E-l 

3.0238 BfO 

3.5794 BfO 

P II 

3.4825 E-l 

3.0779 BfO 

3.5248 BfO 

Na III 

3.5524 E-l 

3.1080 BfO 

3.4873 BfO 

Na IV 

3.6170 E-l 

3.1253 BfO 

3.4595 BfO 

Mg V 

3.6763 E-l 

3.1354 BfO 

3.4378 BfO 

Al VI 

3.7303 E-l 

3.1411 BfO 

3.4204 BfO 

SI VII 

3.7784 E-l 

3.1442 BfO 

3.4061 BfO 

P VIII 

3.8241 E-l 

3.1455 BfO 

3.3940 BfO 

S IX 

3.8648 E-l 

3.1458 BfO 

3.3837 BfO 

21 X 

3.9019 E-l 

3.1453 BfO 

3.3748 BfO 

A XI 

3.9358 E-l 

3.1445 BfO 

3.3671 BfO 

K XII 

3.9668 E-l 

3.1433 BfO 

3.3602 BfO 

Ca XIII 

3.9953 E-l 

3.1420 BfO 

3.3541 BfO 

Sc XIV 

4.0216 E-l 

3.1406 BfO 

3.3487 BfO 

T1 XV 

4.0458 E-l 

3.1392 BfO 

3.3438 BfO 

V XVI 

4.0682 E-l 

3.1377 BfO 

3.3394 BfO 

Cr XVII 

4.0890 B-l 

3.1362 IfO 

3.3353 BfO 

Mn XVIII 

4.1083 E-l 

3.1348 BfO 

3.3316 BfO 

Pa XIX 

4.1263 E-l 

3.1334 BfO 

3.3283 BfO 

Z - 00 

4.6176 E-l 

3.0813 BfO 

3.2504 BfO 
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The hydrogenlc wave functions using the screening parameters 
are compared with the self-consistent field wave functions for neutral 
oxygen (Hartree, Hartree and Swirles^ 8 "^) and negative ion of lithium 
(Roothaan, Sachs, and Weiss ^ B "^) in Figures B-l to B-4. The hydro¬ 
genlc wave functions are given below: 

0 I <ls 2 2s 2 2p 4 X S) Is electron: 42.394 r exp (-7.6592 r) 

0 I (la 2 2s 2 2p 4 l S) 2s electron: 7.8491 r exp (-2.4881r) (l-2.4881r) 

0 I <ls 2 2s 2 2p 4 1 S) 2p electron: 8.3864 r 2 exp (-2.2103 r) 

Li (1s 2 2s 2 ) 2s electron: 0.94079 r exp (-0.66524r) (l-0.66524r) 

The agreement between the two sets of wave functions is 
remarkable for 0 I and is good for Li. 

Total ionization energies based on hydrogenlc wave functions 
are compared with experimental values in Table B-20. The agreement is 
excellent; for Be I it is as good as that for the multiconfiguration 
Hartree-Fock wave functions of Kibartas, Kavetskis and Iutsisf® 32 ^ 
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Figure 1-4. Li” (le 2* ): Wave Function of 2s Electron 


TABLE B-20 


TOTAL IONIZATION ENERGIES (Atomic Unit) 



Hydrogenlc 

Wave Function 

H-F 

Have Function 

Experimental 

0 I ^s) 

- 74.62 


- 74.95 

F II ^S) 

- 98.61 


- 98.92 

Be I ( l S) 

- 14.70 

- 14.64* 

- 14.67 

B II ( l 8) 

- 24.38 


- 24.35 

C III ( 1 8) 

- 36.55 


- 36.55 

N IV ( 1 S) 

- 51.23 


- 51.25 

0 V (*8) 

- 68.40 


- 68.45 


♦Baaed upon multiconfiguration Hertree-Foek calculations of Kibartaa, 
Kavatakla and Iutala> B ’ 32 ^ 





(d) The Dipole Integrals 


The electric dipole Integral 1^ Is given by Equation (B-90) 

In terms of a coefflcientY^n^i n^) defined by Equation (B-89). 
Computations on an IBM 7090 have been performed for tha coefficients 
Y U for all two-electron combinations of practical interest, and the 
results are given In Table B-21 to B-26. In order to reduce human 
errors in the reproduction of a large collection of numbers, the print¬ 
out from the digital computor has been photocopied. P8I 7, for example, 
means'^ 7 . Numbers such as 0.23999999 or 0.24000011 should be rounded 
off to 0.24. 

2 

A few sample values of the transition Integral 7 computed 
with the screened hydrogenic wave function are given In Table B-27, 
and compared with calculations based on Coulomb approximation 
and H-F wave functions. The agreement is setlsfeetory. 


(^»*« I ») n> j v ra-t 
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TABLE *-ll ¥ (OjS.t^p) (continue) 
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HjP) (continued) 
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TABLE B-23 ¥ (n^.^d) (continued) 
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(4,7) (continued) 
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TABLE B-27 


TRANSITION INTEGRAL cr 2 


Atom 

2 

Transition Integral or 

Hydrogenic 

Coulomb 
Approx. 

Hartree-Foek 

Be I 

2.27 

2.03 

1.86 

B II 

0.977 

1.010 


C III 

0.544 

0.573 


N IV 

0.348 

0.368 


P XII 

0.0522 

0.0540 
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APPENDIX C 
D. Layzer 


NEW GENERALIZATION OF THE SCREENING THEORY 
AND THE HARTREE-FOCK THEORY 

C .l. Introduction 


The distinguishing feature of the Z-dependent theory of 
atomic spectra (Layzer/ C * ) henceforth referred to as A) 
is the fact that it allows for certain effects of configu¬ 
ration mixing in the first approximation. This ensures that 
the predicted term energies have the correct linear depen¬ 
dence on Z. Experiment shows that the structure of term 
spectra, in both simple and complex systems,, depends strong¬ 
ly on the linear Z-dependence of the energies. A knowledge 
of the energy coefficients W 2 and in the expansions 

w - w 2 z 2 + w x z + W Q + o(z _1 ) (C-l) 

enables one to describe quantitatively the most conspicuous 
features of term spectra in isoeleotronic sequences (see 
especially Figures 1 and 2 in reference C-l). 

In order to calculate the coefficients (W 2 is given 
by the well-known hydrogenic formula) one needs to know 
nothing about the radial dependence of the atomic charge 
distribution beyond what is provided by the theory of the 
hydrogen atom. That is, all the necessary calculations 
utilize hydrogenic wavefunctions. One can, of course, use 

less simple wavefunctions in these calculations, as was 

(C- i) 

done in the first calculations of this type (Layzei ' 
but such refinements only improve the accuracy of higher 
terms in the expansion. 
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For some purposes, however, one does require a reason¬ 
ably accurate knowledge of the radial charge distribution 
and hence of the radial wavefunctions. One needs this in¬ 
formation, for example, in order to calculate electric- 
dipole transition probabilities. 

The various procedures in common use for calculating 
accurate radial wavefunctions are all based on the varia¬ 
tion principle. They differ through the restrictions they 
impose on the form of the variational functions. Hartree- 
Fock wavefunctions (i.e., exact solutions of the Hartree- 
Fock integro-differential equations) fall at one end of 
the gamut of variational wavefunctions, screened hydro- 
genic wavefunctions at the other. The latter involve only 
a single adjustable parameter, while the former are re¬ 
stricted beforehand only by the usual conditions of inte- 
grability and smoothness. Other types of variational wave- 
functions seek to achieve some of the accuracy of Hartree- 
Fock wavefunctions without sacrificing all the simplicity 
of screened hydrogenic wavefunctions. Which type of varia¬ 
tional wavefunction one chooses to employ will be governed 
by the purpose in view. The Ilartree-Fock description is 
certainly the most accurate; from a qualitative standpoint, 
however, the screening description offers unrivalled ad¬ 
vantages of simplicity and flexibility. 

Any meaningful comparison between different variational 
wavefunctions presupposes a common approximation scheme. 
Thus we may compare wavefunctions calculated with neglect 
of configuration mixing, or with certain specified kinds 
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of configuration mixing taken into account. (In practice, 
the calculations that allow most fully for configuration 
mixing are those based on the simplest wavefunctions, so 
that strict comparisons are difficult.) 

According to the Z-dependent theory, the most important 
kind of configuration mixing is that between configurations 
belonging to the same complex (i.e., the same set of princi¬ 
pal quantum numbers and the same parity). On the other hand, 
the theory of radial wavefunctions has usually been devel¬ 
oped within the framework of the central-field approximation, 
though attempts have occasionally been made to allow for 
effects of configuration mixing. A notable example is the 
work of Hartree, Hartree, and Swirles^ C_3 \ which dealt 
with the mixing of configurations ls 2 2s 2 2p n with configura¬ 
tions ls 2 2p n+2 . However, such attempts have lacked an ade¬ 
quate theoretical framework, comparable in simplicity and 
generality with that provided by the central-field approxi¬ 
mation. The present paper seeks to supply such a framework. 

C .2. The Theoretical Framework 

Let Y denote an antisymmetric N-electron wavefunction 
pertaining to a definite state in a definite electronic con¬ 
figuration. The quantity 

W = <YHY> , (C-2) 

where H is the N-electron Hamiltonian and the brackets < > 
indicate an integration over all 3N coordinates and a sum¬ 
mation over all the spin variables, is the variational ener- 
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gy. It is a functional of the radial wavefunctions. If 
these are specified beforehand, except for the values of 
certain parameters, then W reduces to a function of these 
parameters. W exceeds the true energy of the state under 
consideration provided that f is normalized and is orthogo¬ 
nal to the true eigenfunctions of all states lower than the 
one under consideration. The condition 6W * 0 leads to a 
set of simultaneous integro-differential equations for the 
radial wavefunctions, or to a set of simultaneous algebraic 
equations for the adjustable parameters. 

One can easily generalize the theory just sketched so 
as to include the effects of mixing among a finite number 
of configurations. Let ... , denote variational 

wavefunctions for states with fixed values of the quantum 
numbers SLp in different configurations. Suppose, for the 
sake of definiteness, that ^ corresponds to the deepest of 
the states under consideration. Then the corresponding 
variational energy W 1 is the lowest eigenvalue of the matrix 
|| || . As before, the condition 6W 1 » 0 yields a 

set of simultaneous equations for the unknown functions or 
parameters. (Note that all the radial wavefunctions or 
parameters that figure in the complete set may 

occur in a single equation.) 

An accurate solution of the secular equation for W* 
would lead to a set of variational equations of prohibitive 
conqplexity. The procedure usually adopted consists in 
retaining only first- and second-order contributions to W*, 
in the sense of ordinary perturbation theory. Thus W 1 is 
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given by 


* <¥ x |H|y 


Z 

i/l 


K^IhI^I 2 


(C-3) 


where denotes the mean energy of the ith configuration. 
Even with this simplification, the variational equations with 
configuration mixing taken into account present considerable 
more difficulty than the usual equations. And frequently the 
effects of mixing are so large that second-order perturbation 
theory does not describe them adequately. 


Another important difficulty inherent in the usual methods 
of dealing with configuration mixing concerns the choice of 
interacting configurations. It is of course true that the 
more interactions one takes into account, the more accurate 
will be the theoretical description. It has been amply dem¬ 
onstrated, however, that in order to obtain predictions of 
term spacing in neutral and slightly ionized systems that 
are accurate to within a few percent, one needs to allow for 
a very large number of interactions. Hence a calculation 
that allows for mixing between a small number of configura¬ 
tions would be significant only if the choice of the configu¬ 
rations considered could be justified on grounds other than 
mathematical convenience. 


Hie problem of configuration mixing takes on a differ¬ 
ent complexion if we adopt the Z-dependent approach. Con¬ 
sider the expansion 
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(C-4) 


*■ I Vi • 

i 

where the refer to distinct configurations and the sum 
runs over all configurations connected with configuration 1 
(say) by nonvanishing energy matrix elements, (we assume 
that the state corresponding to ¥ has been nominally assigned 
to configuration 1). 

The numerical values of the coefficients a^ depend on 
the choice of radial wavefunctions. For example, if the 
radial wavefunctions in ^ are solutions of the Hartree- 
Fock equations, then the coefficients pertaining to con¬ 
figurations with the same set of azimuthal quantum numbers 
as configuration 1 will all be very small, because the vari¬ 
ation principle ensures that most of the effects of mixing 
between configurations characterized by the same set of azi¬ 
muthal quantum numbers have already been incorporated in the 
wavefunctions. Now, every coefficient a. is either 0(1) or 

-1 1 
0(Z ). The coefficients of the first type are associated 

with configurations that have the same set of principal 

quantum numbers as configuration 1; those of the second type 

with configurations that have a different set of principal 

quantum numbers. 

The coefficients of the first type can be written in 
the form 

a. = a° + 0(Z _1 ) , (C-5) 

0 

where a^ is independent of Z. It was shown in A that the 
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constants a? do not depend on the form of the radial wave - 
functions . The coefficients and parts of coefficients that 
do depend on the form of the radial wavefunctions are 0(z”*). 

We may conveniently split off the terms on the right side 
of ( C-4) that have constant coefficients: 

*-£•!*!♦ I*i’i • ,c -‘> 

i i 

where 

a^ = a i - a° = OU** 1 ) . (C-7) 

The nonvanishing coefficients a° are finite in number. They 
can be calculated with arbitrary precision, being essentially 
by-products of the calculation that yields the energy coeffi¬ 
cients (see A). 

We propose to replace the usual approximation, which con¬ 
sists in neglecting configuration interaction entirely, by an 
approximation that represents the wavefunction for a given 
state by the finite sum 


T - £ . (C-»» 

where the coefficients may be regarded as Known numbers. (Por 
spectroscopic terms belonging to completes of the type (1 2 2 <1 ) 
and (l 2 2 q )° they are given in Table C-l. 


The approximation just described becomes increasingly 
accurate as z increases along an isoelectronic sequence. 
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Even for moderately ionized and un-ionized systems, however, 
it should prove superior to the central-field approximation, 
since experiment has shown that the linear Z-dependence of 
the term energies usually dominates the structure of the term 
spectrum. 

C. 3. The Generalized Screening Parameters 

Notation. The letters r,t,t* will denote configu¬ 

rations; the symbol r will distinguish different states in the 
same configuration; a will denote an orbital nl; and a will 
denote a complex. Thus the label a(r)T distinguishes the 
state T nominally assigned to the configuration r of the com¬ 
plex a. In the present notation, equation ( C-4) takes the 
form 


T(arr) * £<arr|tr>¥(tr) . (C-9) 

t 

The Hamiltonian H for a many-electron system is given, 
in the nonrelativistic approximation, by 

H » » 0 + V , (C-10) 

where 


H 0 = £ (-p?/2m - ze 2 /r i ) , 


(C-ll) 


V 


F* 

i<j 


(C-1Z) 
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The variational energy in the present approximation is 
defined as 

W = <arr|n|arr> . (C-13) 

We now evaluate W using screened hydrogenic wavefunctions 
for the undetermined radial wavefunctions that figure in 
this formula. 

For the one-electron part of W we obtain 

<arr|H 0 |arr> = £ <arr| tr> <tr|H Q | t ' P> < t T| arr> 
t,t * 

- £ l<arr|tr >| 2 <tr|H 0 |tr> 

t 

- I l<arr|tl '>| 2 t£ qat (-Z 2 + S 2 t )/2n 2 ) , 

t a (C-14) 

where denotes the number of orbitals of type u in the 
configuration t, and s^ t denotes the screening parameter 
that corresponds to the orbital a in configuration t. A 
more complete notation for this screening parameter is 
s z (at|rD. The suffix Z indicates that the screening para¬ 
meter depends in general on Z; s z differs from by terms 
0(Z S . The labels rl’ indicate that the value of the screen¬ 
ing parameter depends on the state under consideration. 


For the two-electron contribution to W we have 


1 


<arr|v|arr> = £<arr|tr> <tr'|v|t *r> <tT|arr> . 

t,t‘ (C-15) 

Hence 

<»rr|H|arI> . £|<arr|tr>| 2 + a 2 t /2n 2 ) + 

t a 


The variational equations are 


I<-S? <tr|v|t ^ ) • 


(C-16) 


3T- ■ - sr; «tr| v| tr> ♦ 2 £ <= t f,rr < tr l v l t ' I » - 0 • 

at n a at tyt (C-17) 

where 

V Z - s at - c tf,rr - <tT|.rr>/<tr|«0 . 

(C-18) 

The diagonal matrix element <tr|v|tr^ may be written in 
the form 


<tr|v|tO - v t + v tr , <c-w) 

where the contribution V represents the average value of 

<tr|v|tr>. Explicitly, 


V 


t 


^<taSL|v|taSL>(2S + 1) (2L ♦ 1) 
aSL 

£ (2S + 1)(2L + 1) 
aSL 


(C-20) 
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The term-dependent contribution v^p vanishes when the con¬ 
figuration t consists of complete groups only, or of complete 
groups plus or minus a single electron. The contribution to 
v fc p from interactions between electrons in a group l n satis¬ 
fies the symmetry condition 

v(l n D = v(l 4l+2 " n r) . (C-21) 


In a first approximation one may wish to neglect the term- 
dependence of the screening parameters. We accordingly define 


n 2 3v 

•toll V}'. 5 ' • (C-22) 

tx at 

If the configuration t contains no more than one incomplete 
group, can be put in the form 


where 


v- Iv ap) ' 

<*P 


rq a\ (a = p) 

12 ^ 



(a ji p) 


(C-23) 


(C-24) 


The two-electron interaction (ap) is defined in A; in a first 
approximation (ap) ~ F^(ap) . The coefficient denotes the 

number of ways in which two electrons labelled by the indicated 
quantum numbers can be selcted from the configuration t. 


With the help of these definitions we can write s(at) in 


the form 


s(at) 


(C-25) 


” I % ‘ V* ( “ |P> ' 

P 

where 

8 (a|P) = (1 - £ 6 a p) n a §7 (<*P) • (C-26) 

The two-electron screening parameter s(a|p) provides an unam¬ 
biguous , qunatitative measure of the screening of an a elec¬ 
tron by a p electron. 

The screening parameters s (a|p) are independent of the 
configuration t, and can therefore be evaluated once and for 
all. 


The exact screening parameters are given by 
2 


s(at|rD 


■ tat > + f St" [v tr 


2 X c tt , ;rr v tt , r 

tyt 


) , 

(C-27) 


where 


v tt r " < tr lv|t*r> • 


(C-28) 


In this equation one is not in general justified in neglecting 
the contribution from configuration mixing. In some instances the 
second term in the brackets ( ) may be considerably larger than 
the first term. This is always the case, for example, when the 

configuration t is composed of complete groups but contains one 

2 2 

or more incomplete shells, e.g.. Is 2s . Then v^ - 0, but v tt ,j, #* 0. 
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The set of equations ( C-2?) for all orbitals a in the 
mixture ( C-9) determines the corresponding screening para¬ 
meters s z (ot|rr). In spite of their complicated appearance, 
these equations are comparatively easy to solve numerically. 

I.. practice one uses an iterative procedure, first evaluating 
the right sides for Z ■ «. This gives the screening constants 
S M (at|rO. In the next approximation one evaluates the right 
sides for finite values of Z, replacing the parameters s by 

the constants s , and so on. 

00 

One can take advantage of the dominant contribution of 
s(at) to s(at|rr) to shorten the iterative evaluation of Z- 
dependent screening parameters. One begins by solving the 
system of equations (C-22) for s z (at). One then evaluates 
the term-dependent correction to s z (at|rO using s z <at) in 
place of s z (at|rO in equation ( C-27). With this method 
one needs to go through the iterative procedure only once 
for a given configuration. 


C. i»—8ffl£.r a Used partyee-Fock yhegpy 


The Hartree-Fock equations may be written in the form 


( - * ' 


Ul + 1) _ Z ~ Y nl 
2r 2 


)P„ 


(C-29) 


The quantities and are functionals of all the radial 
wavefunctions appropriate to the configuration being consid¬ 
ered. The exchange term X ftl has its origin in the antisymmetry 
of many-electron wavefunctions. 


If we take as our starting point the approximation ( c-9) 
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instead of the central-field approximation, we arrive at an 
equation that is formally identical with ( C-29) This equa¬ 
tion differs, in fact, from the ordinary Hartree-Fock equation 

only through the exchange term X ,, which contains additional 

n i 

contributions arising from configuration mixing. These con¬ 
tributions are of precisely the same type as the ordinary ex¬ 
change contributions. Thus the increase in accuracy that re¬ 
sults from allowing for the dominant effects of configuration 
mixing has not been bought at the price of added complexity, 

There exists a very close analogy between the Hartree- 
Fock theory and the hydrogenic screening theory, which it is 
convenient to exploit because the screening theory is so much 
simpler than the Hartree-Fock theory. The generalized Hartree- 
Fock equations are analogous to the screening equations (C-27). 
Just as in the screening theory one is not justified in neglect 
ing the mixing contributions of v tt ,p beside the term-dependent 
contributions v^p, so in the generalized Hartree-Fock theory 
one is not justified in neglecting the mixing contributions to 

X . , 

nl 

Calculations based on the generalized Hartree-Fock theory 
are planned for the near future. 
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APPENDIX D 

n. T„ : ,.. er 

APPLICABILITY OF THE SCREENING THEORY TO HEAVY ATOMS 
D.l. Introduction 


In the following, the term "screening theory" will be 
used to designate the Z-dependent approximation scheme descri¬ 
bed by Layzer and, in a generalized form that allows 

fully for relativistic effects, by Layzer and Bahcall 
The term "Slater theory" will be used to designate the con¬ 
ventional theory of atomic spectra, whose basic approximations 
were introduced by Slater in his classic paper of 1929f D ’ 3 ^later's^ D *^ 
recent two-volume monograph, "Quantum Theory of Atomic Structure", 
describes the modern form of the theory. 

The screening theory was devised to account for certain 
systematic discrepancies between the Slater theory and experi¬ 
ment, pertaining to the centers of gravity of spectroscopic terms. 
These discrepancies are sometimes large enough to reverse the 
predicted order of term energies in a configuration. 

Now, any reasonably general theory of atomic structure muot 
incorporate some fairly drastic simplifying assumptions The 
Slater theory rests on the assumption that the effects of con¬ 
figuration mixing can be neglected in a first approximation. 

The screening theory, on the other hand, rests on the assumption 
that the energy of an atomic system can be expressed as a rapidly 
convergent power series in (Z - s)~ , where Z is the nuclear 
charge and s is a screening parameter. This assumption leads 
naturally to a theory that is in most respects simpler than the 
Slater theory, but which allows for certain effects of configura¬ 
tion mixing in the first approximation. 
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As yet, the screening theory has been applied only to 
light atoms — specifically to atoms with ten or fewer elec¬ 
trons. In these applications it has proved very successful, 
accounting in a quantitative way for features of term spectra 
that the Slater theory describes only semi-quantitatively or 
qualitatively. The question now arises whether the screening 
theory will also provide a useful description of the spectra 
of heavy atoms. The answer hinges on how rapidly expansions 
of the term energies in inverse powers of Z - s converge. 

The fact that Z - s increases with N (the number of atomic 
electrons) argues for the applicability of the theory to 
heavy atoms. On the other hand, one would also expect the 
coefficients in the energy expansion to increase with N, 

Since it is far from obvious to what extent the second effect 
will offset the first, the question cannot be decided with 
confidence on general grounds? specific calculations are 
required. 

The iron-group elements provide ideal subjects for such 
calculations. They are sufficiently heavy to afford a defini¬ 
tive test. At the same time, they are sufficiently light to 
enable relativistic contributions to the energy to be treated 
as perturbations. Finally, the spectra of many atoms and 
ions in the iron group have been well analysed experimentally. 

In evaluating the merits of a theoretical description, 
one must compare it not only with the data it sets out to 
describe but also with other available theoretical descriptions. 
During the past few years, attempts have been made to patch up 
the Slater theory, as it applies to the spectra of elements in 
the iron group, by adding to the theoretical energy formulae 
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an empirical correction term of the form aL(L + 1). That a 
term of this type could account for a fairly large part of the 
discrepancy between the Slater theory and experiment, in spectra 
of the iron group, was first pointed out by Layzer ^ D " 5 ^ and, 
independently, by Racah. ' The Slater theory, amended in 

this way, has since been used extensively by Trees, Rohrlich, 

and Racah and his pupils in semi-empirical analyses of term 

(r-7) 

spectra in the iron group. Racah has suggested that a 

theoretical basis for the L(L + 1) correction is afforded by 
the work of Bacher and Goudsmit/ 1 ^' 8 * 

In the following, we shall re-examine the theoretical 
basis of the L(L + 1) correction. We shall also present the 
results of new calculations, based on the screening theory, 
for complexes with a partly filled M shell. We shall find 
that the modified Slater theory cannot account quantitatively 
for the effects of configuration mixing described by the screen¬ 
ing theory. We shall also find that the latter are of the same 
order of magnitude as the observed discrepancies. This suggests 
that, the screening theory will provide an adequate description 
of these spectra. Unfortunately, sufficient experimental data 
are not yet available to test the theoretical calculations 
presented here, the complexes here studied having been .hosen 
for their comparative mathematical simplicity rather than for 
comparison with experiment. However, work now in progress 
should lead before long to a definitive experimental test of 
the screening theory as applied to atoms and ions in the iron 
group. 
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Theoretical Results for Complexes with a Partly-Filled 
M Shell * 

We have considered the following complexes: (1 2 2 8 3 2 ), 

.,2.8,2,0 .,2,8,16. .,2,8,16,0 .. _ . , 
(123) , (123 ), (123 ). The configurations belong¬ 

ing to these complexes are listed in Table D-l. These are the 
simplest complexes with partly-filled M shells, other than 
those consisting of closed shells plus or minus a single elec¬ 
tron. Complexes containing from three to seven 3d electrons 
present a somewhat more difficult theoretical problem. 

The evaluation of the energy coefficients and the 
unitary matrices that diagonalize the energy submatrices 
V. , proceeds along the lines laid down in the reference 

" pSL (D-l) 

cited earlier (Layzer ). [w^ is the coefficient of Z 

in an expansion of the nonrelativistic part of the energy in 
inverse powers of Z; V ( n )p SL i® a finite matrix whose elements 
connect states labelled by the indicated quantum numbers 
is an eigenvalue of V ( n )p SL > the radial integrals that figure 
in the matrix elements being evaluated with straight hydro- 
genic wavefunctions.] 

The results of the calculation are present in Tables D-2 thru 

D-9. Tables D-2 thruD5give the coefficients W^ as well as the 

values that would have been obtained for these coefficients 

if configuration interaction had been neglected (i e , the 

values of W^ according to the Slater theory) , Tobies D-6 ‘>ru D-9 

give the unitary matrices that diagonalize the various energy 

submatrices V. , . 

(n)pSL 

* The work summarized in this subsection was supported in part 
by a grant from the National Science Foundation. 
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TABLE D-i 


TEEMS AND CONFIGURATIONS FOR COMPLEXES WITH A 
PARTIALLY FILLED M SHELL 


Complex 

Configuration 

Terms 

uW) 

l* 2 2*2p 6 3« 2 

l s 


3i3d 

V J D 


3p 2 

1 S, 3 P, l D 


3d 2 

l S. 3 P. l D, 3 F. l G 

oW‘) 

l* 2 2e 2 2p 6 3e 2 3p 6 3d 8 

l S. 3 P. l D. 3 F, l G 


3. 2 3d 10 

l 8, 3 P. l D 


3«3p 5 3d 10 

l D, 3 D 


3p l 3d 10 

l s 

oW>° 

le 2 2e 2 2p 8 3e3p 

l P. 3 P 


3p3d 

l P. 3 P, 1 D, 3 D, l F. 3 F 

uWV 

le 2 Ze 2 2p 4 3e 2 3p 5 3d 9 

l P, 3 P, l D, 3 D, l F, 3 F 


3. 3p V° 

V 3 p 







TABLE D-2 


v[ in (1 2 2 8 3 2 ). 

Wj is the 

term-dependent 

part of W^; 

W 1 " W 1 + W 1 ' 

( W l>c.f. denotes 
is neglected. 

the value 

of obtained when configuration mixing 

Configuration 

Term 

W i 

< W i>c.f. 

V<Vcf. 

3s 2 


-.00357 

0 

-.00357 

3s3d 

3 d 

-.00228 

-.00228 

- 


l D 

.01581 

.00684 

.00897 

3p 2 

3 P 

-.00481 

-.00432 

-.00049 


l D 

-.00560 

.00432 

-.00992 


x s 

.01748 

.01727 

.00021 

3d 2 

3 f 

-.00564 

-.00564 

- 



.00295 

.00201 

.00094 


3 P 

.00373 

.00322 

.00051 


l G 

.00616 

.00616 

- 


l s 

.02719 

.02382 

.00337 



TABLE D-3 


v[ in (1 2 2 8 3 2 )° (See Table D-2) 


Config. 

Term 

W 1 

(Wj) c.f. 

W 1 -(W 1 )c.f. 

3r3p 

1 P° 

.01801 

.02116 

-.00315 


V 

-.00753 

-.00705 

-.00048 

3p3d 

1 P° 

.02215 

.01900 

.00315 


V 

.00255 

.00207 

.00048 


V 

-.00869 

-.00869 

- 


3 d° 

.00086 

.00086 

- 


1 F° 

.01934 

.01934 

- 


V 

-.00859 

-.00859 

. 
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TABLE D-4 



W l' 

in (1 2 2 8 3 16 ). 

(See Table D-2) 


Configuration Term 

V 

<V>cf 

w 1 , -(w 1 ‘) 

3p 6 3d 10 

x s 

.00262 

0 . 

+.002620 

3a 3p 6 3d 9 

3 d 

-.00228 

-.00228 

.... 



.01926 

.00684 

+.012424 

3s 2 3 P 4 3d 10 

3 P 

-.00406 

-.00432 

+.000258 



-.00763 

.00432 

-.011951 


x s 

.01625 

.01727 

-.001024 

3s 2 3 P 6 3d 8 

3 f 

-.00564 

-.00564 

.... 


X D 

.00154 

.00201 

-.000472 


3 P 

.00296 

.00322 

-.000258 


X G 

.00616 

.00616 

.... 


x s 

.02226 

.02382 

-.001595 
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TABLE D-5 


in (1 2 2 8 3 16 ). (See Table D-2) 


Configuration Term 

V 

<V>cf 

V’< W l'> 

3s 3 P 5 3d 10 

1 P° 

.02281 

.02116 

.001649 


V 

-.00678 

-.00705 

.000274 

3s 1 3p 5 3d 9 

V> 

.01735 

.01900 

-.001649 


V 

.00180 

.00207 

-.000274 


1 D° 

-.00869 

-.00869 

.... 


3 d° 

.00086 

.00086 

.... 


1 F° 

.01934 

.01934 

.... 


V 

-.00859 

-.00859 

.... 


cf 



TABLE D-6 


THE UNITARY MATRICES | |<$rr|tl>| |(1 2 2 8 3 2 > 


^ 1J 

3s 2 


3p 2 

3d 2 

3. 2 

.99031 


.13789 

-.01670 

3p 2 

-.13490 


.98347 

.12079 

3d 2 

.03308 

Is. 

-.11737 

.99254 


3s 3d 


3p 2 

3d 2 

3 83d 

.39843 


.43093 

-.01285 

3p 2 

-.43732 


.89702 

.06405 

3d 2 

.39642 

2l 

-.05193 

.99786 



3p 2 


3d 2 

3p 2 


.99884 


.04807 

3d 2 


-.04807 


.99884 
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TABLE D-7 


THE UD1TAKY MAItICIS | |<Jrr | tl>| |, (l 1 !*! 1 )® 

Je i 


tssa 

hip 

3p3d 

^vOW 

3«3p 

5p34 

uh 

.Mitt 

.12875 

U*9 

.99877 

.04956 

3p3d 

> .12675 

.99168 

3p3d 

.04956 

.99877 


I 

I 


\ 


SOS 







TABLE 0-8 


THE UNITARY MATRICES | |<trr | tf>| Ul 2 2®3 16 ) 


p5E83B 

3p*3> 

3. J 3p*3d 10 

3. S 3p‘3d 8 

3p 6 3d 10 

.99491 

-.09904 

.01833 


.09991 

.99345 

.05552 

3.V3d* 

-.01271 

.05707 

.99829 


i, 



6 9 

3*3p°3d* 

3. 2 3p A 3d 10 

3. 2 3p‘3d* 

3»3p 6 3d 9 

.831773 

.555116 

.000637 

3* 2 3p 9 3d 10 

-.554800 

.831337 

-.032791 

3. 2 3p*3d 8 

-.018733 

.026921 

.999462 


li^N> ld 

3.V3d 10 

3« 2 3p 6 3d® 

3* 2 3p A 3d 10 

.999707 

-:024220 

3« 2 3p 6 3d 4 

.024220 

.999707 
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TABLE D-9 


THE UNITARY MATRICES | |<fc r r | tf> | ,(1 2 2 8 3 16 )° 


.. 

New^° ld 

3s3p 5 3d 10 


3. 2 3p 5 3d 9 

3s3p 5 3d 10 

.99770 


-.06782 

3s 2 3p 5 3d 9 

_ 

.06782 


.99770 

3P° 

^N° ld 

3s3p 5 3d 10 


3s 2 3p 5 3d 9 

3s3p 5 3d 10 

.99960 


-.02818 

3s 2 3p 5 3d 9 

.02818 


.99960 
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We call attention to the following features of the results. 


(i) In the complexes considered here, the effects of con¬ 
figuration mixing within a given complex are comparable in mag¬ 
nitude with the corresponding effects in complexes of the form 
(1 2 2^) studied previously. Thus, the mixing coefficients, which 
indicate the magnitude of departures from the central field 
approximation, are of the order of .1 and the changes in the 
intervals AW^ between adjacent terms are of the order of 10-30%. 

(ii) In a few instances the changes may be much greater 
than the values mentioned in the last paragraph. For example, 
the configurations 3s3d and 3p 2 in the complex (1 2 2 8 3 2 ) are 
very strongly mixed, the mixing coefficients being about .44. 

The same is true of the corresponding pair of configurations 
in the complex (1 2 2 *j 6 ), where the mixing coefficients are 
about .55. In these cases the resulting change in the value 

of W is great enough to change the original order of the 

1 

spectroscopic term (for large Z) . Thus D) actually lies 
below the value of Wj^P) in the configuration 3p 2 of (1 2 2 8 3 2 ) 
when allowance is made for configuration mixing within the 
complex. The change in for this term thus exceeds 100% 
of the separation AW^ * - W 1 ( 3 P) as predicted by the 

Slater theory. 

(iii) It is necessary in general to allow for all effects 
of configuration mixing within a given complex, though in parti¬ 
cular cases some contributions may be negligible. In particular, 
.interactions connecting configurations that differ in the quantum 
numbers of a single electron cannot in general be neglected be- 
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side interactions that connect configurations that differ through 
the quantum numbers of two electrons. The configurations 3s3d 
2 

and 3p discussed in the last paragraph illustrate this point. 

(iv) The effects of configuration mixing are smaller in 
the complex (1 2 2 8 3 16 ) than in (1 2 2 8 3 2 ). This behavior is similar 
to that observed in complexes of the type (1 2 2**) t The effects 
of configuration mixing tend to diminish with increasing q. The 
reason is that the off-diagonal elements of electrostatic inter¬ 
action are essentially independent of the number of electrons in 
the system, whereas the differences between diagonal elements 
of the energy matrix increase with increasing H (total number of 
electrons). Since the mixing coefficients are, in a first approxi¬ 
mation, equal to the ratios between off-diagonal interaction 
elements and diagonal energy differences, they will tend to diminish 
with increasing N. 

p.l. . .Pm Saal-Bapirigal Eftfrod of Bacher and Gpgdamjt 

The method of Bacher and Goudsmit * D " 8 * seeks to relate 
the term energies of an N-electron atom with the term energies 
of its ions, eo that if the appropriate term energies of systems 
(M,Z) with N < N are known (from either theory or experiment) 
one can calculate the term energies of the system (N,Z) to a 
certain approximation. 

The theoretical relations connecting the N-electron ener¬ 
gies with the 1-, 2-, ... , (N - 1)-electron energies are based 
on the independent-orbital approximation. That is, each electron 
is assigned a definite pair of orbital quantum nusibera nl. We 
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now derive these relations, using a more explicit notation 
than that of Bacher and Goudsmit. 

Let denote the number of Is orbitals, the number of 
2s orbitals, q^ the number of 2p orbitals, and so on up to q^, 
so that for an N-electron atom 

q x ♦ q 2 ♦ ••• + q r - * • 

Let r N denote the set of quantum numbers that are required, in 
addition to the q^, to specify the state of an N-electron atom. 

Consider a given state (q^)^* An M-electron subsystem, 

consisting of the nucleus and M < N electrons, will not in 

1 M i 

general be in a definite state (a ... a The a denote 

pairs of orbital quantum numbers, not necessarily all distinct, 
selected from the set of N pairs that characterise the N- 

electron state under consideration. We shall use the abbrevi- 
1 M 

ation (a ... a ) » a M « Thus we may use the symbol in 

place of (q x ... q r ) r „- With every state we may associ¬ 

ate a definite probability p I a M r ' M l a I | r N l • More generally, we 
consider the set of all probabilities of the type p ( a A r A l° B r B J» 
where A < B < N and the set a A is contained in the set a B . 

These probabilities are not all independent, since they evi¬ 
dently satisfy all relations of the form 

Z P(r A |r B ),tr B |r c ) " P(r A |r c> ' <»- l l 

B 

(For simplicity we have omitted the labels a.) 

By means of what are now standard methods in the theory of 
complex spectra, one can evaluate any probability p (° A r A l® B r B )• 
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£ Pfa M r M |a N r N )W<a M r M |a NV 

r 


is the mean energy associated with the subconfiguration (* M 
when the N-electron system is in state r^. Note that the ener¬ 
gies w ( a M r M l a N r N ) do not coincide with the energies of the cor¬ 
responding M-electron ions of the N-electron system which would 
be denoted by w » because the electrons in the sub¬ 
system move in the combined field of a nuclear charge Z and N-l 
orbital electrons, while the electrons in the ion move in the 
combined field of the same nuclear charge and M-l orbital elec¬ 
trons. Thus the energy levels of the subsystem lie above the 
corresponding levels of the ion. However, in practical appli¬ 
cations of the method to be described, Bacher and Goudsmit used 
empirical ionic energies as approximations to the corresponding 
subsystem energies. The latter cannot, of course, be obtained 
directly from experiment. We shall return to this question 
later. 


Following Bacher and Goudsmit, we write the energy of an 
M-electron subsystem as a weighted sum of contributions w ( a K r K )* 

K < M, where cl stands for a set of K labels a 1 selected from 
K 1 M 

the set of M labels ci M = a ... a . The energy contributions 
w are defined implicitly by the set of equations 

w< w - I w(a) ♦ ♦ ••• + 
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(»- 2 ) 



The mean energy contributions W are defined in analogy with 
the Mean energies Wt 

“ (0 K la M r M> ‘ Z Pt “K r K l “M r M )w<a K r K ) 
r K 

The set of all these equations, for all states with N < II, 

has a unique solution for the energy contributions w in terms 
of the energies Ws 


“W - "'Vh 1 - I W(a M-llVi + 

°M-1 °M-2 

- ... - <-l) M £w<a> 
a 


In this equation all energies pertain to subsystems of the 
M-electron system in the state The proof of the formula 

is toy Mthematical induction. For M - 1, the formula reduces to 


W(a) « w(a) 


which coincides with equation (D-2) for M ■ 1. Suppose now 
that the formula holds for M » 1, ... , n. We now show that 
it must then hold for M - n+1. From equation (0-4) we have. 


- I'< a „ r „l“n + I r n + 1> * 

“„ r n 

x <“«% r n> ' I P ' a n-i r „-ll a „ r n»"«*n-i r n-X> + - 1 


(D-3) 


(0-4) 


(0-5) 
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L Pla n-i r n-ll a n+i r n+l } 


X f"«Vi r n-l> - Z I> ta n - 2 r n- 2 l“„ r n ,W( “n-2 r „-2» 


Using the probability-composition rule (0-l}-« we can put the 
right side of this relation into the form 


“‘“n+lW + 


n+l-k^n+l^+l* * 


The number of individual contributions to the kth summand is 
C k ( n JJ 1 ) . This number is also given by 

k k+1 

j-0 j-0 

Using the addition theorem for binomial coefficients, we obtain 


I <- i,3 < n £ 3 >C- 3 > * <- i > k+i O - 0 • 

j-0 


Sc - <“*> 

which completes the proof. 


Bacher and Ooudsmit assume that, for a given state o^r^, 
the energy contributions 

I-KW 


]... . 

(D-6) 

(D-7) 


( D-8) 


(D-9) 

( “ 10 ) 
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diminish rapidly with k. Consider the expansion 


W( “N r N> ' * F W<a 2 |a H r H> + ••• + F ' W(0 kl“N r N> + 

“ “2 “k 

+ ... H- w(a H r M ) , 

which is a special case of equation (D-2). By omitting the 
last term on the right we obtain an approximate formula for 
the energy of the N-electron system in terms of quantities 
pertaining to k-electron subsystems with k ^ N. From equa¬ 
tion (D-4) it follows that the approximation w ( a N r 1 |) * 0 is 
equivalent to the approximation 

W( Vh> - - I“ (a N- 2 IW + ---tZ” la ' a N r N ) • 

"n- 1 a N-2 “ 

Thus the energy of the N-electron system is expressed as a 
linear combination of energies of k-electron subsystems with 
k < N. 


More generally, we obtain the kth approximation in the 
sense of Bacher and Goudsmit. when we retain the first k+1 
terms on the right side of equation ('0-11). I-i particular, 
the first approximation (which is exact for 2-electron sys¬ 
tems) is defined by the equation 

w(1>(a N r N > " Z W(a> + * 


£w(a 2 |a N r N ) - (M-2) £*(<*1*^) 


(/0-U) 


(0-12) 


(0-13) 
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Me) 

With the help of equation (D-4) we can express W , for 
any k < N, as a linear combination of energies pertaining 
to (k+1)-electron subsystems. 

In practical applications of the Bacher-Goudsmit approxi¬ 
mation method one uses experimental values of the ionic ener¬ 
gies w ( a ] c r ] c l ct j c r ] c ) in Place of the subsystem energies w ^ ct j c r j c l a N r N ) • 
In order for this to be a good approximation, the electrons 
that are present in the N-electron system but not in the k- 
electron subsystem must make only a minor contribution to the 
electrostatic field seen by the remaining electrons. Since 
the removal of "core electrons" (electrons in closed shells) 
does substantially alter the effective electrostatic field of 
the atom, one should always regard the core electrons and the 
field they produce as fixed. That is, one should reckon all 
energies from the core energy as a fixed zero-point. The 
formulae already developed apply without change to these 
relative energies, since they do not depend on any assump¬ 
tions regarding the nuclear contribution to the electrostatic 
field except the assumption of radial symmetry. Since the 
core, by definition, is composed of complete groups of elec¬ 
trons, its contribution to the field will have spherical 
symmetry and may therefore be lumped with the nuclear con¬ 
tribution . 

_Ilie-Bacher-Goudsmit Mfthgd apd ^gnliguration^Lxina 

The independeht-orbital approximation is just as basic 
in the theory of Bacher and Goudsmit as in the theory of 
Slater. The basic equations of the theory, (D-2) and (D-4), 
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presuppose the possibility of assigning every subsystem of an 
N-electron system to a definite electronic configuration, and 
the coefficients in these equations are to be evaluated by 
means of a theory based on the same approximation. Neverthe¬ 
less, the Bacher-Goudsmit approximation does, in a sense, 
allow for certain effects of configuration mixing, because 
it makes use of empirical ionic energies, since these include 
contributions from configuration mixing, so must the "predicted" 
energies. Now, the approximation of greatest practical and 
theoretical interest in the Bacher-Goudsmit scheme is the first 
one, which yields the energies W^ defined by equation (V-13). 
The question therefore arises: Which effects of configuration 
mixing are included in when the two-electron energies 

w^^V are taken from experiment? 

Bacher and Goudsmit discuss this question briefly and 
schematically, using ordinary perturbation theory. They con¬ 
clude (incorrectly, as we shall see) that the first approxi¬ 
mation includes all effects of configuration mixing with 
doubly excited configurations (i.e., configurations differ¬ 
ing from the original one in two pairs of one-electron quan¬ 
tum numbers). They also conclude (correctly) that the first 
approximation does not include all effects of configuration 
mixing with singly excited configurations. Before turning 
to a detailed consideration of these questions, however, let 
us examine the problem from another angle. 

As has previously been pointed out (Layser the 

Slater theory predicts certain qualitative regularities in 
the structure of term spectra, which do not accord with the 
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observed regularities. For example, the Slater theory pre¬ 
dicts that all configurations consisting of n or 6-n equiv¬ 
alent p electrons, together with any number of complete 
groups, have the same structure apart from a scale factor 
Specifically, if n = 2, 3, or 4, the two independent inter¬ 
vals between the spectroscopic terms with L * 0, 1, 2 are 
predicted to satisfy the relation SD/DP = 3/2. The Z- 
dependent approximation scheme predicts different quali¬ 
tative and quantitative regularities. It predicts, for 
example, that the interval ratio SD/DP approaches a con¬ 
stant value with increasing Z in a given isoelectronic 
sequence, and that this constant (which can be calculated 
from the theory) depends both on the principal quantum 
numbers of the p electrons and on the complete groups pres¬ 
ent in the electronic configuration How does the Bacher- 
Goudsmit theory stand on this matter? 

2 2 2 

Consider a simple example, the configuration Is 2s 2p 
in the carbon sequence. We regard the K-electrons as com¬ 
posing a fixed core, and proceed to calculate the terra 
energies relative to this core in the approximation (D —13) 

We obtain 

(1) 2 2 2 — 2 
W VA '(2s*2p SL) * W(2s ) + 4W(2s2p) + W(2p SL) - 4W(2s) - 4W(2p), (D-14 

where we have omitted unnecessary labels on the fight side 
Note that the mean two-electron energy W(2s2p) does not depend 
on the quantum numbers SL, because the two 2s electrons form 

a complete group. From this equation we see at once that the 
2 2 

configuration 2s 2p is predicted to have the same structure 
as the configuration 2p*. Actually, the interval ratio SD/DP 
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is substantially less than 3/2 in all systems belonging to 

2 2 

the beryllium sequence (i.e., in all configurations Is 2p ). 
The difference between the two limiting ratios is due to 
the mixing of configurations belonging to the same complex; 
both limiting ratios are correctly predicted by the Z-depen- 
dent theory. We conclude, therefore, that the Bacher-Goudsmit 
method does not in general allow correctly for this type of 
configuration mixing — the most important type according to 
the Z-dependent theory — in the first approximation. 

Now, in the example just given the "perturbing" con¬ 
figurations that need to be considered if .the limiting inter¬ 
val ratios are to be correctly predicted are 2p (interacting 

2 2 2 2 
with 2s 2p ) and 2s (interacting with 2p ). These are 

"doubly excited" configurations in the sense explained above 

(though the configuration 2s 2 should perhaps be described 

as doubly de-excited relative to 2p 2 ). It follows that the 

first approximation in the Bacher-Goudsmit scheme does not, 

in fact, allow fully for the effects of mixing with doubly 

excited configurations. 

The argument given by Bacher and Goudsmit fails because 
it does not allow properly for the Pauli exclusion principle. 
Bacher and Goudsmit tacitly assume a strict correspondence 
between doubly excited N-electron configurations and doubly 
excited configurations in which the principal quantum num¬ 
bers of the excited electrons differ from those of the orig¬ 
inal electrons. It does not, however, extend to doubly 
excited configurations with the same principal quantum num¬ 
bers as the original configuration. Thus, the de-excitation 
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of the 2p electrons in the configuration 2p 2 to form the 
2 

configuration 2s is permitted by the exclusion principle; 

2 2 

but in the configuration 2s 2p , de-excitation of the 2p 

4 

electrons would yield the non-permitted configuration 2s . 

The predictions of the first-order Bacher-Goudsmit 
theory differ from those of the first-order Z-dependent 
theory in other important respects. For example, we have 
seen (see subsection 2) that the effects of mixing with 
singly excited configurations may be very large. We have 
also seen that in some instances the effects of configura¬ 
tion mixing are too large to be adequately described by 
second-order perturbation theory Finally, we have seen 
that the effects of configuration mixing are systematically 
greater in almost-empty shells than in the corresponding 
almost-fullshells. All these features of the Z-dependent 
theory are at variance with the first-order Bacher-Goudsmit 
theory. 

We conclude that the Bacher-Goudsmit theory does not 
provide a satisfactory framework for dealing with the most 
important effects of configuration mixing, even semi-em-. 
pirically. 

P .5. The Modified Slater Theory 

The theory of Bacher and Goudamit is held by Racah * D ‘ 7) 
to provide a theoretical justification for the 
L(L + 1) correction in the following sense. By means of 
equation (D-13) one can express the energies of an N-elec- 
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tron system in terms of the energies of two-electron sys¬ 
tems. Let us admit, for the sake of the argument, that the 
dominant effects of configuration mixing for the N-electron 
system are included in the energies calculated with equation 
(D-13) if one uses the true values of the energies 
in evaluating the right side Instead of using experimental 
values for the last-named energies, one could use an interpo¬ 
lation formula containing a number of adjustable parameters 
Suppose we are considering configurations of the type d n 
The configuration d 2 has five terms, so that if one uses an 
interpolation formula with five adjustable parameters, the 
values of the parameters will be uniquely determined by 
experiment. Equation (D-13) then enables one to express 
the energy of any term belonging to the configuration d n 
in terms of these five parameters * According to Slater's 
theory,, the mean energy of all the terms in a given con¬ 
figuration labelled by given values of S and L is a linear 
combination of certain radial integrals, the Slater P and G 
integrals. For configurations d n there are three such 
integrals; F°, F 2 and F 4 . However, even if one regards 
these as adjustable parameters, one cannot obtain a very 
good representation of the experimental data. The addition 
of two empirical correction terms, aL(L ♦ 1) * 0Q, the second 


*In practice the parameters are not regarded as having 
fixed values, corresponding to the configuration d 2 , but 
as functions of the effective central field. One can 
evaluate the parameters for n > 2 by equating the pre¬ 
dicted energies to their experimental counterparts and 
solving the resulting equations by the method of least- 
squares . 
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of which was suggested by Racah^" 6 ^, makes it possible to represent the 
2 

term energies in the configuration d exactly and to represent the term 
energies in other configurations that have been examined with fairly good 
accuracy. 

From the discussion in the preceding subsections, however, we know 
that equation (D-13) does not provide a satisfactory framework for dealing 
with the dominant effects of configuration mixing. It is true that the 
empirical description of spectra in the iron group based on the Slater- 
Racah interpolation formulae has been remarkably successful. However, 
guided by the considerations of the present paper, one can easily choose 
examples of configurations for which the method will give poor results. 

A simple example is the configuration ls 2 2s 2 2p 6 3p 2 . 
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APPENDIX E 
D. Layzer 

RELATIVISTIC EFFECTS IN HEAVY ATOMS 

E. 1 Introduction 

In heavy, as in light, atoms the chief groupings of energy levels 
have their source in the nonrelativistic (electrostatic) interactions 
between the electrons. These groupings are characterized by sets of 
prinicpri quantum numbers (•»*). Adjacent groups are separated by 
intervals of the order of Z^ cmhc^, where Z g is the effective nuclear 
charge for the outermost electrons. 

The structure of a given group, however, has an entirely different 
character in heavy atoms rnd in light atoms. In the latter, each group 
(or complex ) splits up into spectroscopic terms aSL consisting of levels, 

dlsti iguished by the quantum number J. The spacing between terms is of 
2 2 

the order of Z g a me , while the spacing between levels in a given term 

4 4 2 2 3 

is of the order of Z g a me . In light atoms, by definition, a « 1, 

so that the fine structure is clearly differentiated from the coarse 

structure. The fine structure originates in relativistic interactions 

between the orbital electrons and the nucleus, t'> a course structure in 

nonrelativistic interactions between the electrons. 

2 3 

In heavy atoms, on the other hand, the coupling parameter a Z # is 
not small. Consequently the relativistic interactions cannot adequately 
be treated by first-order perturbation theory, as they can be in light 
atoms. The structure for the major groupings (complexes) is correspondingly 
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complicated, in general, there is no clear-cut distinction between 
coarse structure and fine structure; the two kinds of interactions con¬ 
tribute comparably to the spacing between levels in a given complex. 

3 2 

When Z # a » 1, the relativistic contributions become dominant; 
the complex breaks up into sub complexes , characterized by quantum 
numbers (n*jS- Within a given sub complex the spacing of levels is 
due primarily to the nonrelativistlc (electrostatic) interaction. 

This state of affairs is characteristic of x-ray spectra. Here Z g 
may be comparable to Z itself because the missing electrons belong to 
inner shells. 

The conventional method for describing the structure of energy 

2 3 

spectra in the intermediate-coupling regime (a z @ 1) is set forth in 
(E-l) 

Theory of Atomic Spectra. It is based on the central-field ap¬ 

proximation and allows for the leading relativistic contribution to 
the energy - the spin-orbit interaction. This is evaluated for electrons 
moving, not in the Coulomb field due to the nuclear charge Z, but in an 
effective central field due to the nuclear charge plus the electrons. 

The conventional method has several drawbacks. 

It is better suited to a semi-empirical description than a 
purely theoretical one. This is because the radial integrals (|) 
that determine the magnitude of the spin-orbit interaction cannot be 
accurately evaluated by purely theoretical means. They ere best 
evaluated by treating them as adjustable parameters, to be assigned 
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value* that will make the theoretical formulae fit the data as well 
as possible. This Is nearly always the adopted In practice. 

The spln-orblt Interaction appears as a separate term In the 
Haailltonlan for a many-electron atom only In the Pauli (low Z) ap¬ 
proximation. Whan Z Is not small, It is by no means clear that the 
spin-orbit Interaction adequately approximates to the relativistic 
part of the Hamiltonian. 

An alternative description of relativistic effects In the 

structure and spectra of many-electron atoms is provided by recent work 
(E-2) 

of Layzer and Bahcall This Is an extension of the nonrelatlvlstlc 

Z-dependent theory^ 8 " 3 ^. The principal features of the Z-dependent 
relativistic theory are as follows: 

(a) The Hamiltonian eaployed is shown to contain all known 

relativistic interactions up to (but not including) those of order 
4 2 2 

a Z me . The largest known Interactions not Included In the Hamiltonian 
are those responsible for the lamb shift. 

(b) The spln-orblt Interaction Is allowed for exactly, In 
the seroth approximation - Just as in the Dirac theory of the hydrogen 
atom. 

(c) The theory enables one to evaluate all contributions of 
order Z** compared with the spln-orblt interaction. This Is not pos¬ 
sible with the conventional theory - and herein lies the crucial dif¬ 
ference between the old and the new methods. It is the possibility 

of precisely evaluating the next-to-the-leading relativistic contributions 
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that renders the new method quantitative where the old one was merely 
seal-quantitative. 

Although the next-to-the leading relativistic contributions can be 

evaluated in principle, the actual calculations are long and difficult. 

2 

Ivan for such coaparatively slaple electronic configurations as la 2s 
and ls 2 3s, considered by Layser and Bahcall^ the evaluation of those 
terns required a long calculation. For heavy atoas the corresponding 
calculations are at present prohibitively difficult: It is therefore of 
Interest to inquire whether there is an approxiaatlon that will enable 
one to estlsute the contributions in question with fair accuracy. We 
deal with this question in Section E.2. In Section E.2, we consider a 
related problem (as it turns out to be)-the relativistic generalisation 
of the nonrelativlatlc screening paraaeters defined by Layser* 8 ’ 3 *'. 




2 4 2 

The coefficient of a Z ac in the expression for the energy of an 

atoaic atata, according to the theory of Layser and Bahcaii**’ 2 *, 

2 3 

is a function of the coupling paraaater X ■ a Z . In the extreme rela¬ 
tivistic case, this coefficient, which we denote by a^, is given by a 
direct generellsatlon of Soaaarfeld'a formula: 



In the opposite Halt (X - 0) a^ does not approach the value that would 
be predicted on the basis of tha central-field approxiaatlon} to obtain 
ths correct limiting valuq one swat allow for tha effects of configuration 


m 


mixing. (This can be done most simply by using, in place of the central* 
field wave function \|r, a finite linear combination of wave functions 
£*1^, the constants a^ being given by the screening theory; see 
Appendix C. 

In the general case a^(X) can be calculated, without excessive 
labor, to any desired accuracy by methods described in the paper by 
Layser and Bahcall cited earlier. Me may, therefore, consider it to 
be known. 

4 3 2 

It is the coefficient a^(X) of the term a Z me that presents the 
real problem, and which we would like to approximate. 

Now in the standard theory of one-electron spectra (including x-ray 
spectra), the relativistic contribution to the energy is usually 
written in the form 

a 4 (Z - o) 4 amc 2 + 0(a 6 Z 6 ) (E-2) 

where a^ is given by Sommerfeld * s formula, 

*':b(r**-*) <E - l,) 

and o is a suitably chosen constant. Evidently 



Thus a knowledge of o is equivalent to a knowledge of a^. 

Now some, though not all, of the contributions to a^ can be inter¬ 
preted physically in terms of screening. For exasple, the spin-other- 
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orbit interactions, which contribute to a^, tend to reduce the magnitude 
of the (true) spin-orbit interaction in much the same way as the inter- 
electronic electrostatic interactions reduce the magnitude of the electro¬ 
static inter .ctions between the nucleus and the electrons. On the other 
hand, such two-electron interactions as the spin-spin interaction have 
no one-electron counterpart. The interpretation of a as a screening 
parameter is therefore an oversimplification. Nevertheless, it may 
furnish a useful means of estimating the value of a^ through Equation (J>3). 

The simplest assumption one can make about o is to equate it to S, 
the screening parameter that figures in the screening theory. A com¬ 
parison of emprical values of o for hydrogenic spectra with theoretical 
values of S shows that this a a fair approximation, though not as good 
as one might wish. 

Fortunately, the fractional errors due to inaccuracies in the 
estimated value of o tend to become progressively smaller with increasing 
Z because the quartic term becomes progressively larger in comparison 
with the cubic term. Thus, the screening approximation to a^ should 
be considerably better for Uranium, say than for Lithium. 

A semi-empirical study to improve the estimate of a might proceed 
along the following lines: 

The two-electron screening parameter S(a/p) [representing the 
screening of an a electron by a p electron] is given in a first ap¬ 
proximation by the formula^ 6 
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S<Of/p) 


I dr j Pp (r 1 ) dr 1 


(E-4) 


This can be written as 


S(a/p) «j w(r 2 ) P < rjj dr 2 (E-5) 

i.e. as a certain weighted average of the probability P that 

a 0 electron will be found closer to the nucleus than r , the distance 

a 

of an a electron. The weight w is given by 

“< r n> • »a< r a> r a ' < E ' 6 > 

This weighting factor is appropriate for electrostatic screening, but 
not necessarily for other kinds of screening. A more general law than 
(E-6) la: 

“k (r a> ' p a < r a> 'a'" < E - 7) 

It would be interesting to test whether, by an appropriate choice of k, 
one could arrive at a formula for a in better agreement with experiment. 
[The connection between the total screening constants S Q and the two- 
electron screening parameters is discussed in the paper by Layzer^® ^ 
cited earlier. ) In the meantime, the old definition (S - o) can be 
used with a reasonable expectation of fair accuracy. 
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E.3 Relativistic Screening Parameters 


The screening parameters that we have dealt with so far are defined 
by a nonrelativistic theory. The question naturally arises as to what 
constitutes the proper relativistic generalization of the screening 
parameters. 

To begin with, the relativistic screening parameters must depend 
on the coupling parameter X as well as on Z. They must, in fact, have 
the form 

S(X,Z) - S q (X) + S_ 1 (X)Z" 1 + 0(Z" 2 ). (E-8) 

In the limit X -* 0, S(X,Z) -* the nonrelavistic screening parameter S(0,Z). 
The coefficients S^(X) are all bounded from above and below. In the 
relativistic limit X -* oo, 

S(X,Z) - S q (oo) (E-9) 

since the remaining terms on the right side of (a) vanish in this 
limit. 

It is possible to develop a detailed theory of the X-dependent 
screening parameters along the lines of the variational nonrelativistic 
theory. The main result is that the contributions to the screening 
parameters from electrons in closed shells are independent of X. 

Now the X-dependence of the screening parameters could be significant 
only in heavy atoms. Since most of the electrons in a heavy atom are 
nonsally found in closed shells, the amplitude of variation of the 
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screening parameters with X must be small, so that the nonrelativlstlc 
approximation to the screening parameters should be adequate for most 
practical purposes. 
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APPENDIX F 

D. Layzer and A. Naqvi 

STATISTICAL DESCRIPTION OF COMPLEX SPECTRA 

F. 1. Mean Energies of Configurations and Complexes 

In any partial description of complex energy spectra the 
mean energies of configurations and complexes play a central 
role. Even if no other information about the energy spectrum 
is available one can still obtain from the mean energies, with 
the help of Saha's equation, useful estimates for the concen¬ 
trations of the various atomic and ionic species in a mixture 
of gases. Again, in constructing radial wavefunctions for 
use in the evaluation of absolute transition .probabilities — 
either by the screening approximation or by a more refined 
method — one may not wish to use a different set of wave- 
functions for each distinct level. In a first approximation 
one could then use a single set of wavefunctions correspond¬ 
ing to the mean energy. 

If we neglect entirely the effects of configuration mix¬ 
ing, we can write the mean interaction energy of any configu¬ 
ration as a sum of two-electron contributions (nl,n'l'). If 
the atom has N electrons the mean interaction energy consists 

of (j) » N(N - l)/2 3uch contributions. For example, the 

* 2 2 5 

mean energy of the configuration Is 2s 2p is 

(Isis) + (2s2s) + 10(2p2p) + 4(ls2s) + 10(ls2p) + 10(2s2p) . (F-l) 

Here N ■ 9 and the number of two-electron interactions is 
9 X 8/2 » 36. 

The two-electron interactions can be expressed in terms of 
the Slater F and G integrals. The appropriate formulae are 
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given by Layzer * F " 1 ^ and Slater ^ F ' 2 \ but are implicit 
in, earlier work by Condon and Shortley and Short ley. ^ F "*^ 


To calculate the total mean energy of a configuration, 
including the one-electron part, we may conveniently employ 
the Z-dependent screening theory. The mean energy is then 
given by the formula 

5 - - I *nl <I n/ JnJ ' 

Here is the number of electrons of the type nl in the 
configuration under consideration and 

Z , a Z - 8 , . 

nl nl 

The screening parameters s^ are defined in terms of the 
two-electron interaction energies: 

•nl ■ I'Vv - 6 »l,nT ,,(nl|, ' l '> ' 
nr 


n 2 d{nl,n' l' )/d* nl (nljfn'l') 

s(nl|n’l') * 

£n 2 d(nl,nlJ/dz^ (nl -n'l') 

The two-electron screening parameters depend on Z, approach¬ 
ing constancy as z •* m. 


(F-2) 


(P-3) 


(P-4) 


(P-5) 
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The preceding discussion can easily be extended to com¬ 
plexes. The advantage of doing this is that the resulting 
formula for the mean energy of a complex will have the cor¬ 
rect linear dependence on Z. On the other hand, if one is 
dealing with a complex that contains a large number of con¬ 
figurations, the extra work required to allow properly for 
ell the interactions may offset the resultant gain in accuracy. 

F. 2. The Coefficients for Complex Configurations 

Instead of working with the energies themselves, one may 
wish to consider the statistical properties of the energy co¬ 
efficients W , defined as the coefficient of Z in an expansion 
of the energy in powers of Z . In this subsection we discuss 
briefly some relevant properties of the W^, before returning, 
in the following subsection, to the main stream of our discus¬ 
sion. The considerations of the subsequent subsections will 
apply equally to the energies themselves and to the coeffici¬ 
ents W 

The coefficient is composed of two parts. The two- 
electron part can be calculated from the two-electron 
interactions and is the same for all levels belonging to a 
given configuration. The remainder of W^, which we denote 
by W', results partly from LS-dependent interactions within 
the configuration and partly from interactions connecting 

terms belonging to different configurations in the same com- 

(F*l) 

plex. It has been shown (Layzer ) that in atoms with a 


233 



partly-filled L shell these two contributions to are often 
comparable in magnitude Some recent, as yet unpublished, cal¬ 
culations show that the same is true of atoms with partly-filled 
M shells Hence the effects of configuration mixing within a 
complex cannot in general be neglected beside the LS-dependent 
contributions from interactions within a given configuration. 

As a preliminary to calculations of thai allow properly 
for the effects of configuration mixing, and as a temporary sub¬ 
stitute for such calculations we may conveniently evaluate the 
contributions to from interactions within the configuration. 
This is by no means a trivial matter for complex configurations, 
since such configurations may contain several distinct spectro¬ 
scopic terms characterized by given values of S and L. In the 
representations normally used the energy operator is diagonal 
with respect to S and L but not with respect to the quantum 
numbers that serve to distinguish spectroscopic terms with the 
same S and L Hence the energy matrix, instead of being dia¬ 
gonal. will in general consist of multi-dimunsional blocks 
This is the form in which energy calculations have usually been 
left in the past. In order to complete the diagonalization 
one needs to know the ratios of the .various Slater integrals 
that figure in the energy matrix elements In the standard 
theory of complex spectra these ratios are regarded either as 
parameters to be determined semi-empirically, or as theoretical 
quantities whose evaluation requires the solution of a set of 
Hartree-Fock equations On the other hand m the Z-dependent 
scheme one evaluates the Slater integrals using straight hydro- 
genic wavefunctions The contributions to from interactions 
within a given configuration have definite numerical values 
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which can be calculated once and for all, as can the elements 
of the transformation matrices that diagonalize the blocks re¬ 
ferring to spectroscopic terms of given S and L. The trans¬ 
formation to a representation in which has no nondiagonal 
elements connecting terms in the same configuration is the 
first step toward the ultimate objective of a representation 
in which the energy matrix has no nondiagonal elements connec¬ 
ting terms belonging to the same complex. 

P. 3. The Variance of the Energy 

The variance of the energy (in a given configuration) is 
defined by 

varW =<(W - W) 2 > = o 2 . (F-6) 

Moszkowski * F ‘ 5 ^ has written down formulae that give o w for 
an arbitrarily complex configuration in terms of the o^s for 
two-electron configurations. The latter can easily be evalu¬ 
ated by direct calculation. Proofs of Moszkowski‘s formulae 
have not as yet appeared; but for the purposes of this report 
we shall assume that the formulae are correct — or at least 
that can be found for any configuration of interest. 

The parameter o w measures, roughly, the width of a con¬ 
figuration. Knowing o w , one can estimate (from Boltzmann's 
formula) the variation in the population of levels over a 
given configuration. More generally, one can estimate the 
error that results, in any particular calculation, from 
assigning the same energy to all levels in a given configu¬ 
ration. From one can also derive a rough estimate for 
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the average spacing between levels in a configuration. Thus 
if one assumes that, say, 50% of the levels lie in the range 
W ± o w , one can reasonably expect the resulting estimate of 
the level spacing to be accurate to within a factor of 2. 

F. 4. The Distribution Function of the Levels in a Configuration 

Knowing o we know the width of the probability distri¬ 
bution, but not its shape. The question now arises: What 
is the best assumption as to the shape of the distribution 
that is consistent with our limited knowledge of the energy 
spectrum? 

In addition to W and o we know from direct calculations 
w 

that the distribution of energy levels in a configuration is 
strongly skewed. This fact has suggested (unpublished Convair 
report, 1961) that the distribution can be represented by an 
exponential function of the energy. Once the exponential form 
has been chosen, the two parameters on which such a distribu¬ 
tion depends are uniquely determined by the values of W and o^. 

In fact, the distribution is given by 

N exp [-(w - W)/o w ] (W > W - o w ) 

P(W) = (F-7) 

0 (W < W - o w ) 

Thus the energy spectrum is supposed to have a sharp cut-off 
a w units below the mean energy W. The probability distribution 
achieves its maximum value at the cut-off and then falls off 
exponentially. 


236 



Before this law, or any such law, can be compared with experiment 
we need to have some reasonable and consistent way of passing from the 
theoretical or experimental energy spectrum, which is discrete, to a con¬ 
tinuous probability distribution. One such way is the following. We 
assign to each level the appropriate statistical weight (2S + 1). Then 
we divide the energy range into a convenient number of equal intervals, 

bounded on the left by the lowest energy in the configuration, and 
determine the total statistical weight of all levels falling in each 
energy interval. We then plot a histogram of the running means of the 
statistical weights against the energy. This corresponds to smearing 
out the discrete energy levels, which occur in bundles, into a somewhat 
continuous distribution of energy. As long as the intervals used to 

smooth the data are reasonable small compared with the total energy 
range, the exact manner of smoothing is of no importance. One should 
bear in mind that, as discussed in Subsection F.2, the distribution of 
energy levels as calculated with neglect of configuration mixing is at 
best semi-quantitative, so that much of the detail present in an un- 
smoothed plot of the distribution of energy levels is spurious anyway. 

Figures F-l thru F-6 Illustrate the smoothed distributions of 

energy levels in configurations consisting of equivalent d and f elec- 

(F-8) 

trons. Experimental data from the compilation by Hiss Moore , are 
taken for the d-electron configurations and theoretical energies are 

taken for the f-electron configuration from the paper by Elliot, Judd 
(F-9) 

and Runciman. If in the experimental data all theoretical 

terms are not identified, the missing terms are given in each 
figure. These plots exemplify the skewed character of energy dis¬ 
tribution, but the skewness appears to be much less marked 
than that associated with the exponential 
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Figure F-l. Smoothed Distribution of Experimental Energy 
Levels in V II (Configuration 3d^) 
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Figure F-2. Smoothed Distribution of Experimental Energy 
Levels in Cr II (Configuration 3d^) 
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Figure F-3. Smoothed Diecribucion of Experimental Energy 
Levels in Nb II (Configuration 4d^). 
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distribution. In every case the main peak of the distribution 
lies significantly to the right of the lower cut-off. 

Another defect of the exponential distribution consists in 
its not having the lower cut-off in the right place. The lower 
cut-off has in fact been arbitrarily placed units below the 
mean energy W. This assignment is usually in error by a quan¬ 
tity not much smaller than o itself. Thus the main features 
J w 

of the distribution are not as well represented by the expon¬ 
ential law as one would like. 

We have sought to improve the statistical description in 
two ways. First, we have added a third parameter, which can 
in every case be calculated exactly with comparatively little 
labor. Secondly, we have chosen the distribution function in 
accordance with an objective principle — the principle of 
maximum entropy inference — which ensures that the probability 
distribution incorporates the maximum possible uncertainty 
consistent with having the correct mean, variance, and lower 
cut-off. These two improvements are discussed in detail in the 
following subsections. 

F. 5. The Lower Cut-off 

The lower cut-off is simply the energy of the ground 
state of the configuration. Now, according to a well-known 
rule of Hund, the lowest spectroscopic term in a given con¬ 
figuration is the one with the highest value of S; or, if two 
or more terms have the same value of S, the one of these with 
the highest value of L, One of the levels belonging to this 
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term, namely, the level with M g = S and = L, is charac¬ 
terized not only by definite values of M_, S, M_, and the one- 

i L S 

electron azimuthal quantum numbers l , but also by definite 
values of the one-electron magnetic quantum numbers m*. Hence 
the corresponding wavefunction is a Slater determinant, and 
the corresponding energy can be obtained by a straightforward 
calculation. As an example we calculate the energy of the 
ground state of the configuration f 7 . (See Section F-9) 


F. 6. Maximum Entropy Inference 

Consider the problem of determining the probability dis¬ 
tribution p(x) associated with a given random variable X. 

Suppose that the possible values of X are know to lie in the 
range (a,b) and that the first k moments of X are also knownt 

b 

Hi - f x i p(x)dx (i - 1, ... ,k) . (F-8) 

a 

p(x) must satisfy the further conditions 
b 

P(x) >0, J p(x)dx - 1 . (F-9) 

a 

According to the principle of maximum entropy inference (Jaynes 
the best choice for the function p(x) is the one that 
maximises the entropy 

b 

H - - J p(x) in p(x)dx , (P-10) 
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subject to the k + 1 auxiliary conditions expressed by the 
above integral equations. Using the method of Lagrange multi¬ 
pliers, one easily finds that p(x) must have the form 


2 k 

P ( x) = N exp [a^x + a 2 x + ... + a^x | . 

The constants a. . . . a, , N, can be found from the k + 1 inte- 
1 k 

gral conditions. 

The entropy of a probability distribution is a measure of 

(F-7) 

the uncertainty inherent in the distribution. Shannon 
showed that no other function can consistently bear this inter¬ 
pretation. Thus the principle of maximum entropy inference 
may fairly be described as the correct method for choosing a 

probability distribution when partial information about the 

(F-6) 

distribution exists. As Jaynes has pointed out, it is 

a quantitative embodiment of Laplace's principle of insuffi¬ 
cient reason. 

F. 7. Application to the Di st ribution of Energy Levels in a Con ¬ 
figuration 

We now have three parameters at our disposal — the mean 
energy W, the standard deviation of the energy o^, and the 
energy of the ground state, which we denote by W Q . Let 

x = (W - W Q )/e , 


(F-ll) 


(F -12) 


\ 


where c is an as yet unspecified unit of energy. The variable 
x is a dimensionless measure of the energy relative to the 
ground state. 


Applying the principle of maximum entropy inference, we 
find that the probability distribution has the. form 

p = N exp - b)“i] (x > 0) . 

The energy unit e has been chosen so as to make the coefficient 
2 

of x in the exponential -1. By integrating the last equation 
over the range of x we obtain a formula for the normalization 
constant N: 

N = 2tt” s [ 1 + erf (b) J" 1 
where erf(x) is the error function 

x X 2 

erf (x) = 2 ir~ 2 e" Y dy . 

0 

Thus p(x) depends on the single parameter b. 

Let 


. - (W - w Q )/e , a « aje j 

U and o are the mean energy and the standard deviation in 
dimensionless units. 

The integral conditions which ensure that the probability 
distribution has the required mean and variance yield two inde¬ 
pendent relations from which b and e can be evaluated in terms 
of u and a. The following form of these rSlatioi>s is conveni¬ 
ent for calculation: 

M- - b + |p(0) , 

(o/ur - u •• hp(o) ]/2u 2 . 


(F-13) 


(P-14) 

( F-15) 


(F-16ab) 


(F-17) 
(F-18) 
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By substituting the first equation into the second, we obtain a trans¬ 
cendental equation for b in terms of the ratio (exAO. The numerical 
solution of this equation for a range of values of oVu can be obtained 
easily. Table F-l gives or/u for various values of b, and the seme is 
plotted in Figure F-7. Having obtained the value of b from this 
table, one can evaluate the energy unit e from equations (F-17) and 
(F-16a). 

The agreement between the predicted and observed distributions 
seems to be satisfactory. If future studies should disclose a 
systematic discrepancy between the predicted and observed distribu¬ 
tion, then the very existence of such discrepancies would represent 
fresh information about the probability distribution and could be 
incorporated into the conditions of the problem. Moreover, it is 
likely that any such systematic discrepancy would have an assignable 
physical cause. Thus the principle of maximum entropy inference 
appears to be a powerful method for isolating those properties of com¬ 
plex spectra that are susceptible to interpretation in reasonably 
simple terms. 

gtitiitisil StmlBilan at tht flafelaafctea gatsm 

Consider a transition array, defined as the aggregate of lines 
resulting from all possible transitions connecting two 





TABLE F-l 


Z 
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-- 

Jd... 

0-0 

.75122 

C.l 

.70515 

0.25 

.63495 

0.5 

.59937 

0.75 

.57272 

1.0 

.53512 

1.25 

.48832 

1.5 

.43851 

1.75 

.39107 

2.0 

.34742 

2.5 

.28235 

3.0 

.23567 

3.5 

.20203 

4.0 

.17677 
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subject to the k + 1 auxiliary conditions expressed by the 
above integral equations. Using the method of Lagrange multi¬ 
pliers, one easily finds that p(x) must have the form 

2 k 

p(x) = N exp fa^x + a 2 x + ... + a^x v ] . 

The constants a. ... a, , N, can be found from the k + 1 inte- 
1 k 

gral conditions. 

The entropy of a probability distribution is a measure of 

(F-7 ) 

the uncertainty inherent in the distribution. Shannon 
showed that no other function can consistently bear this inter¬ 
pretation. Thus the principle of maximum entropy inference 
may fairly be described as the correct method for choosing a 
probability distribution when partial information about the 
distribution exists. As Jaynes 1 has pointed out, it is 

a quantitative embodiment of Laplace's principle of insuffi¬ 
cient reason. 

F. 7. Application to the Distribution of Energy Levels in a Con ¬ 
figuration 

We now have three parameters at our disposal — the mean 
energy W, the standard deviation of the energy o w , and the 
energy of the ground state, which we denote by Let 

x » (W - W 0 )/e , 

where e is an as yet unspecified unit of energy. The variable 
x is a dimensionless measure of the energy relative to the 
ground state. 
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(F-ll) 


(F -12) 



0ICB43- ?T0 
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Figure F-7. Solution to Equations F-17 and F-18. 



given configurations. What can we say about the probability 
distribution of line frequencies in such an array, given the 
probability distributions for the energy levels of the two 
configurations? 


Let denote the energy of the ith level in the configu¬ 
ration a. The frequency corresponding to a transition ai * bj 
is given by Bohr's frequency condition, 

hv(ai,bj) -W al -»W bj . 

We assume for the sake of definiteness that configuration a 
lies above configuration b, so that the frequencies defined 
by this equation are always positive. 

If departures from LS coupling are not too great, the 
combination spectrum [the spectrum of frequencies v(ai,bj)) 
will have a clumpy structure. Transitions between levels 
belonging to the same pair of spectroscopic terms give rise 
to multiplet8 . and in LS coupling "the separation between 
multiplets is large compared with the separation between 
lines in a multiplet. 

The distribution of lines in a multiplet and the dis-. 
tribution of multiplets in a transition array present distint 
problems. In a first approximation we may assume that the 
two distributions are stochastically independent — that is, 
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(F-19) 



that the form of the probability distribution for lines in a 
multiplet does not depend on the position of the multiplet in 
its transition array. Similar remarks apply to the case of 
pure jj coupling. 

On the other hand, if the coupling scheme is not close to 
either LS or jj coupling, neither the individual energy spectra 
nor the combination spectrum will exhibit any marked fine struc¬ 
ture. in this respect the case of intermediate.coupling is 
easier to deal with than either extreme LS coupling or extreme 
jj coupling. 

There is a second important respect in which the stochastic 
description of the combination spectrum is simpler for inter¬ 
mediate coupling than for LS coupling. In LS coupling, transi¬ 
tions connecting levels belonging to terms of different multi¬ 
plicity (different S) are strictly forbidden. This well-known 
selection' rule substantially reduces the number of lines ,in a 
transition array that occur with nonvanishing intensity. More¬ 
over, it forces one to deal separately with the various classes 
of transitions corresponding to the various possible values of S. 
In jj coupling transitions connecting terms differing by more 
than one quantum number j are forbidden. In intermediate coup¬ 
ling, on the other hand, intersystem transitions are not only 
allowed, but often have high probability, since intermediate 
coupling prevails in heavy atoms, we may reasonably assume that 
the choice of an initial level and the choice of a final level 
are stochastically independent events. The remainder of this 
discussion proceeds from this assumption. 



T 


\ 


If we regard the quantities v, as random variables 

[whose possible values are the numbers v(ai,bj),W ai , W^, re¬ 
spectively] , then the assumption just stated means that the 
random variables V» a and are stochastically independent, so 
that the random variable v is the difference between two stoch¬ 
astically independent random variables. This fact enables us 
to write down the probability distribution of v in terms of the 
probability distributions of and According to a well- 

known theorem in probability theory, the first distribution is 
just the convolution of the probability distributions of W a and 
W^. Let P a (x) denote the probability distribution of and p^ 
denote the probability distribution of Then the probability 

distribution p ab of 


(F -20) 


is given by 


P ab (y) " / p a (x)p b (y “ x)dx • (F “ 21) 

From this equation, or, more directly, from the definition of 
stochastic independence, one can show that the expectation and 
the variance of W ab are given by 


<*- 22 ) 


Var W . - Var + Var W. . 

SD ft D 


(F -23) 


As an illustration of these rules we consider the example 
of two normal distributions. By an appropriate choice of sero- 
points we can always ensure that * a ■ ■ l# ab « 0, so that we 
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may put 


p (x) * (27^)’^ exp [-x 2 /2a 2 ] 


P b (x) = (2t ro 2 )" 2 

Since the normal probability distribution is symmetric the 
probability distribution p fe of -W^ coincides with p^ in the 
present example? but this is not true in general. A short 
calculation, in which one performs the integration by com¬ 
pleting the square in the exponent, gives 

p (x) = ( 27 tct 2 )" 2 exp [-x 2 /2o 2 ] , 


4 




the probability distribution of the difference W ab is also 
normal and its variance is the sum of the variances of and 
W^. This is the well-known "reproductive property" of the 
normal distribution. 

We assume that the probability distributions p # , p^ have 
the form derived earlier in the present report. It is conveni¬ 
ent to introduce the function 

1 (x > 0) 

e(x) - 

0 (x < 0) 

With the help of this function we can write the probability 
distribution for If in the form 


P_(x) « <p(x>o ,B )€(x) , 


(F-24) 

(F —2 5) 


(F-26) 

(F-27) 


( F “28) 

(F-29) 
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where 


<p(x;o a £ a ) = N a exp [-(x - &^) 2 /2a*] . ( 

The definition of is exactly analogous. Note that 

P a (x) - P a (-x) ( 

We therefore obtain for the probability distribution of W afa 


p ab (y) " J - y>o b iP b )e(x)e(x - y)dx 


Prom the definition of the function €(x) it follows at once, that 


€ (x) e (x - y) 


e(x - y) (y > 0) 


e(x) (y < 0) 


p ab (y) " J " yj0 b ,& b )dx (y < 0) ; 

0 

00 

p ab ( y ) » /*<x,o a .P a )*(x - yjOjj^hJdx 

y 

00 

J <p(y + z;c a ,P a )<P(*f 0 b ,P b ) d* (y > 0) . 


The integrations are straightforward (ms the Mathematical Mote at 
the end of this appendix). The final result is: 


F—30) 

P-31) 

F-32) 

P-33) 

P-34a) 

P-34b) 
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Pab<y> - P ab<*> 


( y < 0) 


/Vb 2 + ®b a * 2 + y (T * 2 ' 1 

1 - erf f -*-**- b * * ■ 


(F-35a) 


■ P. b (y) 


where P ab <y) « \ 


b g a + y g a \ 

< 

/^ 2 + v , 2 - y V \1 

\ a a a b V2 (o- a 2 + <T b 2 ) yfj (y > 0) 
(F-35b) 

-. 2 - b 2 J * b L ^(-. 2 - b 2 ) 


and erf(x) has been defined in equation (F-15) 


F.9. Ground State Energy for a Configuration of Equivalent Electrons 

He consider a configuration of equivalent electrons 4 q . The energy 
of the completely filled shells is common to all the terms of the con* 
figuration, i q . It can be calculated, if desired, but is of no conse* 
quence for our purpose. Furthermore, the average energy of such a 

(F*2) 

configuration can be calculated in a straightforward manner (Slater ). 

W. v (i,) * W av (i2) (P ' 36) 

24 

and » av (i 2 ) -J] <*.*) ?k <ni, n4) (F-37) 

k 

|l 

where the coefficients f for configurations of p, d, f and g electrons 
u 

are given in Table 3*1 and F are the Slater F integrals defined in 
Equations (B-26) and(B-27). 
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Us shall describe below a method of calculating the energy of the 
ground term. 

It is well known that a subshell l can contain at most 2(2i + 1) electrons 
and is therefore exactly half filled when q • 21 + 1. It is also well 
known that the term energies (with respect to the average energy) of a 
configuration with n electrons are the same as those of a configuration 
with [2(2i + 1) - q] electrons. It is, therefore, necessary only to 
consider configurations which are half or less than half filled, i.e. 

< 21 + 1 . 

The ground term of any configuration has the largest value of S 
for that configuration. If several terms have the same largest value of 
S, say S M|t , it is the one with the largest value of L, say L^, 
among them. Each term is highly degenerate, but one of the SL Mg 
states belonging to the term (S |Mv , L mv ) Is such that 

"s- s -x 

"l * L «u.. 

The zero-order state, which belongs uniquely to the above set of values 
of the quantum numbers can be found easily. Mg * S Mav arises when all 
spins are parallel and in this case m # of all n electrons is +| (.r -i ) 

Therefore, 

Mg - - \ q (F-38) 

and the multiplicity of the ground term is (2 S mv + 1) or (q + 1). 
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The largest value of consistent with the above requirement on the 
nig and Mg Is obtained by a simple application of the exclusion principle. 
The in^ values of the q equivalent electrons must be 

1, . (i - q f 1) (P-39) 

For this case 

\ - f - 1 * » < p -«» 

The corresponding zero order state Is designated in the usual way by 

(i) + u- D + . a - q + 1) + 


where the quantity in parenthesis is the value of m^, of individual 
electrons and the + sign indicates that m g - t- 1/2. Hie ground terms of 
a few configuration of type 2° are given in Table F-2. For a shell which 
is exactly half filled (q - 22 + 1), it is readily seen from Equation 

(F-5) that L n>ax « 0 and the term is + ^S. 

TABLE F-2 

Ground Terra for Configurations of Equivalent Electrons 

Configuration Multiplicity L Cround Term 

1= q + 1] U | (2/-q+l)J 


P 

d 

f 

8 


2 

3 

7 

5 


3 

4 
8 
6 


3 

0 

10 


P 

'*F 

«S 

6 ki 
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The energy of the ground term in the absence of external fields is 
the same as that of any one of the zero order states of this term, say 

of the zero order state mentioned above. This energy can be found by 
(F-3) 

standard methods (Condon and Shortley ), and is given by, 

W mV V I a k (/ mj, / m^) - b 1 (i mj, / m^)| F k (ni, ni) 
a>& k* *• ■* 

(F-41) 

where the first summation is over all two electrons combinations and 
the second summation is over k which takes on all even values from zero 
to 21. The coefficients a k and b k are defined by Equations 8 6 (14) and 

8* (6) of Condon and Shortley, and are tabulated by them in Tables 
6 6 

1 and 2 for s, p, d and f electrons. 

^ As an example we calculate below the energy of the ground term 

| of configuration f 7 . The zero order state of lowest energy is 

3 + 2 + 1 + 0 + -1 + -2 + -3 + 

and uniquely belongs to the 8 S term. Table F-3 gives the appropriate 
k k 

e and b for f electrons. From the susm of these quantities, given in 
the last row of this table, it readily follows that the energy of the 
ground term, *8, of the configuration f 7 is 

I 
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TABLE P-3 


* k *nd b l ‘ (/ mjj, t m^) for f-electrons 

Note: *° ■ 1 for *11 cases 


b° * 0 If mjj 4 


Electrons 

a p 

n 

* 4 

x 1089 

Bfl 

; b 2 

i 

x 225 

b* 
x 1089 

b 6 

x 7361.64 

3 + 2* 

0 

-21 

-6 

25 

30 

7 

1 + 

-15 

3 

15 

10 

54 

28 

0 + 

-20 

18 

-20 

0 

63 

84 

-1 + 

-15 

3 

15 

0 

42 

210 

-2 + 

0 

-21 ‘ 

-6 

0 

0 

462 

-3 + 

25 

9 

1 

0 

0 

924 

2 + 1 + 

mm 

n 

-90 

15 

32 

105 

0 + 



120 

20 

3 

224 

-1 + 

■fl 

H 

-90 

0 

14 

378 

-2 + 


i: 

36 

0 

70 

504 

-3 + 

H 

1 

-6 

0 

0 

462 

l + 0 + 

12 

6 

-300 

2 

15 

350 

-1 + 

9 

1 

225 

24 

40 

420 

-2 + 

0 

-7 

-90 

0 

14 

378 

-3 + 

-15 

3 

15 

0 

42 

210 

0* -1 + 

12 

6 

-300 

2 

15 

350 

-2 + 

0 

-42 

120 

20 

3 

224 

-3* 

-20 

18 

-20 

0 

63 

84 

-1 + -2 + 

0 

-7 

-90 

15 

32 

105 

-3 + 

-15 

3 

15 

10 

54 

28 


mm 

-21 

-6 

25 

30 

7 

IS 

m 

-77 

-462 

168 

• 

616 

5544 


200 - 

























F. 10. Mathematical Mote 


From Equations (F-30) and F-34a) we obtain, for y < 0, 


00 

r 

(X - P a ) 2 (X - y - P b ) 2 ■ 

• h . s j “p 

, 2 (J 2 2 a . 2 

o 

La b J 

00 

r i 

-"b J “» 

-(a * 2 + 2 qx + r) dx 

<■ J 

o 



dx <F-43a) 


(P-43b) 


<r 2 + a 2 

, _fl_fe_ 


(F-44a) 


■ P a a h 2 + (y + fi b ) cr a 2 

‘ -.V 

P „ 2 g b 2 ♦ (y ♦ ft b > 2 tr a 2 

1 -.V 


(F-44b) 


(F-44c) 


The Integral on the right hand aide of Equation (F-43b) can be evaluated 
by Mane of formula (6b), Table 314, page 65 of the coapllatlon of 
definite Integrals by GrSbner and Hofrelter. Applying this foraula 
we obtain Equation (F-35a). 


For y > 0, we combine Equations (F-30) and (F-34b) end get 



? f <* - <-y + 

p.b (1,> ■"b J • xp -,71 

o L a 


(2 - P b >' 


(F-45) 


which is Che same as Equation (F-42a) if we interchange the subscripts 
a and b and replace y by (~y). Performing the integration as before 
we obtain Equation (F-35b). 
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APPENDIX G 

H. E. Stubbs and G. A. Victor 

THE EFFECT OF LINE WIDTH, STRENGTH, AND SHAPE 
ON RADIATIVE TRANSFER IN AN EXPANDING GAS 

G.l Introduction 

The details of line width, shape, and strength have a very important 
role in the calculation of radiative transfer phenomena in non-expanding 
gas clouds. This dependence on details of line shape introduces very 
serious complications in most astrophysical applications^ 6 since 

reliable information on line shapes is available only for simple gases 
as hydrogen and helium^ 6 ” 3 ^. The object of this appendix is to investi¬ 
gate a simple model of a rapidly expanding gas which results in expressions 
for line absorption which are largely Independent of detail line shape. 

The origin of the simolification in this model is found in the fact that 
the Doppler ahift due to the velocity of expansion becomes the dominant 
factor, and that the details of line shape due to collision and thermal 
Doppler broadening become unimportant. 

Thia appendix will be divided into three sectiona. Part G.2 will 
be devoted to some simple estimations of line widths due to the thermal 
Doppler broadening, and to compare it with the Doppler shift due to the 
expansion velocity. Part G.3 will show that resonance absorption is 
independent of line shape. Part G.4 will present a simple calculation 
for an expending hydrogen gas for the transmission of the 82, 260 cm"* 
line. 
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G.2 Doppler Shift Due to Velocity of Expansion 


I 


We shell construct the following idealised model for our expending 
cloud: 


The number density of hydrogen atoms is given by: 

p (r,t) - 1.5 x 10 3 /t 3 particles/cm 3 r < R o (t) 

p (r,t) - 0 r > R o (t) 

R o ( t) - V # t 

V Q - 1 x 10’ 2 c (6-1) 

t is the time in seconds 
c is the speed of light - 3 x 10 l0 cm/sec. 


the relative velocity between two particles is along their line of 
centers and is given by: 

*»i -'-T- ( °- l > 

Where r^ 2 ia the perticles seperetion: 

The relativistic Doppler formula is: 


X » X o 0 cos 6 (0-3) 

8ince the velocity is along the line of centers we can set cos 0 ■ 1. 
Since V rc i « C we have the relation - .02. therefore: 

SlAX " ,02 \> 

The formula for the thermal Doppler width is 
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*X T - X 


(G-3) 



for Hydrogen we cen conetruct the following table: 


3/2 kT 

5 S 

V W 

100 ev 

3.2 x 10“ 4 X 

' 

1.6 x 10‘ 2 

10 ev 


5.0 x 10 -3 

1 ev 

3.2 x 10 -5 X 

0 

1.6 x i0’ 3 

.1 ev 

io"\ 

5 x 10’ 4 


o 



From Table G-l we see that the Doppler width due to thermal Motion is 
always much smaller than the maximum Doppler shift due to the velocity 
of expansion. The results of the rest of this appendix depend on this 
fact. 


G.3 Independence of Resonance Absorption on Line Shape 

Let a photon of frequency u be emitted at r ■ 0. Let the 
probability of emission of a photon of energy between ftu and If (w + dw) 
per unit time per particle be given by P(u). The probability that a given 
atom will absorb a photon of energy Vu is given by, 

p(w) - k P(u) (0-6) 

from the principle of detail balance. 
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The probebility chat a photon o£ frequency u will tranaverse a die- 


tance r before being absorbed Is given by: 

J A(|,w)d| 

I (u,r) * e ° 


(G-7) 


with 


A(|,u) - p(|) p (w'[r,u ]) 


(G-8) 


where p(() is the density of particles which can absorb this resonance 
radiation. I(u,r) can therefore be written as: 


f 


P(i) p ( 41 - 


(G-9) 


I(u,r) ■ e 
where we have used the Doppler formula: 

Au * (G-10) 

and evaluated this at a particular time t Q . 

We Intergrate this over all frequencies, weighted with P(w) to get the 
probability that any photon will go a distance r. 


T(r) 


•£ 


+oo 

du P(u) e 


• J, p 


which will be re-written as 


I(r) 


r - k Jf 

• du P(w) e 


(l) K41 - tH) 


(G-li) 


( 0 - 12 ) 


267 



If p(|) it • constant out to a distance r, which ia conaiatant with our 
cloud aodel, we have: 


T(r) 


r°° >f‘ 

■ I du P(u) » + ^ Jy 
J~OD 


cj-Aoj(r) 

P(l)d| 


(G-13) 


where: 1C' » kp Au(r) » - 


Finally 


W (l- £> 


{ * *L K'POi'Jdu 1 

P(u)du e Ju 

•oo 


Using the result of part G.2 we can calculate this for r • oo 
_ -oo 

<-+00 

T(oo) - / P(u)due* 


£ 


Jd£ P(u')du' 


if wa let t) •/ Pt-j'Jdu' 

then dt) - P(u)du 


(C-14) 


(G-1S) 


T(oo) 




(G-16) 


Proa Equation (G-16) we see that for an infinite cloud 
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resonance absorption is independent of line shape. For a finite cloud, 
it is only the photons thatsre emitted near the edge that will show resonance 
absorption that is dependent on the form of P(u). The ratio of the thick¬ 
ness of this edge to the cloud radius is essentially We there¬ 

fore conclude that resonance absorption is controlled primarily by line 
strength, for this simple model of an expanding cloud. The object of 
this appendix is to investigate a simple model of a rapidly expanding gas 
which results in expressions for line absorption which are largely 
independent of detail line shape. 

Arguments similar to those .given above for resonance absorption can 
be given for line absorption due to two close, but non-overlapping lines. 

It is easy to see that an absorption line centered at can only absorb 
radiation from lines centered at such that: 

\ - h < 

It is again obvious that this absorption is independent of line 
shape except for lines where the equality holds in Equation (G-17). The 
effect of a finite cloud on the absorption of non-overlapping line radia- . 
tlon does little but introduce some simple geometrical factors. 

G.4 Iransmlfslon of, Hydrogen 82,260 cm" 1 feinj 

The line of interest is the 82,260 cm* 1 line which is emitted by a 
transition from the first excited state to the ground state. We will 
assume that the gas has a temperature history as shown in Figure 0-1. 

Xn addition, we assume that the gas is at thermal equilibrium at all times. 
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TEMPERATURE (•*) 


010 * 40-1390 



Figure G-l. Teaperature During Expansion. 
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At equilibrium, the population of states is given by the Saha equation 


N n “ N p N e n 2 h 3 [2™kT]‘ 3/2 exp (G-18) 

Up ■ number density of protons 

N e ■ number density of electrons 

* r ■ -onisation potential ( r degree) 

n - principle quantum number 

N r - number density of atoms in state n 

he number density of particles which can absorb the 82,260 cm -1 line is 
given by 

p(t,T) - p o (t)f(t,T) (G-19) 

Njd) 

«•>«. £(t.T) - -j c y 

Istimations of f(t,T) from Equation(C-18)lead to the results shown in the 
following table. 


Time 

fiams,its] 

°l c » 3 J 

R [cm] 

K t 

f t 

r 

CD 

IQ' 6 

1.5 x 10 21 

2 x 10 2 

7 x 10 7 f 

5.7 x 10 4 

2.5 x 10' 5 

10* 5 

1.5 x 10 18 

2 x 10 3 

7 x 10 5 f 

1.08x 10* J 

1.3 x 10' 1 

10* 4 

1.5 x 10 15 

2 x 10 4 

7 x 10 3 f 

.26 

1.4 x 10’ 4 

w- J 

1.5 it 10 33 

2 x 10 5 

7 x 10 l f 

1 

1.4 x 10* 2 

10* 2 

1.5 x 10 9 

2 x 10 6 

7 x 10“ l f 

1 





From Equation (G-16) we have: 


T t <»> ■ 7: 

t 

where the aubscript indicatea time dependence, 
ing expreaalon: 


(G-20) 

la given by the follow* 


R , . p „ £ (t.Tj.A.t (G . 

C 8n v 3 V 8« v 3 

o 

A la the ginatein coefficient #*'5 x 10 8 sec" 1 . The reaulta for K* 
et verioua timea ia shown in Table G-2. FiguresG-2 end G-3 ahow the 
varietion of T(ao) with time for thia simple model. 
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FRACTION TRANSMITTED TO 00-(T) 



Figure G-2. Tranealeeion in Expanding Medium. 


1 
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